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Abstract Low energy cosmic rays (up to the GeV energy domain) play a cru- 
cial role in the physics and chemistry of the densest phase of the interstellar 
medium. Unlike interstellar ionising radiation, they can penetrate large col- 
umn densities of gas, and reach molecular cloud cores. By maintaining there a 
small but not negligible gas ionisation fraction, they dictate the coupling be- 
tween the plasma and the magnetic field, which in turn affects the dynamical 
evolution of clouds and impacts on the process of star and planet formation. 
The cosmic-ray ionisation of molecular hydrogen in interstellar clouds also 
drives the rich interstellar chemistry revealed by observations of spectral lines 
in a broad region of the electromagnetic spectrum, spanning from the sub- 
millimetre to the visual band. Some recent developments in various branches 
of astrophysics provide us with an unprecedented view on low energy cosmic 
rays. Accurate measurements and constraints on the intensity of such par- 
ticles are now available both for the very local interstellar medium and for 
distant interstellar clouds. The interpretation of these recent data is currently 
debated, and the emerging picture calls for a reassessment of the scenario in- 
voked to describe the origin and/or the transport of low energy cosmic rays in 
the Galaxy. 
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1 Introduction 


The formation of stars is a central question in astrophysics (Shu et al. 
2019) 
While it is certain that star formation takes place inside interstellar molecular 
clouds (MCs) as the result of the gravitational collapse of their dense cores, 
the details of such process remain quite uncertain. 


MCs are cold, dense, magnetised, and turbulent (Crutcher} |2012 
nebelle and Falgarone} |2012} |Heyer and Dame} |2015). While the first two 


features in the list tend to favour the formation of stars, as they imply low 
thermal pressure support and short gravitational free-fall time, respectively, 
the latter two oppose to it, as they provide a non-thermal pressure support 
against gravity. The relevance of the magnetic pressure support and the impact 
of magnetohydrodynamical turbulence on the dynamical evolution of clouds 
depend on how tightly the gas and the magnetic field are coupled. 

Ultimately, the level of coupling depends on the gas ionisation fraction: a 
neutral gas would not be affected at all by the presence of a magnetic field 
and, conversely, a magnetic field would be frozen into a fully ionised gas. The 
ionisation fraction found in dense MCs is at the level of ~ 1077, which is quite 
small, but nevertheless much larger than what one would expect in a very 
cold (~ 10 K) gas which, due to its large cloud column density, is protected 
by external sources of ionising radiation 1989). It follows that an 
additional source of ionisation must be present inside MCs. 
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A similar line of reasoning can be pushed forward also in connection with 
the formation of planetary systems, which is intimately connected to star for- 
mation. This is because the dynamics and evolution of protoplanetary disks 
depend, again, on the gas ionization level, which dictates the effectiveness of 
mechanisms of magnetic field transport such as magneto-rotational instability 
or magneto-centrifugally launched winds (2011). 

Finally, a surprisingly complex chemistry has been revealed by spectro- 
scopic observations of interstellar clouds from the submillimetre to the visible 
band. Such chemistry is made possible by the large gas column density of 
clouds, which protects molecules from interstellar radiation and allow them to 
survive and proliferate. However, the rate of neutral-neutral chemical reactions 
is way too slow in the diluted and cold interstellar medium (ISM), implying 
that some level of ionization of the gas is mandatory in order to allow faster 
ion-neutral reactions to build up the molecules we observe 
2013). Likewise, the chemical rich- 
ness of planetary systems is inherited from the processes taking place in the 
parent MC, and is further affected by the ionization state of protoplanetary 
disks p20). 

What said so far clearly indicates that a central question in astrophysics 
is: what keeps the densest phase of the ISM slightly ionised? 

Once stars are formed from the gravitational collapse of MC cores, nuclear 
fusion reactions begin to operate in their hot interiors, and to create elements 
heavier than hydrogen in a process called stellar nucleosynthesis (Burbidge 
1957] (Cameron) [1957] [Trimble}[T994} [Wallerstein et al} 1997) Woosley 
et al.||2002). The very early Universe emerging from Big Bang ae 
was only made of hydrogen, some helium, and a very small fraction of lithium 
(Coc and Vangioni| (2017). It is thanks to stellar nucleosynthesis that the great 
variety of elements that we observe today was generated. Notwithstanding 
the great success of both Big Bang and stellar nucleosynthesis theories in 
predicting the observed abundances of elements, it became soon clear that a 
substantial fraction of the lithium and the totality of beryllium and boron 
found in the present day Universe must have been produced in another way. 
A peculiar origin of lithium, beryllium and boron (LiBeB) was suggested by 
the fact that their cosmic (solar system) abundances are much smaller (many 
orders of magnitude) than those of their neighbours in the periodic table (see 
Fig. g. [Burbidge et_al.][1957] pointed out that LiBeB are fragile, and are 
destroyed very effectively in the thermonuclear reactions taking place in the 
hot and dense stellar interiors. Their origin, then, had to be searched elsewhere. 
In the absence of a viable mechanism, they invoked an unspecified x-process, 
operating in an environment where both temperature and density must be 
low, as the responsible for the synthesis of light elements. 

A second fundamental questions then emerges: what produces the light el- 
ements (LiBeB) observed in the Universe? 

Remarkably, both the crucial questions asked above have the very same 
answer, that is, the interactions of low energy cosmic rays (LECRs) with in- 
terstellar matter. 
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Fig. 1 Cosmic abundances of elements in the solar system (black) and in cosmic rays 


(green). Figure from |Tatischeff and Gabici] (2018). 


Cosmic rays (CRs) are mostly made of energetic atomic nuclei (mainly 
protons). They fill the entire Galaxy and reach the Earth as an isotropic flux 
of particles. Unlike ultraviolet radiation, the intensity of CRs is not com- 
pletely attenuated by large column densities of matter. Therefore, CRs are 
the only ionising agents able to penetrate the depths of MCs and maintain 
there the gas ionisation fraction at the observed level of ~ 1077. Moreover, 
secondary electrons produced in ionisation events cool by transferring their 
energy to the ambient gas, providing the source of heating needed to explain 
the gas temperatures measured in clouds. Not all CRs contribute significantly 
ooo o o the n o ae ae of the gas, but ie ae rae (< ; Tae ones 


eo link ee e and the aes 5 eee — ae ev- 
ident after the intensity of CR nuclei heavier than hydrogen was measured, 
revealing a striking difference with respect to cosmic abundances. The ratio 
between the abundance of CNO and LiBeB nuclei is very large (~ 10°) in the 
solar system but only of the order of a few for CR particles (see Fig. (i). This 
puzzle was solved in the early seventies, when it was understood that the syn- 
thesis of LiBeB mainly results from spallation of CNO nuclei of interstellar gas 
by LECRs. Spallation occurs when a nucleus in the interstellar gas (for exam- 
ple carbon), struck by a CR particle (for example a proton), ejects a number 
of lighter particles and transforms into a LiBeB isotope. To this direct channel 
one should also add the contribution from the reverse process, i.e. spallation 
of CNO LECRs by interstellar nuclei, and from a+ qa reactions, where two 
helium nuclei collide to synthesize lithium. These three main reactions involv- 
ing LECRs and interstellar matter account for most of the LiBeB isotopes 
produced by the x-process, which is now called spallogenic nucleosynthesis 


(Reeves et al} [1970 1971). 
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LECRs also impact on the diffuse interstellar gas on large scales. This fol- 
lows from the fact that the interstellar energy densities of CRs, magnetic field 
and thermal and turbulent gas are in rough equipartition, with LECRs provid- 
ing a sizeable contribution to the total energy density of CRs. Therefore, CRs 
may affect the dynamics and the characteristics of diffuse interstellar matter in 
various ea es as oy the pressure support needed to launch galactic 
winds or exciting jae ee turbulence as they 
stream across Be ea plasmas (Wentzel Wentzel} {1974} [Zweibel] (2017). 

Finally, understanding where and how LECRs are accelerated, how they 
are transported from their sources to the Earth, and what is their final fate 
is important per se: revealing the origin of CRs is one of the main open ques- 
tions in high energy astrophysics. Supernova remnant shocks, that form in the 
ISM as the result of supernova explosions, are often invoked as the sites where 
the bulk of galactic CRs are accelerated 2013). However, 
two things should be kept in mind. First of all, the supernova remnant origin 
remains to date an hypothesis (see for a critical review), 
and alternative acceleration sites have been proposed, including stellar wind 
termination shocks, stellar clusters, or superbubbles, i.e. the cavities inflated 
in the ISM by the collective effect of stellar winds and recurrent supernova 


explosions in star clusters So and Montmerle| |1983} [1983] [Bykov] [2014] [Lin-| 
[genfelter| [2018} 2019). Second, at present it is not quite 


clear whether LECRs, which are g acterized by sub-GeV particle energies, 
have the same origin as higher energy ones (see for example the discussion 
on protostellar shocks as LECR accelerators in [Padovani et al.| 2020), nor if 
the intensity of LECRs measured within the solar system is representative of 
the entire Galaxy or simply reflects some very local ambient conditions (for 
example the presence or absence of nearby CR sources, see 
or the fact that we live inside an interstellar cavity called the local bubble, see 
e.g. [Silsbee and Ivlev|2019). All in all, the origin of galactic CRs is still not 
well understood, and this is particularly true for LECRs. 

Addressing all the issues raised above goes beyond the scope of this review, 
which will be focussed on the impact that LECRs have on the densest phase 
of the ISM. This choice is motivated by a number of recent developments in 
pee exploration (the Voyager probes crossing the heliopause, 
3] [2019), laboratory astrophysics (the measurement of the dissociative re- 
ate Pa 5 Ae , a pivotal molecule in interstellar cloud chemistry, 


McCall et al. and CR astrophysics (the advent of a precision era in 


ao energy fie L einen from space, [Boezio et _al.|/2020} [Aguilar et al] 


. These activities triggered a renewed interest in both observational and 
rot otical studies of interstellar clouds irradiated by LECRs. An attempt to 
review this field seems therefore timely. 

Before proceeding, I provide here a list, certainly incomplete, of review 
articles covering aspects of the physics of LECRs which are not treated here. 
The reader interested in the CoA nucleosythesis of light elements is 


referred to the reviews by |Reeves| ), (Vangioni-Flam et al. (2000), and 
Tatischeff and Gabicij (2018). The oe a the isotopic composition of CRs 
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is cose in FRE aed Soil CHG [Wink a] OT) 
The extended review by |P E treats, among others, topics 
which are not covered by the present review, including the impact of LECRs 
in circumstellar disks, stellar cosmic rays, LECR acceleration in protostellar 
jets, and extragalactic studies of LECRs. Finally, a discussion on the effect 
that CRs might have on the origin of life can be found in [Dartnell| (2011). 

The remaining of the paper is structured as follows. After an introduction 
on CRs (Sec. B), we will review the difficulties encountered by direct (Sec. 
and indirect (Sec. [4]and|5) attempts to measure the intensity and spectrum of 
LECRs. In these two Sections, we will also highlight the recent observational 
breakthroughs that radically changed or improved our knowledge of LECRs. 
Sec. [6] will be devoted to a discussion on the transport of LECRs in and 
around MCs. To this purpose, we will consider both isolated MCs and clouds 
located in the proximity of powerful CR sources. A list of open issues in LECR 
astrophysics will be provided in Sec. {7| and we will conclude in Sec. [8| where 
future perspectives will be also outlined. 


2 What are cosmic rays? 


CRs are energetic charged particles that reach the Earth’s atmosphere from 
outer space. The distribution in the sky of the arrival direction of CRs is 
remarkably isotropic. The vast majority of CR particles are atomic nuclei, 
with a contribution from electrons at the percent level, and an even smaller 
one from antimatter (positrons and antiprotons). Amongst CR nuclei, protons 
largely dominate, with helium and heavier nuclei (metals) amounting to ~10% 
and ~1% of the total number of particles, respectively (for reviews see |Gaisser| 
020). 

With the exception of anomalous CR¢"| they are originated outside of the 
heliosphere. The vast majority of them are in fact originated within our Galaxy, 
and only the highest energy particles, that won’t be discussed here, likely 
have an extragalactic origin as they can hardly be confined by the Galactic 
magnetic field. Understanding the origin of the bulk of CRs is one of the most 
fundamental and long standing questions in high energy astrophysics. The 
CR elemental composition and energy spectra are now measured with great 
accuracy and carry crucial information about where CRs are accelerated and 
how they traveled from their sources to us. 

Fig. [2| shows the local spectra of CR nuclei of atomic number from Z = 1 
to 8 (top to bottom). Data have been collected by detectors operating at the 
top of the Earth’s atmosphere: the satellite borne PAMELA (blue) and the 
Alpha Magnetic Se station (Kikiant aa (AMS- 02, = shaded regions) mounted on the 


International Space Station [Adriani et al.| {2014} [AMS Collaboration et. al. 2014} |AMS Collaboration et al. 


1 Anomalous CRs are neutral atoms in the local ISM that, due to the motion of the Sun, 
enter the heliosphere, become partially (mainly singly) ionised due to charge-exchange, solar 
radiation, or electron impact, and are eventually accelerated in the solar wind up to energies 
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Fig. 2 Spectra of CR nuclei of atomic number from Z = 1 to 8 (top to bottom) as a function 


of rigidity. Data from AMS-02 (red) and PAMELA (blue). Figure from (2019). 


(2002). For a better representation, spectra have been multiplied by the par- 
ticle rigidity] R to the power 2.7. The suppression observed in all spectra at 
rigidities smaller than a few tens of GV is due to the effect of the solar wind 
on low energy particles, and will be discussed in detail in Sec. 


Two things should be noted. First of all, for particle rigidities above ~ 10 
GV, the spectra of lithium, berillium, and boron are markedly steeper than 
those of the other elements, which are roughly flat once multiplied by R27. 
Second, while the Solar abundances of LiBeB are many orders of magnitude 
smaller than that of carbon, such a difference is way smaller for CRs (see 
Fig. (ip. As seen in the Introduction, this peculiarity in the abundances of light 
elements fits well within a scenario where such elements are produced as the 
result of the interactions of LECRs in the ISM. 


2 The rigidity of a fully ionised nucleus of momentum p is R = pc/Ze, where e is the 
elementary charge. If pc is expressed in eV and the particle charge in natural units (e = 1) 
the rigidity has units of Volts. Particles of equal rigidity have the same gyration radius 
around a magnetic field of strength B, rg = R/B, and therefore follow the same trajectory. 
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2.1 The confinement of cosmic rays in the Galaxy: boron-over-carbon ratio 


To understand more quantitatively the peculiarities in the intensity and spec- 
tra of light elements, consider the production of CR boron due to spallation 
of heavier CR nuclei by interstellar matter. In such process, an energetic nu- 
cleus (most likely of carbon or oxygen, which are the most abundant elements 
heavier than boron, see Fig. [1) hits a nucleus of interstellar gas (most likely an 
hydrogen nucleus), loses one or few nucleons in the impact and transforms into 
boron. Let us call no( E) and no(£) the average number density of CR carbon 
and oxygen nuclei in the Galaxy having an energy per nucleon E£, and ocg 
and o9_,g8 the appropriate spallation cross sections for boron production. In 
first approximation, the energy per nucleon is conserved in spallation reac- 
tions, and at large enough energies the spallation cross sections are roughly 
energy independent (see e.g. Fig. 4 in (2018). There- 
fore, for a fixed value of E, the production rate of CR boron in the Galaxy 
is gp ~ n{$Mu(oc Bnc + Go>BNo), where niP™ is the hydrogen density of 
the interstellar gas and v is the velocity of the incident nucleug?| Noting from 
Fig. [2| that the measured abundances and spectral slopes of CR C and O are 
almost identical, no(E) ~ no(£), the expression for the production rate of 
CR boron can be simplified to qg ~ ni?“ unc(oc-.B + cog). By doing so, 
we are implicitly assuming that the spectrum of CRs observed locally is rep- 
resentative of the entire Galaxy, i.e. that there are no large spatial variations 
of the CR intensity in the Galactic disk. 

After being produced, CR boron nuclei spend some time Ty5)4 wandering in 
the ISM before escaping the Galaxy, or they spallate on interstellar hydrogen 
to transform into lighter nuclei in a typical time tg = (n}?“ogv)—!. Here, 
og is the total spallation cross section for CR boron. The balance between 
production and escape/destruction leads to an equilibrium density of CR boron 
equal to ng = qBTeff, Where Toft =Tigu + Tp” is an effective timescale. 

To proceed further, it is convenient to introduce a quantity called gram- 
mage, defined as the mean amount of matter traversed by CRs in a given time 
Tj, aS Xj = mpn Mvyr;, where mp is the hydrogen (proton) mass, and where 
the presence of helium and heavier elements in the ISM is neglected. Then, 
Xzsm(£) represents the mean amount of matter traversed by CRs before es- 
caping the Galaxy, while Xg corresponds to that traversed by CR boron before 
spallating to transform in a lighter element. Note that while the former may 
be an energy dependent quantity, the latter is not, as long as particles of large 
enough energy are considered (as og is roughly constant). Rearranging all the 
equations derived above one finally gets an expression for the CR Boron-over- 


Carbon ratio (B/C), which is an observable quantity (Gaisser et al. |2016): 


NB Xism 0C0+B+90>B (1) 
no 1+ Aye Mp ` 


3 Nuclei characterised by the same energy per nucleon move at the same speed. 
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Fig. 3 Boron over carbon (B/C) ratio. Data from AMS-02 (blue points [Aguilar et al.|2016} 
and Voyager 1 (red points, [Cummings et al.|[2016). The dashed line represents the power 
law E~1/3, The grammage derived by data and the corresponding CR residence time in the 
ISM are shown in the inset (solid line) together with the power law E~1/2 (dashed line). 


Eq. [1] illustrates well how the measurements of the CR B/C ratio provide us 
with an estimate of the residence time of CRs in the ISM. In particular, the en- 
ergy dependence of the B/C ratio is fully determined by that of X;gu&7Trsm, 
as all the other quantities on the right hand side of the equation are roughly 
energy independent. 

The observed B/C ratio is shown in Fig. |3| as measured at the top of 
the Earth’s atmosphere by AMS-02 (blue points), and in the local ISM by 
the Voyager 1 probe (red points). The two datasets do not seem to connect 
smoothly, and exhibit quite different behaviours. The B/C ratio measured by 
AMS-02, after a plateau, decreases steadily with particle energy, while that 
measured by Voyager 1 is virtually energy independent. Here, we will limit the 
analysis to the higher energy dataset (AMS-02), which can be interpreted in 
a rather straightforward way, and we avoid a discussion on low energy data, 


whose interpretation is very uncertain (see e.g. Fig. 1 in|/Tatischeff et al./2021 
or Fig. 9 in [Cummings et al./2016). 


The first thing that can be inferred from the high energy observation of 
the B/C ratio is that the residence time of CRs in the ISM is indeed energy 
dependent, and it decreases with energy. The absolute value of the grammage 


Xrsm can be read directly from Fig. |3| This can be done by using Eq. |1} with 
the measured values oc_,p + 00-.B ~ 100 mb and Xp ~ 7 g/cm? (Gaisser 


2016), which gives Xrsm œ% 13 and ~ 7 g/cm? for a particle energy 
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of ~ 3 and 10 GeV/nucleon, respectively. Such an estimate is quite close to 
results from accurate studies (Jones et al.|2001] obtained X; ~ 12 g/cm? at ~5 
GeV/nucleon). The grammages derived here correspond to residence times for 
CRs in the ISM equal to Trsm ~ 8 and 4 Myr for particle energies of 3 and 10 
GeV/nucleon, respectively, for a typical density of the ISM ni?” ~ 1 cm73. 

The decrease of the observed B/C ratio with particle energy implies that 
well above 10 GeV/nucleon, Xzsm becomes much smaller than Xg and Eq. 
asymptotically reduces to ng/ncxXrsm.|Aguilar et al.| found that B/C 
data can be fitted with a single power law in energy E~’ with 6 ~ 1/3 for 
particle energies larger than ~ 30 GeV/nucleon. It is sometimes erroneously 
claimed that such measured value of 6 also describes the energy scaling of the 
grammage Xyg mMocE~, However, such claim is incorrect because the parti- 
cle energies probed by AMS-02 are not large enough to satisfy the condition 
Xrsu > Xp. This is shown in the inset of Fig. where the value of the 
grammage inferred from data and from Eq. |l}is plotted for particle energies 
exceeding 3 GeV/ nucleorf{] The decrease of the grammage with particle en- 
ergy is in fact faster than E~'/?, and rather closer to ~ E~'/?. This is in 
broad agreement with more sophisticated studies (the reader interested in an 
accurate derivation of the grammage is referred to the recent works by [Evoli] 
and [Génolini et al 2019) 

Now that the grammages traversed by CRs of different energies have been 
estimated from data, it is worth comparing them with the gas surface den- 
sity of the Galactic disk X4isk, which is of the order of few milligrams per 
squared centimeter [2001). The latter is a proxy for the minimum 
possible average grammage, i.e., that accumulated by a CR produced in the 
disk and moving away from the production site along a straight line. Strikingly, 
the accumulated grammages are about three orders of magnitude larger than 
the disk surface density. This remarkable difference can be interpreted in two 
ways: either CRs remain somehow (diffusively?) confined in the gaseous disk 
for a time Trg before escaping the Galaxy, or they are confined in a much 
larger volume for a time Tese > Trsm and during this time they cross the disk 
= Xrgm/Xaisk = O(10° — 10*) times. The problem is that the measurements 
of the CR, B/C ratio can only constrain the product between the typical con- 
finement time and the average gas density felt by CRs (i.e. the grammage) and 
not the two quantities separately. Therefore, in order to distinguish between 
the two scenarios, an additional observable is needed. As the mean density 


of the Galactic disk is very well measured (e.g. 2001), what is much 


needed is an independent measurement of the confinement time. 


2.2 Radioactive isotopes as cosmic ray clocks 


Short lived radioactive isotopes are present in the cosmic radiation, being 
produced in the spallation of CR nuclei by interstellar matter. They can be 


4 At smaller particle energies the effect of Solar modulation becomes significant. 
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used as cosmic ray clocks, as their abundance relative to stable isotopes may 
depend quite strongly on the escape time from the Galaxy Tese. This is true 
for isotopes whose decay time Traa does not exceed Tesc, otherwise they would 
behave exactly as stable isotopes. Also, Traq should be long enough to allow 
the intensity of radioactive isotopes to be detectable by available instruments. 
10Be, with a lifetime (at rest) of Trag ~ 2 Myr (half-life equal to Ty/2 ~ 1.4 
Myr) i is therefore an ideal cosmic ray clock (Ptuskin and Soutoul| (Ptuskin and Soutoul} {1998} 1998). In the 
observer rest frame, the lifetime of a !°Be CR isotope has to be multiplied by 
its Lorentz factor y due to relativistic time dilation. The measured abundance 
of !°Be with respect to its stable companion ?Be provides, to date, the most 
stringent constrain on the escape time of CRs from the Galaxy. 
Unfortunately, at least three complications arise when studying '°Be in 
CRs. First of all, measurements of such isotope are available only at parti- 
cle energies smaller than ~ 2 GeV/nucleon (e.g. |Nozzoli and Cernetti) [Nozzoli_ and Cernettil [2021] 
and references therein P} At these energies, ionisation losses in the ISM can- 
not be neglected, as it was done, implicitly, in writing Eq. [2] ( see Sec. 
for a discussion on energy losses). Second, besides spallation of CR C and 5 
by interstellar matter, other reactions produce a significant amount of !°Be 
(e.g. (2018). Predicting the abundance of such isotope in CRs 
therefore would require the solution of a network of reactions. Third, at low 
particle energies, where measurements are available, the particles’ Lorentz fac- 
tor is y S 2 and the lifetime of !°Be in the observer frame is a factor of a few 
shorter than Trsm, which is in turn smaller (or at most equal, in the limiting 
case when the confinement volume coincides with the Galactic disk) than Tee. 
Assuming that the confinement of CRs in the Galaxy is diffusive, and that 
the Galaxy is a flattened system of height H, one can estimate the particle 
diffusion coefficient as D ~ H?/Tesc. It follows that the CR isotopes of !°Be 
observed at the Earth must have been produced within a distance l ~ y DTrad, 
i.e. within a volume significantly smaller than the entire confining volume, as 
l~ H4/Trad/Tese < H. To interpret correctly the measurements of !°Be, then, 
it is mandatory to know with some accuracy the spatial distribution of matter 
in the local ISM. This is of particular relevance as the Solar system is known to 


be embedded in a large cavity in the ISM, called local bubble (see 
2001 2002; and the discussion in Sec. |7| 


Despite all these complications, the way in which the measurements of the 
abundance ratio !°Be/®Be in CRs can be used to constrain their confinement 
time in the Galaxy can be understood in a qualitative way if one considers 
only the most important difference between the behaviour of the two iso- 
topes °Be and !°Be: the former is stable, the latter is not 
(1958). Assume then that the CR isotopes !°Be and °Be are produced at a 
rate qio and q9 and removed from the Galaxy in a typical time Traa and Tese, 
respectively. In the absence of any other effect, at equilibrium one would get 
n(?°Be)/n(°Be) = (q10/99)(YTrad/Tesc)- The isotopic ratio at production q10/q9 


5 At the International Cosmic Ray Conference that took place (virtually) in Berlin in July 
2021, the AMS-02 collaboration presented preliminary results on the measurements of the 
isotopic ratio 1°Be/9Be up to energies of ~ 10 GeV/nucleon. 
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is known, as it depends mainly on the production cross sections. Therefore, 
the observed abundance ratio n(‘°Be)/n(°Be) provides an estimate for the CR 
escape time from the Galaxy Tego. 

The comparison of model predictions with the observed n(!°Be)/n(° Be) 
ratio indicates that the residence time of CRs in the Galaxy Tese is significantly 
longer than the time spent into the ISM 7;sy,, and is found to be of the order of 
~ 100 Myr for particle energies < 1 GeV/nucleon (e.g. |Ptuskin and Soutoul| 
1998). This implies that the confinement volume of CRs is larger than the 
Galactic disk (where grammage is accumulated), and that CRs spend most 
of the time outside of the disk, in a larger region called Galactic halo, before 
leaving the Galaxy. The existence of an extended Galactic halo is also required 


by radio observation of the Galaxy at high latitudes (Beuermann et al.||1985). 


2.3 Diffusive models for the transport of cosmic rays 


The transport of CRs in a disk-plus-halo system is often modelled in a simpli- 
fied way by taking the Galaxy as a box inside which CRs move at a speed v and 
undergo a random walk in space with mean free path A. In this scenario CRs 
diffuse with a diffusion coefficient D ~ Av, and occasionally cross an infinitely 
thin disk where interstellar matter is concentrated (e.g. 2001). If 
H is the height of the halo, it can be shown that a CR particle of velocity 
v will cross the disk Neross ~ Hvu/D ~ H/A » 1 times before leaving the 
Galaxy. The grammage determined from the B/C ratio is then ~ NerossXdisk 
and therefore constrains the ratio D/H only, and not D and H separately. On 
the other hand, the measurement of the abundance ratio '°Be/*Be is sensitive 
to the escape time only, which is Tese ~ H?/D. Then, the simultaneous knowl- 
edge of these two abundance ratios (B/C and !°Be/*Be) allows us to estimate 
both the diffusion coefficient (D ~ 107° cm?/s at ~ 1 GeV/nucleon) and the 
size of the Galactic halo (several kiloparsecs). For an accurate derivation of 
these two quantities see e.g. and references therein. 

Finally, the energy dependence of the diffusion coefficient can also be 
constrained within the framework of the simple diffusive scenario outlined 
above. Indeed, for relativistic particle energies (v ~ c) the scaling simplifies to 
Dol1/XrguCl/TeseXE® with 6 ~ 0.5, as determined from Fig. |3| Accurate 
studies provide a range of values 6 ~ 0.3...0.6 compatible with observations 
(Strong et al. |2007). In fact, it would be more appropriate to think in terms 
of particle rigidity R rather energy per nucleon E (see footnote 2), but for rel- 
ativistic particles one simply has R ~ (A/Z)E ~ 2E, and Dx R°. This shows 
that particles of larger rigidity diffuse faster and leave the Galaxy earlier than 
particles of lower rigidity, and this fact can be used to constrain the spectrum 
that CR sources, whatever their nature, must inject in the ISM. 

To do so, let us assume that CR sources located in the Galactic disk in- 
ject in the ISM primary nuclei (H, He, C, etc.) characterised by an energy 
spectrum qp inj R71. The energy dependent escape of CRs from the Galaxy 
will steepen the spectrum towards an equilibrium particle distribution func- 
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tion NpX dp,injTesc% Qp,inj/DoR-°+ = R75, where s ~ 2.7 is a measured 
quantity (Fig.[2). Secondary CR nuclei (e.g. LiBeB) are produced by spallation 
of heavier primaries by interstellar matter, and their spectrum is steep as it 
mimics the equilibrium spectrum of primaries: gs inj%np. A further spectral 
steepening is induced by CR escape, to give an equilibrium spectrum of secon- 
daries NsLGs, inj TescXNp/ DRT), In fact, this latter result applies only at 
very large energies (larger than those shown in Fig. [2]and |B). The steepness of 
Ns is less pronounced, but still clearly visible at lower energies. These simple 
considerations explain the difference in the spectra of primary and secondary 
CR nuclei shown in Fig. |2| Moreover, knowing that s ~ 2.7 and 6 ~ 0.3...0.6 
one can deduce that at relativistic particle energies the injection spectrum of 
CRs in the ISM must have a slope a ~ 2.1...2.4. This is a remarkable result, 
as it has been obtained without making any assumption on the nature of CR 
sources, except for the fact that they have to be located in the Galactic disk. 

It is now well understood that the confinement of CRs in the Galaxy is due 
to the scattering off the turbulent interstellar magnetic field, which deflects 
the trajectory of CRs and prevents them to escape quickly from the Galaxy 
(e.g. (2013). The scattering is so effective to keep the CR particle 
distribution function very close to isotropy, in agreement with what is observed 
at the Earth’s location. The isotropization of the trajectories of CRs is the main 
obstacle in the search for the sources of such particles, as the arrival direction 
of a CR does not point to the position of its accelerator. As a consequence, 
the production of CRs at many discrete Galactic sources results in an almost 
isotropic flux of particles at the Earth. As we will see in the following (Sec.[4-Ip, 
indirect observations of CRs have to be employed to overcome this problem. 
Instead of observing CRs, one can search for neutral particles, most notably 
photons, which are produced by CR interactions with matter or radiation 
at (or in the vicinity of) the acceleration site. As neutral particles are not 
deflected by magnetic fields their arrival direction pinpoints the position of 
the CR source. Unfortunately, despite many decades of indirect searches, we 
still do not have a conclusive answer to the problem of the origin of CRs. We 
postpone a discussion of this crucial issue to Sec. [4.1.3] 


3 Difficulties in the direct observations of the local interstellar 
spectrum of low energy cosmic rays 


The brief introduction on CRs provided in the previous section was mainly 
focussed on particles of energy = 1 GeV/nucleon. We discuss now the main 
observational difficulties encountered in the study of particles of comparable 
and lower energy, which here are referred to as LECRs. To do so, it is manda- 
tory to distinguish between measurements of the local (near Earth) and remote 
(at some given location in the Galactic disk) intensity of low energy particles. 
As we will see in the following, recent progresses in diverse branches of astro- 
physics and space sciences allowed us to overcome, or at least strongly mitigate 
such difficulties and to obtain an unprecedented view on LECRs. 
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3.1 Solar modulation 


The spectrum of CR nuclei observed at the top of the atmosphere exhibits 
a broad peak at transrelativistic (~ GeV) particle energies, and falls very 
steeply at higher energies (see coloured data points in Fig. (4). For this reason, 
ultra-relativistic particles contribute very little to the CR population both in 
terms of number of particles and of total energy. While the observed intensity 
of ultra-relativistic CRs is remarkably stable in time, that of lower energy par- 
ticles is not. The variation in the LECR intensity exhibits a main periodicity 
of 11 years, indicating a clear link to solar activity. This effect is illustrated 
in Fig. |4| which shows near-Earth spectra of CR protons measured close to 
minima (blue), maxima (red) and intermediate levels (yellow data points) of 
solar activity. The highest (lowest) LECR intensities are measured in corre- 
spondence of solar minima (maxima). This effect is called solar modulation, 
and is interpreted as follows: CRs are originated outside of the solar system 
and, in order to reach the Earth, they have to overcome the solar wind, which 
is an outward flow of magnetised and turbulent plasma emanated from the 
Sun. While high energy CRs can reach virtually undisturbed the Earth, the 
presence of the magnetised solar wind prevents low energy charged particles 
of any species to penetrate freely into the inner heliosphere 
(2013). It can be seen from Fig. [4] that solar modulation is irrelevant 
for ultra-relativistic particles, and becomes important in the trans-relativistic 
and non-relativistic energy domain. 

It follows that the local interstellar spectrum of LECRs cannot be de- 
termined by means of near-Earth observations and this has been, until very 
recently (see Sec. the main obstacle in the study of these particles. CR 
spectra measured at the top of the Earth atmosphere need to be demodu- 
lated to give the local interstellar spectrum. Unfortunately, the demodulation 
of LECR spectra is not at all a straightforward task due to our limited knowl- 


edge of both heliospheric physics (e.g. 1999) and cosmic ray transport 
in turbulent magnetic fields (e.g. 2020). 


3.2 The transport equation 


Technically, a quantitative understanding of solar modulation can be achieved 


by solving the CR transport equation, that was first derived by [Parker] (1965). 
It is a partial differential equation describing how CRs are transported in the 


six-dimensional phase space (x,p). Here, we provide an heuristic derivation 


of such equation following quite closely the approach by (2013). For 
more formal derivations of the transport equation the reader is referred to the 


(1968). 


and to the more recent monographs by 
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Fig. 4 CR proton spectra measured by the Voyagers (white triangles, 


2016| and 2019) and by other instruments at different times. Red, green, and 


blue data points refers to epochs of maximum, intermediate, and minimum solar activity. 
The black line is an analytic representation of the unmodulated CR spectrum, while the 
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Consider a volume V containing N CRs. The rate of variation of N with 
time is given by the continuity equation: 
dN 
—=-0S-d 2 
w=- fs-da+Q (2) 
V 


where the first term on the right is the flux of particles S integrated across a 
close surface containing the volume V, and Q is a source term representing, 
for example, an accelerator operating within the volume V and producing CRs 
at a rate of Q particles per second. The continuity equation can be recast in a 
more familiar form after using the divergence theorem, §S -da = | V - S dV, 
and considering volume densities rather than total number of particles, i.e. 
N = În dV and Q = Sq dV, to obtain: 
on 

or +V-S= (3) 
which can be solved provided S and q are known. Being interested here in the 
study of the heliospheric transport of CRs which are produced outside of the 
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Solar system, we can set q = 0 and impose that at very large distances from 
the Sun n has to be equal to the interstellar density of CRs. 

In the expressions above, n represents the number density of CRs of any 
energy. In fact, what we measure are particle spectra, and so we should 
rather seek for a transport equation for the particle distribution function 
f(p), where p is the momentum of the CR particle. Since the heliospheric 
magnetic field is turbulent, CRs are scattered by magnetic irregularities. The 
scattering is strong enough to keep the particle distribution function very close 
to isotropy, so that we can neglect small deviations from isotropy and write 
n ~ 4r | dp p° f (p) (a discussion on the accuracy of this approximation can be 
found in e.g. [Kirk|/1994). 

Now we have to derive an expression for the differential flux of particles 
S(p) in the solar wind. To a first approximation, the wind can be treated as 
a spherically symmetric outflow of plasma moving at a speed uy. We already 
said that the plasma is magnetised, and that the magnetic field is turbulent. 
The magnetic turbulence can be described as a superposition of magnetohy- 
drodynamic (MHD) waves (for example Alfvén waves). Such waves move with 
respect to the plasma at a velocity v4 « uy, so that they can be considered to 
be virtually at rest in the rest frame of the wind. CRs are coupled to the out- 
ward moving plasma because they scatter off those MHD waves, and this orig- 
inates an outward advective flux of CR particles Sa(p) = f(p)uw. In addition 
to that, particle scattering also results in a spatial random walk of CRs, which 
can be described by a diffusive component in the flux Sa(p) = —D- V f (p), 
where we made use of Fick’s law and introduced the diffusion tensor D. 

However, we notice that while S,(p) represents the advective flux in the 
fixed (solar system) rest frame, Sa(p) describes the diffusive flux in the wind 
rest frame, because it is in that frame that the scattering centres are at rest. 
Thus, it is convenient to introduce a mixed frame, where spatial coordinates 
are specified in the fixed frame, and particle momenta in the wind frame. In 
this mixed frame we can write Stot(p) = Sa(p)+Sa(p) = f(p)uw —D-Vf(p). 

We can now proceed substituting n with f(p) and S with Stot(p) in Eq. 
but this would not result in the correct CR transport equation. The reason is 
that the additional term >> $ (p?<p)f) must be added to the left hand side of 
the transport equation, to account for the fact that particles can gain or lose 
momentum. Here, <p) represents the average rate of momentum gain/loss. 
As pointed out by (2013), the additional term is the divergence of 
the particle flux in momentum space, which is analog to the advective part 
of V-S (with the substitutions uy —> <p), and V > > gp). As noted 
by Parker, in the heliosphere the most important contribution to <p) comes 
from the adiabatic cooling of CRs in the expanding solar wind, at a rate 
(D) = —5V- uw. 

After performing these substitutions and additions, and some simple ma- 
nipulations, the transport equation finally becomes: 


of of 


p 7 
ar Tew VI-V-(D-VF)— 5 (V- uw) 5 = 0 (4) 
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and can be further simplified by setting of = 0, as the typical time scale of 
variation due to solar modulation (a fraction of the solar cycle) is much longer 
than the CR propagation time through the heliosphere (less than a year). 


3.3 Solutions of the transport equation 


The simplest solution of the transport equation can be obtained by neglect- 
ing adiabatic losses (V - uw = 0), by imposing spherical symmetry, and by 
assuming an isotropic diffusion for CRs (the diffusion coefficient in this case 
becomes a scalar function of position and momentum). Under these simplify- 
ing assumptions, and recalling that there are no sources/sinks of particles, the 
transport equation reduces to the diffusion—advection equation: 


Lm 


Uw f — D> 


0 (5) 
where r is the radial coordinate (the distance from the Sun). This simplified 
transport equation describes a situation where an outward advective flux is 
perfectly balanced by an inward diffusive flux. The solution is then: 


f(r,p) = frrs(p) O™ = frrs(p) exp |- i ar | (6) 


where we introduced a non-dimensional modulation parameter M. The local 
interstellar spectrum (LIS) of CRs (f = frrg) is recovered at the border of the 
heliosphere (r = rzrg), while an exponential suppression appears for any r < 
rprs. Under most circumstances, the diffusion coefficient D is a monotonically 
increasing function of particle momentum, and it is then possible to define a 
critical momentum that satisfy M(r,p,) ~ 1. For any fixed value of r, CR 
particles characterised by momenta » p, are unaffected by the presence of 
the solar wind, and can penetrate undisturbed into the inner heliosphere. On 
the contrary, CR particles of low momenta cannot, and for this reason the 
spectrum of CRs observed at the Earth is strongly suppressed with respect to 
the local interstellar one for p « px. 

Let us now proceed a step forward and discuss the most widely used ap- 
proximate analytic solution of the transport equation in the heliosphere. To 
do so, we should remind that f(r,p) is the isotropic part of the particle dis- 
tribution function computed at a position r (in the fixed frame) and at a 
particle momentum p (in the wind frame). When seen from the fixed frame, 
the particle distribution function is no longer isotropic, because both the en- 
ergy and the arrival direction of particles are different in that rest frame. The 
anisotropy adds a contribution to the the advective flux in Eq.|5} which should 
be therefore corrected to give: 


BT 


Cuwf — De 


0 (7) 
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where C = -4S is called Compton-Getting coefficient (derivations of C 


can be found in e.g. |Compton and Getting} |1935} (Gleeson and Axford} |1968 
1970). After substituting C into Eq.|7|we obtain an expression which 


is identical in form to the Liouville equation in a conservative field (e.g. |Fisk 


etal) [1973] (Quenby] 1984): 


of of 
v— +p =0 8 
or Pap (8) 
where v is the particle velocity. The effective “force” p = “8% accounts for (in 


an approximate way) the combined effect of diffusion, advection, and energy 
losses. When written in this form, the transport equation is called force field 


equation (Gleeson and Axford} |1968). 


Another way to interpret the force field equation is to rewrite Parker’s 


transport equation (Eq. |4) in the equivalent form (Moraal} |2013): 


1 
2 pe Ww Vf) =0 (9) 
and notice that the force field equation is obtained after equating to zero 
the term in square brackets. This can be obtained by setting, as done above, 
of = 0, and by further imposing that the last term in the left hand side of the 
equation should be small when compared to both the diffusive and advective 
fluxes. The two conditions to be fulfilled are then: i) uyr/D « 1 which is 
always satisfied provided sufficiently large particle energies are considered (D 
is a growing function of particle energy), and ii) V f/f « C/r, which is satisfied 
more easily in the inner heliosphere, in the vicinity of the Sun, where r is small. 
In other words, the force field approximation applies to particles of high energy, 
that suffer only moderate solar modulation, and therefore are characterised by 
a mild spatial gradient (2004). 
Eq. [8] can be solved by the method of characteristics, that states that the 
particle distribution function is constant along the curve defined by: 


+V-[Cuyf—-D-Vf] + 


dp _ UwP 

dr 3D 
At this point, we should recall that the diffusion coefficient is the product 
between the particle velocity and the mean free path for spatial scattering. 
A convenient way to write it is: D = 35(p)X(p, r), where 8 = v/c is the 
particle speed in units of the speed of light, and À the mean free path multi- 
plied by c (so that it has the same units as D). If the latter is separable, i.e. 
A(p,r) = àp(p)à (r) then we can integrate along the characteristic (Eq. 
and introduce a modulation parameter: 


TLIS PLIS 
dsf arf ape (11) 
r r p 


(10) 


where prrs is the momentum a CR particle had when it entered the helio- 
sphere. The solution of the force field equation can then be written as: 


F(p, r) = frrs(pris(p,r)) - (12) 
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A simple analytic expression for such solution can be found if the diffusion 
coefficient is proportional to particle rigidity A, ~ R = 4, where Z is the 
charge of the particle in units of the (positive) elementary charge e. With this 
assumption, the parameter ¢ has the dimensions of an electric potential (V) 
and is therefore called modulation potential. Moreover, the second integral in 
Eq. [L1]can be performed analytically and gives ¢ = fai -E, where Ezrs and 
E are the kinetic energies of the CR particle when it entered the heliosphere 


and reached the position r, respectively. 

We notice now that it is more convenient to express Eq. in terms of 
the CR intensity, which is the physical quantity which is actually measured. It 
represents the number of particles of a given kinetic energy E passing across a 
unit surface per unit time and unit solid angle, and is connected to the particle 
distribution function through j(F) = p? f (p). It is a simple exercise to obtain 


the very well known result (Gleeson and Axford} |1968): 


E(E + 2Am,) 
(E + 0)(E+6+2Am,) 


H{E,r) = jiis(E + 8) (13) 


where A is the atomic mass number of the CR particle, mp is the proton 
mass, and (r) = Zed(r). The expression above has been widely used in 
the literature, mainly due to its simplicity. If we set r equal to the distance 
between the Sun and the Earth, then only one parameter, ¢, is needed in order 
to connect the local interstellar spectrum of CRs, jzr1s(£), to the one observed 
at the top of the Earth atmosphere, je(F). 


Caballero-Lopez and Moraall (2004) investigated, among others, the limi- 


tations of the force field approximation. Interestingly, they found that while 
this approach provides a good description of the effect of solar modulation 
suffered by LECRs in the inner heliosphere (see also e.g. [Ghelfi et al./[2016), 
the diffusion—advection equation is a more appropriate description at large 
distances from the Sun, where energy losses become less important. In re- 
cent times, continuous observations of CR intensities of unprecedented pre- 
cision became available, thanks especially to space borne detectors (for re- 
views see and [Boezio et al.|[2020). The interpretation of 
such data calls for more sophisticated approaches, involving either modifi- 
cations/generalisations of the modulation potential 
et al.| Gieseler et al./2017), or extensions of the force-field model ( 
and Mertsch| |2019), or fully numerical approaches (Vos and Potgieter} 
Boschini et al 2019). 


Despite the continuous progresses made in the study of particle transport in 
the heliosphere, the local interstellar spectrum of LECRs, as derived through 
a demodulation of the spectra observed at the top of the atmosphere, is in- 
evitably affected by uncertainties. As we will discuss in the following Section, 
direct measurements of LECR spectra beyond the heliopause became recently 
possible, providing us with an unprecedented view on such particles. 
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Fig. 5 Left: integral count rate of CR nuclei and electrons of energy = 70 MeV /nucleon and 
15 MeV, respectively, as measured by the CRS onboard Voyager 1. Figure from [Cummings] 
[et_al.] (2016). Right: integral count rate of CR protons of energy above 211 MeV (A) and 
lower energy CR nuclei, protons, and electrons (B) measured by the LECP instrument 


onboard of Voyager 1. Figure from (2013). Time is on the x-axis. 


3.4 A breakthrough in the measurement of the local interstellar spectrum of 
low energy cosmic rays: the Voyager mission 


On September 5th 1977, the Voyager 1 probe was launched from the Cape 
Canaveral Station. It was preceded by the twin probe Voyager 2, launched on 
August 20th. The goal of the Voyagers was to explore the solar system and, 
following two distinct paths, reach the outer heliosphere and possibly beyond, 
before the end of operations. The study of Galactic CRs and of energetic par- 
ticles in the heliosphere was one of the scientific objectives of the programme 
977). 

More than 40 years after launch, the Cosmic Ray Systems (CRS) onboard 
of the Voyagers are still collecting and sending us data. Each system is com- 
posed of three particle detectors: the High-Energy Telescope System, the Low- 
Energy Telescope System, and the Electron Telescope. These detectors are 
sensitive to nuclei of charge Z = 1 to 30 in the energy range 1-500 MeV 
for protons, 2.5-500 MeV/nucleon for iron, and to electrons of energy 3-110 
MeV (1977). These detectors are complemented by the Low En- 
ergy Charged Particle instrument (LECP), sensitive to particle energies in the 
range 15 keV < E < 40 MeV/nucl, and designed to study energetic particles 
in planetary magnetospheres and interplanetary space [1977). 

The left panel of Fig. |5| shows the integral (nuclei plus electrons) count 
rate as measured over decades by the CRS onboard Voyager 1. Time is on 
the x-axis, and spans from few days after the launch to 2015. Later times 
correspond to larger distances between Voyager 1 and the Earth. The 11 year 
periodicity induced by solar modulation is clearly visible in earlier data and, 
as expected, it becomes less evident at later times, when the probe reaches 
larger distances from the Sun, where the impact of Solar modulation is weaker. 
Three milestones can be identified in this epic journey. On February 17th 1998, 
Voyager 1 reached a distance from the Sun of 69.4 AU and, overtaking Pioneer 
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10, became the farthest human-made object in space. On December 16th 2004 
the probe crossed the solar wind termination shock, at a distance of ~ 94 AU 
from the Sun (label TSX in Fig. B), and entered the heliosheath, which is the 
outermost layer of the heliosphere (Stone et al.|{2005). Finally, on August 25th 
2012 Voyager 1, at a distance of 121.6 AU from the Sun, passed beyond the 
heliopause ( (abel HPX in Fig. |5) and entered interstellar space 
[2013}/Krimigis et al.} [2013). After ae the CR count rate became manly 
stable in time (see inset in Fig. B), as one would expect for CRs unaffected 
by solar modulation. Few years later, on November 5th 2018, the twin probe 


Voyager 2 also left the heliosphere (Stone et al. |2019 2019). 


The scenario just described is very widely, though not universally accepted. 
claimed that the Voyagers might still be inside the 
heliosphere. Moreover, some residual modulation could possibly affect mea- 
surements even outside of the heliosphere, due to the presence of a bow shock 
or of an adiabatic flow ahead of the heliopause (2011). Finally, 
the data collected by the Voyagers revealed that the structure of the helio- 
spheric boundary is more complex than previously thought 2015). All 
these caveats should be kept in mind while reading the rest of this review, 
where we will assume, as a working hypothesis, that the Voyagers are indeed 
located beyond the heliopause and are currently measuring the pristine local 
interstellar spectrum of CRs. 

The working hypothesis adopted here was formulated following the unex- 
pected measurements reported by instruments onboard Voyager 1 in the sum- 
mer of 2012 (Stone et al.| (2013). In particular, the right 
panel of Fig.|5}shows the surprising behaviour of particle count rates measured 
by the LECP instrument for epochs around the presumed heliopause crossing. 
The intensity of low energy ions and electrons of solar origin suddenly dropped 
by more than a factor of 10° on August 25th 2012 and eventually disappeared 
(bottom panel), while the count rate of higher energy CR protons of Galac- 
tic origin simultaneously increased by about 10% (top panel). A very similar 
behaviour was observed by Voyager 2 on November 5th 2018, when the probe 
was at a distance of 119 AU from the Sun. Data may be interpreted by saying 
that, for an observer that crosses the heliopause, the measured intensity of 
energetic particles of solar origin drops as such particles stream away in in- 
terstellar space, while the intensity of Galactic CRs increases as they become 
virtually unaffected by solar modulation. 

Several analytic fits or representations of the local interstellar spectrum of 


CR protons can be Suli in the literature ( Schlickeiser et al. et al. Schlickeiser et al.|[2014}|Ivlev et al. Ivlev et al./ et al. 
[2015} [Phan et al.} [2018] . We adopt here a doubly broken power law functior°| 


12.5 BE 0.35 
jiis(E) = (sev) k m`ĉ2s —1 sr~tMeV~! (14) 


[1+ i) | [1+ aa) ] 


6 Which provides a better representation of data with respect to the singly broken power 
law spectra often adopted in the literature. However, this is just a convenient descriptive 
expression, being neither a physically motivated choice, nor a formal fit to data. 
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where the parameter s can be tuned to vary the shape of the spectrum around 
the spectral break at 2.2 GeV/nucleon. Smaller (larger) values of s make the 
break sharper (shallower). The parameter has been introduced because the lo- 
cal interstellar spectrum of CRs is not constrained in the GeV energy domain, 
as direct near Earth observations are heavily affected by solar modulation. 
Eq. [14]is plotted as a solid black line in Fig.|4} and provides a convenient de- 
scription of both interstellar measurements of CR protons at low particle en- 
ergies (Voyager data, white triangles in Fig. 5), and near-Earth ones (coloured 
data points in Fig.|/4) at high enough particle energies, where solar modulation 
has no effect. To describe the intermediate energy regime, we set s = 2.1, a 
choice that will be motivated in Sec. 


On the contrary, near-Earth proton spectra (coloured data points in Fig. A) 
are shaped by solar modulation for particle energies below few tens of GeV. As 
seen in Section [3.1] the effect of solar wind on LECRs can be accounted for in 
an approximate way by means of Eq. In that approach, the modulation po- 
tential (¢) alone suffices to account for all of the relevant heliospheric physics. 
Before the Voyagers entered interstellar space, the only measured quantity in 
Eq. [13] was the near-Earth spectrum of CRs, jẹ(E). Therefore, a value for ¢ 
had to be chosen/estimated in order to derive the local interstellar spectrum 
of LECRs. Now, after the two Voyager probes have crossed the heliopause, 
both the near-Earth and local interstellar spectra of CRs are measured, and 
Eq. can be used to infer ¢. The results of such a procedure are shown in 
a the dotted blue, green, and red curves represent expectations of 
the near-Earth CR spectra characterised by values of the modulation potential 
equal to ¢ = 0.6, 0.8, and 1.4 GV, respectively. 


What said above illustrates one of the many reasons why the LECR data 
collected by the Voyagers are so important: they provide us, for the very 
first time, with precious and direct constrains on the physical processes that 
regulate CR modulation. More in general, data collected by the Voyagers had 
a tremendous impact on the study of the heliosphere and of its boundary. 
As this aspect goes beyond the scope of this review, the interested reader is 


referred to (2015). 


The Voyagers measured also the local interstellar spectra of LECR nuclei 
heavier than hydrogen. The most abundant amongst them is helium, whose 
spectrum is shown in the top-left panel of Fig. [6] In the Figure, Voyager 
and AMS-02 data are shown as red triangles and blue circles, respectively. 
noticed that the spectra measured by Voyager for 
hydrogen and helium overlap almost perfectly when plotted as a function of 
the particle energy per nucleon, and when the former is divided by a factor 
of 12.2. On the other hand, the spectrum of helium measured by AMS-02 
at high energies is appreciably harder than that of hydrogen (both can be 
described by power laws in energy, the difference between the spectral slopes 
being ~ 0.07). Following these indications, the black line in the top-left panel 
of Fig. [6] has been obtained by dividing Eq. [14] by 12.2 and substituting [1 + 
(E/2.2 GeV)@-9/2-D]2-1 with [1 + (£/1.6 GeV) (-83/1:8)]1.8. The reason why 
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Fig. 6 Top left: CR helium spectrum as measured by Voyager (Cummings et al. 2016 
red triangles) and AMS-02 (Aguilar et al.| [2017| blue circles). Bottom left: CR electron 


spectrum, as measured by Voyager (Cummings et al.| [2016| red triangles) and AMS-02 


(Aguilar et al.| {2019| blue circles). Top right: CR intensity for protons (blue), helium 
g ( 


reen), and electrons (red curve) in the local ISM. Bottom right: CR energy densities in 
the local ISM. Same color code as in top right panel. 


the spectra of protons and helium nuclei are different in the multi-GeV energy 
domain is not understood (see discussion Secret 

Finally, the bottom-left panel of Fig. |6| shows the local interstellar CR 
electron spectrum. Also in this case, red triangles and blue circles show Voyager 
and AMS-02 data, respectively. As for protons, a doubly broken power law 
function 


5 x 10 (2) " 


0.6/6; | 1.38/59 102 
[i+ (a) | + (are | 
(15) 


with 6; = 0.2 and 62 = 0.5 provides a good description of data and is plotted as 
a black solid line in the Figure. We should stress here that Voyager and AMS- 
02 data alone would not require the presence of two breaks: a singly broken 
power law would provide an equally good fit to data. The choice made in Eq, [15] 
will become clear in Sec. when we will discuss indirect measurements of 
the local CR electron spectrum. 

In order to evaluate the impact that LECRs exert on the ISM, it is in- 
structive to compute and plot the quantities R;(E) = 47 Ej;(E) and w;(£) = 
E?nj(E) = 4rE*j;(E)/v(E), where n;(E) is the particle number density, v(£) 


m’s~!sr7!Mev7! 


jiis(£) = 
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the velocity of a particle characterised by an energy per nucleon E, and the 
suffix į may refer to protons, helium nuclei or electrons. These two quantities 
are shown in the top-right and bottom-right panel of Fig. a respectively. 

LECRs interact with the interstellar gas in many ways. If o;(E)xE~* is 
the cross section of a given physical process j (for example, the ionisation cross 
section), then the quantity R;(£)o;(£) represents the rate of interactions per 
interstellar atom suffered by LECR particles of energy per nucleon within 
a given logarithmic interval around E. It can be seen from Fig. [6] that the 
curves for R(E) relative to protons and helium are very peaked around a 
particle energy equal to Epeak < 1 GeV/nucleon. For particle energies below 
the peak R;(E)ocE'?°, which implies that for “well behaved” cross sections 
characterised by s < 1.35, the total rate of interactions (integrated on all 
energies) due to the process j is dominated by particles of energy around 
~ Epeak- As the condition on s is satisfied by ionisation cross sections (Cion © 
1/E, see Sec.|5.3), the CR ionisation rate in the local ISM is well constrained 
by Voyager and AMS-02 data for CR nuclei (which are largely dominated 
by protons and helium). Unfortunately, this is not the case for electrons, for 
which R;(£) decreases monotonically in the energy domain probed by available 
observations. The CR ionisation rate is therefore not well constrained for CR 
electrons, as it depends on the behaviour of their spectrum for energies smaller 
than few MeV, where no observations have ever been performed. In this case, 
only a lower limit for the ionisation rate in the local ISM can be obtained. 

On the contrary, as shown in the bottom-right panel of Fig. [6] the total 
energy density of CRs in the local ISM is well constrained for both nuclei and 
electrons. The largest contribution to the local energy density is provided by 
particles of energy = 1 GeV/nucleon, and integrating over all energies one 
gets energy densities of ~ 0.55, 0.16, and 0.03 eV/cm® for CR protons, helium 
nuclei, and electrons, respectively. The energy density of nuclei heavier than 
helium is roughly twice that of electrons (2016), giving a 
total CR energy density in the local ISM of the order of ~ 0.8 eV/cm?. 

The statements just made are valid under the (very plausible, see Sec. 
hypothesis that there are no additional components in the CR spectrum that 
emerge at very low energies and dominate the total CR energy budget in 
the local ISM. The presence of such hidden components in the CR spectrum 
has been sometimes invoked in the past in order to explain some presumed 


anomaly in the cosmic abundance of light nuclei (Meneguzzi et al. |1971) 


or, more recently, as an additional source of ionisation of the interstellar gas 
(2016). Before the Voyagers entered the local ISM, such 
ad hoc low energy component was assumed to be present in the MeV energy 
domain, and to be hidden by solar modulation, while now that measurements 
of the pristine interstellar spectrum of CRs extend down to few MeV, an 
hidden component can only show up in the sub-MeV domain. 

Before concluding this Section, few words on heavier CR nuclei are in order. 
Here, we will not discuss in great detail the spectra of primary LECRs heavier 
than helium, but we will just highlight the two most important observational 
findings. First of all, the striking similarity of the spectral shape of protons 


Low energy cosmic rays 25 


and helium CRs observed by the Voyagers does not extend to heavier nuclei, 
which exhibit much harder spectra in the MeV domain 
(2016). Second, the spectra of CR carbon and oxygen (the two most abundant 
nuclei heavier than helium) have been measured by the AMS-02 and found 
to be very similar to that of helium [2017). Therefore, in the 
AMS-02 energy domain, the most abundant primary CR nuclei seem to have 
very similar spectra, with the exception of protons, whose spectrum is slightly 
softer. Conversely, in the lower energy domain probed by the Voyagers, proton 
and helium have very similar spectra, which are in turn different by those of 
heavier nuclei. Both these differences are not understood, and the interested 
reader is referred to[Gabici et_al.J2019]and[Tatischeff et_al.J2021]and references 
therein for a discussion of this very puzzling issue. As particle acceleration 
mechanisms are not expected to distinguish between species when shaping 
their energy spectra, it is tempting to interpret such differences as the result of 
a different origin for different species (e.g. acceleration taking place in different 
phases of the ISM, or in different astrophysical objects, or at different evolutive 
stages of the same objects, etc.). 

Finally, the measurements of the local interstellar spectrum of LECRs re- 
viewed in this Section are also very important in connection with the problem 
of the origin of CRs. With this respect, the main goal is to understand how 
the features observed in the local interstellar spectrum of LECRs could be 
explained as the result of the injection of CRs in the ISM, followed by their 
transport in the turbulent interstellar magnetic field (Strong et al.}[2007). This 
is a very complex problem, whose solution requires a knowledge of the nature 
and spatial distribution of CR sources in the Galaxy, of the topology of the 
interstellar magnetic field over a huge range of spatial scales, of the details of 
the transport of CR particles in turbulent fields, and on the spatial distribu- 
tion of interstellar matter. In the simple model of CR transport in the Galaxy 
described in Sec. the problem was brutally simplified by adopting mean 
particle densities or residence times in the Galaxy. Spatial variations in the 
CR intensities cannot be neglected in more realistic models, and is therefore 
important to review the current status of measurements, inevitably indirect, 
of the intensity of CRs in regions far from the Solar system. This will be done 
in the next Section. 


4 Indirect measurements of the remote interstellar spectrum of low 
energy cosmic rays: gamma-ray and radio observations 


4.1 Hadronic gamma rays 


In 1952 Hayakawa first proposed that, if CRs fill the entire Galaxy, the galactic 
disk should shine in gamma rays. Such gamma rays are generated by the 
decay of light mesons, which are produced in the inelastic collisions between 
CR nuclei (mainly protons) and nuclei in the ISM (mainly hydrogen). In this 
context, the dominant channel for gamma-ray production consists of two steps. 
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First, an energetic CR proton interacts with a proton at rest in the ISM, to 
generate a neutral pion, which immediately decays into two gamma rays: 


ptp—X4+n° , P —y++y. (16) 


In the expression above, X represents all the other products of the interac- 
tion, possibly including more pions. The process is characterised by an energy 
threshold Etn, which can be easily estimated by assuming that the final prod- 
ucts of Eq. are two protons and a neutral pion at rest, i.e. the minimum 
energy configuration. It is straightforward to show that only CR protons of 
kinetic energy larger than Eth = 2m,0 +m? /2Mp © 280 MeV are able to gen- 
erate neutral pions (here and in the following of the Section, c = 1). Neutral 
pions are also produced in collisions involving nuclei heavier than hydrogen, 
which are present in both CRs and in the ISM. Taking this contribution into ac- 
count enhances the predictions of the gamma-ray signal by a factor of ey < 2, 
which depends weakly on the spectrum of the incident CR nuclei 
Kafexhiu et al] [2014). 

As gamma rays can pass through large column densities of matter with- 
out being absorbed, the implications of Haykawa’s proposal are far reaching: 
gamma-ray observations allow us to observe, indirectly, CR protons and nu- 
clei of energy larger than few hundreds MeV located in remote regions of 
the Galaxy. Hayakawa’s conjecture was confirmed in the late sixties, when 
gamma rays of energy exceeding 100 MeV were observed from a band in the 
sky spatially coincident with the Galactic disk (1968). The diffuse 
emission from the Galactic disk is, in fact, the most prominent feature in the 
gamma-ray sky in the GeV domain, and has been measured over the years 
with ever-improving accuracy, most recently by the Fermi/LAT space tele- 
scope (2012). The spectral shape of the emerging spectrum 
of gamma rays from inelastic proton proton interactions can be computed from 
pion production models and/or fits to accelerator data (see e.g. 
and the comprehensive bibliography in{Kafexhiu et al.[2014), and a number of 
convenient analytic parameterisations of gamma-ray production spectra can 
be found in the literature (e.g. 
‘Koldobskiy et al.||2021). Remarkably, the main features of such 


gamma-ray spectra can be inferred from the following two simple arguments. 

First of all, in the rest frame of the neutral pion, the two gamma rays 
produced in the decay will have the same energy EX =m,0/2 ~ 70 MeV and, 
to conserve momentum, they will move along opposite directions. After Lorentz 
transforming to the lab frame, the photon energies are boosted to E+ = (EX an 
Bp, cos 0*), where y and 6 are the Lorentz factor and velocity of the pion, and 
p>, and §* the momentum of the photons and the decay angle in the pion rest 
frame. In the absence of any preferred direction, the distribution of the decay 
angles in the pion frame has to be isotropic. This implies that, for an ensemble 
of pions of a given energy, the emerging spectrum of gamma rays in the lab 
frame is flat, i.e. dn,/dE,«adn,/dcos ¥* = 0, in the energy interval defined by 
cos 0* = +1. If £7” and ae are the maximum and minimum energies of the 
photon distribution, corresponding to cos?* = +1 and -1, respectively, then 
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Fig. 7 Hard-X and gamma-ray emission from the inner Galaxy (|I| < 30° and |b| < 10°). 
Curves show the expected contribution from sources (lower dot-dashed magenta line), dif- 
fuse emission (blue line), and total (upper dot-dashed magenta line). The diffuse emission 
is the sum of the contributions from 7° decay (red), inverse Compton scattering (green), 
Bremsstrahlung (cyan), isotropic extragalactic (black). Data are from Fermi (black), INTE- 


GRAL (magenta and blue), and COMPTEL (green). Figure from|de Angelis et al.| (2018). 


their geometric mean is always equal to 4/ Emer Emin = m,0/2, regardless 


of the pion energy. It follows that, for pions characterised by an arbitrary 
distributions of energies, the gamma ray spectrum will be the superposition of 
many flat spectra which are all symmetric, when plotted versus log(E,), with 
respect to the logarithm of the geometric mean log(m,o/2). This symmetry is 
a distinctive feature of gamma-ray spectra from pion decay, and the peak that 
invariably appears in the spectrum at E = m,o/2 ~ 70 MeV is called pion 
bump. 

The second important aspect of proton-proton interactions is that, even 
though a large number of neutral pions can be produced in each collision, for 
large enough proton energies (> Eth) one leading pion carries away the main 
fraction (of the order of ~ 20%) of the CR proton kinetic energy. Moreover, in 
the same limit, the cross section of inelastic scattering depends quite weakly 
on the energy of the incident particle. Therefore, if the energy spectrum of CR 
protons is a power law Ep ”, then also the pion and the gamma-ray spectra 
will be power laws Ery” with similar slopes sp ~ 8, ~ sy. In fact, a more 
accurate analysis shows that the gamma-ray spectrum does not follow a pure 
power law, but progressively hardens above ~ 100 GeV, where s œ% sp — 0.1 


(Kelner et al.||2006). Finally, well above the threshold (E, > Ef"), the energy 


of the gamma ray is linked to that of the parent CR proton as Ey ~ 0.1 x Ep. 
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Fig. [7|shows the diffuse emission observed from the inner part of the Galac- 
tic disk (data points) over a very broad energy range, spanning from the hard 
X-ray to the high-energy gamma-ray domain. The emission has been extracted 
from a region around the Galactic centre, defined by Galactic coordinates 
|b] < 10° and || < 30°. Data points come from observations performed by 
INTEGRAL and COMPTEL in the hard X-ray/soft gamma-ray domain (ma- 
genta, blue, and green points) and by Fermi (black points) in the high energy 
gamma-ray band 2011). The observed diffuse emission has been fit- 
ted by a multi-component model, that will be described in more detail below. 
For the moment, we only anticipate that in the GeV domain the emission is 
largely dominated by photons from neutral pion decay, represented by the solid 
red line: a striking validation of Hayakawa’s early prediction! In the Figure, 
the spectrum has been multiplied by the photon energy squared, to show the 
spectral energy distribution. In this convenient representation, the spectrum of 
gamma-ray photons from proton-proton interactions is no longer symmetric 
with respect to EY}, and its peak is shifted to £ GeV energies. 

Besides the spectrum, also the morphology of the GeV diffuse emission 
from the disk carries important information, as it is shaped by the spatial 
distribution of both CRs and interstellar matter. As the latter can be derived 
by a number of observations (Ferriére| [2001), the former can be constrained by 
gamma-ray observations. In the following, we will briefly describe how this can 
be done, starting with a review of the available estimates of the local (within 
~ 1 kpc) intensity of CRs, and proceeding then with a discussion of various 
approaches aimed at probing the remote regions of the Galaxy. 


4.1.1 Constraints on the cosmic ray intensity from diffuse gamma rays 


The expected local gamma-ray emissivity due to neutral pion decay can be 
estimated starting from the CR proton spectrum measured in the local ISM 
(Eq. [14p. It represents the number of gamma-ray photons emitted per unit 
energy, time, and per interstellar hydrogen atom. It is defined as 


1971): 


o0 
g15(B,) = an | dEyo(By|Ep)in1s(Ep) (17) 
Ey 


where o(E,|E,) is the differential cross section (parameterised in, e.g. 
describing the probability to have an interaction where a proton 
of energy Ep produces a gamma ray of energy E,. In general, the intensity of 
CR protons at a position r in the Galaxy, j(Ep,r), will differ from the local 
one (r = 0), and the gamma-ray emissivity q,(E,,r) will change accordingly. 
The observed intensity of gamma rays from a given direction in the sky is then 
obtained after integrating along the line of sight: 


(Ey: b) = E | ar (Eypr)nasa(r) (18) 


where l and b are the Galactic longitude and latitude, respectively, and nzsm(r) 
is the hydrogen (atomic plus molecular) density at a distance r along the line 
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of sight. The spatial distribution of interstellar matter nzsm(r) can be de- 
termined, with some uncertainties, from astronomical observations 
(2005). Then, if the intensity of the diffuse gamma-ray emission 
jy(Ey, 1, b) is measured, constraints on the spatial distribution of CRs through- 
out the Galaxy can be obtained. 

In order to check the prediction given in Eq. it is convenient to observe 
the diffuse gamma-ray emission at large Galactic latitudes. This is because, due 
to the small thickness h of the gaseous Galactic disk (few hundreds parsecs), 
the gamma-ray emission from such latitudes is dominated by the contribution 
of CR interactions taking place in the solar neighbourhood, i.e. within a dis- 
tance ~ h/(2sinb). Under these circumstances, Eq. can be simplified by 
making the substitution q,(E,,r) —> ce °(E,), which is acceptable provided 
that there are no large spatial gradients in the very local distribution of CRs. 
This gives: 


5 € E€ 
JEt b) x Takis (E,) | ar mrsm) = Zaks (E Nu) (9) 


where Ny is the hydrogen column density. This approximate expression shows 
how, from the measurement of the diffuse gamma-ray emission and of the gas 
column density, it is possible to obtain an observational estimate of ae (Ey) 
that can be in turn compared with its predicted value (Eq. [17). This can be 
done to a good accuracy because, as we said above, the neutral pion decay 
contribution to the diffuse emission from the Galactic disk exceeds that from 
other radiation mechanisms for photon energies exceeding ~ 100 MeV (Fig.[7p. 

analysed Fermi data to extract the diffuse gamma-ray 
emission at large Galactic latitudes in the range 10° < |b| < 70°. Most of 
this emission is generated by proton-proton interactions taking place within a 
distance of roughly a kiloparsecs from the Sun. The contribution to the diffuse 
gamma-ray emissivity coming from molecular, atomic, and ionised hydrogen 
can be separated as the gas column densities of these three components of the 
ISM are known with fair accuracy. Once this is done, the derived interstellar 
gamma-ray emissivity for atomic hydrogen is taken as the most reliable, as the 
spatial distribution of such component of the ISM, traced by the 21 cm line, is 
known with high precision. An estimate of the local spectrum of CRs can then 
be obtained from such emissivity. Besides the complications in the determina- 
tion of the gas column densities, such method also suffers from uncertainties 
in both the knowledge of hadronic cross sections, and in the determination 
of the contribution from leptonic processes (see Sec. to the gamma-ray 
emission. While different values for the gas column densities of the phases of 
the ISM would shift the overall normalisation of the resulting CR spectrum 
in an almost energy independent way, the latter two effects have a greater 
impact at low energies (see the recent review by for a 
more extended discussion of these issues). 

The gamma-ray-based approach just described is important as it probes 
the local spectrum of CRs in the energy region where direct observations are 
not available. For CR protons, this spans from the highest energy data point 
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Fig. 8 Left: Direct measurements of CR protons (white triangles: Voyager 1 and 2; green 
circles: AMS-02) are compared with indirect constraints from gamma-ray observations of 
the local (~ 1 kpc) ISM (the red solid line, red dotted line, and yellow shaded region are the 
best fit, the best fit divided by 1.35, and the uncertainty, respectively, as in [Strong]2015). 
Right: Direct measurements of CR electrons (red triangles: Voyager 1; blue circles: AMS- 
02) are compared with constraints from radio and gamma-ray observations of the local (~ 
1 kpc) ISM and from local CR observations (red dotted line from 2018). In both 
panels, the black line represents the parameterisation adopted in this review. 


obtained by the Voyagers, at about 300 MeV, to few tens of GeV, i.e. the 
minimum energy for which CR measurements performed inside the heliosphere 
are unaffected by Solar modulation. Remarkably, a comparison between direct 
and indirect measurements of the CR intensity can tell us whether a gradient 
in the spatial distribution of CRs exists or if CRs are uniformly distributed 
in a kpc neighbourhood of the Solar system. A strong spatial gradient would 
indicate the presence of one or more local sources of CRs, and would put into 
question the scenario for CR transport and escape from the Galaxy developed 
in Sec. |2} which is based on the assumption that the local spectrum of CRs 
can be taken as a proxy for the typical one in the entire disk. 


Building on the early work by {Casandjian| (2015), improved 


the analysis technique and adopted more updated cross sections, and {Orlando| 
revised the estimate of the leptonic contribution to the gamma-ray 
emission. The results from these studies indicate that at high enough particle 
energies, where the effects of Solar modulation are negligible and hadronic 
models are most accurate, the spectrum derived from gamma-ray observations 
exceeds by ~ 30 — 40% direct measurements. This is shown in the left panel of 
Fig.[8} where direct observations of CRs are compared with indirect ones. In the 
plot, the local interstellar spectrum of CRs measured by the Voyagers (white 
triangles) is shown together with the near Earth observations by AMS-02 
(green circles). The solid red line represents the best fit CR spectrum derived 
by the local (< 1 kpc) gamma-ray emissivity, with its uncertainty shown as a 
yellow shaded region 2015). If the best fit model is shifted downward 
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Fig. 9 Top: Normalisation of the CR proton spectrum at 30 GeV as a function of Galac- 
tocentric distance as derived by gamma-ray observations. Red, green, and blue data points 


are from 2016 2021) and 2020| respectively. The 


y. 
horizontal dotted line indicates the value measured locally by AMS-02 (Aguilar et al.|/2015). 
Figure adapted from (2021). Bottom: CR ionisation rate versus Galactocentric 


distance. Data are from (2016) and (2019) (green point), |Neufeld 
and Wolfire| (2017) (blue points). The yellow shaded region shows the expected value for the 


local interstellar spectrum of CRs, after the effects of the penetration into a diffuse cloud 
are taken into account (Phan et al.|!2018). The red bar shows the position of the Sun. The 


dashed black line shows the profile predicted by|Wolfire et al.| (2003) renormalised to fit the 


blue data points. 


by a factor of 1.35 (red dotted curve), it smoothly connects the Voyager data 
points to the ones from AMS-02. The black solid line in the Figure represents 
the expression given in Eq. [14] for s = 2.1. For such a choice of the parameter s, 
the expression constitutes a good description not only of direct measurements 
at low (Voyagers) and high (AMS-02) energies, but also reproduces the correct 
shape of the spectrum (though with a slightly different normalisation) in the 
GeV domain, as derived from gamma-ray observation of the local ISM. 


At low Galactic latitudes, the interpretation of the observations of the 
GeV diffuse gamma-ray emission is less straightforward, as in this case the 
approximation made to derive Eq. [19] is no longer valid, and the full expres- 
sion given in Eq. [18] must be adopted. Nevertheless, both the diffuse emission 
jy(E,,1,b) and the space-dependent interstellar gas density nzsm(r) can be 
determined from observations, turning Eq. [T8]into an integral equation for the 
CR intensity j,(E,,r). The inversion of the integral equation is customarily 
performed dividing the Galactic disk in concentric rings around the Galactic 
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centre, and imposing that the intensity of CR nuclei is spatially uniform within 
each ring. The normalisation and spectral shape of the CR intensity are left 
as adjustable parameters and may differ in each ring. Their best fit values 
are obtained by comparing the expected and the observed maps of the diffuse 
gamma-ray emission by means of a likelihood analysis, whose details can be 
found io, e, (2016) [Yang ct al] (2016) and [Pothast et al] (2073) 
The results obtained by (2016) for the density of CR protons as 
a function of the galactocentric radius R are shown as red data points in the 
top panel of Fig. p The values of the CR density obtained in this way for 
R g 15 kpc differ from each other by less than a factor of ~ 2.5, with the sole 
exception of the data point centred at R = 3 kpc, which stands a factor of 
~ 3 above neighbouring points. It has been suggested that this excess might 
be related to an enhanced CR production in the molecular ring, a region posi- 
tioned roughly half-way between the Galactic centre and the Sun, where both 
the distribution of interstellar molecular hydrogen and the star formation rate 
per unit surface of the disk peak (2016). Note also that the CR 
density inferred within the ring that contains the Solar system (that centred 
at R = 9 kpc) exceeds by a factor of ~ 1.6 the value of the local CR den- 
sity as measured by AMS-02, thereby confirming the earlier claim by |Strong 
(2015) that we discussed above. It seems, then, that large scale variations in 
the density of CRs are quite mild in the Galactic disk. 


4.1.2 Gamma rays from MCs 


The procedure of inversion of Eq. just described can only provide the in- 
tensity of CRs averaged over quite large volumes of the disk, for example a 
ring centred at a Galactocentric radius R and of width AR of the order of a 
kiloparsecs or more (see the red error bars in the top panel of Fig. [9). Inter- 
estingly, gamma-ray observations of MCs can be used to probe the density of 
CRs within much smaller volumes, as these objects’ typical radii range from 
parsecs to tens of parsecs (Heyer and Dame} {2015). The other distinctive char- 
acteristic of MCs is their large mass, which provides abundant target for CR 
proton-proton interactions. For this reason, MCs were proposed as potential 
discrete sources of gamma rays in a seminal paper by [Black and Fazio| (1973). 
Of particular interest are giant MCs, i.e. those with masses Mea exceeding 
~ 104Mo. Their masses and radii correlate in such a way that such objects 
are characterised by a typical column density Ny = 107? cm~? 
[Ostriker} [2007). This is of the same order of the total average gas column den- 
sity measured along lines of sight in the direction of the inner Galaxy, which 
explain why the gamma-ray emission from giant MCs is expected to be seen 
as an excess above the diffuse interstellar emission, provided that the intensity 
of CRs is similar in the diffuse ISM and inside clouds. 

The gamma-ray luminosity from a MC of mass Mea can be computed by 
integrating Eq. [18] over the volume of the cloud, and its flux then reads: 


cl 
q (E~) Maea 
FOE.) = si 2 
y (E) = EN rym, (B ~ 
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where da is the distance to the cloud, and yz ~ 1.4 accounts for the presence of 
helium in the cloud. In deriving the expression above, we implicitly assumed 
that the spatial distribution of CRs does not vary significantly over the scale of 
the cloud, that is to say, the intensity of CRs is not attenuated due to ionisation 
energy losses as they go through the cloud’s large column density (note that 
the entire Sec. [6.3] will be devoted to the issue of CR penetration into MCs). 
We also assumed that the hadronic emission dominates the gamma-ray flux, a 
condition easily satisfied if photon energies larger than few hundreds of MeV 


are considered (see e.g. |Gabici et al.| [2007| |2009). Under these conditions, 


qs! (Ey) is the hadronic gamma-ray emissivity at any location inside the MC. 


Eq. [20]shows that if the mass and distance of a MC are known, the gamma- 
ray emissivity per hydrogen atom qS (Ey) can be directly estimated from its 
gamma-ray flux BOB i The density and energy spectrum of CR protons 
inside the cloud can be in turn derived from the gamma-ray emissivity. This is 
why MCs have been sometimes referred to as cosmic ray barometers: we can 
use them to probe the density and pressure of CRs in remote regions of the 
Galaxy (Issa and Wolfendale] [981] [Nharonian| [1901] [Casanova etal 2010) 
The uncertainties in the measurements of the cloud mass and distance are the 
main limitations of this approach. 


MCs are almost entirely made of molecular hydrogen, H2, which unfor- 
tunately is difficult to detect{"| For this reason, the mass of a cloud is often 
estimated from observations of the second most abundant molecule in space, 
carbon monoxide (CO). The CO emission line corresponding to the rotational 
transition J = 1 — 0 falls in the radio domain (A = 2.6 mm) and is the 
prime tracer of molecular gas (2015). As the transition is 
optically thick, the CO brightness temperature integrated over the line profile, 
Wco, is not proportional to the CO column density, but rather to the cloud 
velocity dispersion. For virialized clouds, the velocity dispersion is related to 
their total mass, which is dominated by Hə (Scoville and Sanders} {1987). The 
proportionality constant between the Hy column density N(H2) and the CO 
integrated line intensity Wco is usually called the Xco factor and has been 
determined empirically, though with a significant average (mediated over large 
scales) uncertainty of ~ 30%, with individual clouds fluctuating around that 
value (2013). The infrared emission from dust grains mixed 
with the gaseous ISM is another tracer of interstellar matter, providing com- 
plementary estimates of the cloud masses. The advantage of this approach 
is that dust emission is optically thin, and therefore directly traces the total 
gas mass of the cloud. However, it also suffers from uncertainties and cloud- 
to-cloud variations (e.g. (2021). Finally, the 
observed frequency of the CO line (or any other line) is doppler shifted due to 
the differential rotation of the Galaxy, a thing that can be used to estimate the 


T The main problem is that even the lowest (rotational) excited energy level is too far 
from the ground state to be significantly populated, given the cold environments where H2 


is generally found (see|Stahler and Palla} [2004| for a detailed discussion). 
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distance of a cloud, provided that it is located within the solar circle 
2o07). 

A large number of studies on local MCs were aimed at determining their 
CR content by combining gamma-ray observations and mass estimates (see e.g. 
for an early attempt, and the list of recent references 
compiled by Tibaldo et al./2021). The gamma-ray emissivities derived in this 
way for clouds located within few hundred parsecs from the Solar system 
show a moderate scattering around that derived by (2015) using 
diffuse gamma rays (see e.g. Fig. 10 in or Fig. 2 in 
fet al.[2021), or around the emissivity computed from AMS-02 data (see blue 
data points around R ~ 8.5 kpc in the upper panel of Fig. [9] or [Aharonian| 
et al./2020). Given the intrinsic uncertainty in the determination of the matter 
distribution in our neighbourhood, it is not possible to claim whether direct 
and indirect measurements of the local interstellar spectrum of CR protons 
agree or not. However, if a difference exists, it must be quite small, as data 
constrain local spatial variations within few tens percent nhac hecho 

The blue and green data points in the upper panel of Fig.|9}show the density 
of CRs as derived from observations of both local and remote MCs 
fet al.| 2020} {Peron et al.| {2021}. Overall, they broadly agree with the results of 
studies of the diffuse emission (red data points in the same Figure), except for 
the region around the molecular ring (green points in the Figure). Values of 
the CR density extracted from MC observations are closer to direct AMS-02 
measurements with respect to those derived from the diffuse gas. In fact, this 
might not be a problem, as MCs probe the intensity of CRs at very specific 
locations, while studies of the diffuse gas provide a volume average quantity 
(for an extended discussion of this issue see {Peron et al.| (2021). 

In any case, independently on all these complications, a claim that can 
be certainly made is that the spatial distribution of CR protons of multi-GeV 
energies varies very little on large spatial scales. More quantitatively, variations 
are constrained to be at most a factor of two from the value measured locally 
by AMS-02, with the only exception of the molecular ring. However, also in 
that case the variation would be of at most a factor of ~ 5. 


4.1.8 The Galactic origin of cosmic rays: the supernova remnant paradigm 


Having determined the spatial distribution of CR nuclei in the ISM, we can 
now estimate their total energy in the Galactic disk. As seen in the previous 
Section, the gradient of the large scale distribution of CRs is quite modest. 
Therefore, an order of magnitude estimate of the CR energy in the disk is 
given by Wor ~ wh x Vaisk © 2 x 10°° erg, where wt’ ~ 1 eV/cm? is the 
energy density of CR nuclei in the local ISM (Sec. (3.4), and Vaisk ~ 400 kpc? 
is the volume of the Galactic disk of radius 15 kpc and thickness 500 pc. 
The largest contribution to the CR total energy comes from nuclei of energy 
= 1 GeV/nucleon (Fig. (6). Before escaping the Galaxy, particles having such 
energy spend Tzsm < 10 Myr in the ISM (Sec. 21). To compensate for the 
escape of CRs and maintain stationarity, sources located in the disk have 
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to inject in the ISM a power Por ~ Wer/trsm S 104! erg/s in form of 
energetic particles. This estimate is admittedly very rough, but provides a 
result consistent with accurate studies (see e.g. (2010). 

Following a seminal idea proposed by (1934), it was 
soon realised that Galactic supernovae could provide the amount of energy 
needed to explain CRs (1950). The rate vgn of supernova explosion 
in the Galaxy is equal to about three per century, and each explosion releases 
Esn ~ 10°! erg in form of kinetic energy of the stellar ejecta. Therefore, 
the total power injected by supernovae in the Galaxy is Psy = vgsnEgn ~% 
10* erg. It follows that about < 10% of the energy from supernovae has to 
be converted into energetic particles in order to explain CRs. 

With the advent of radio observations it was realised that the remnants 
of supernova explosions, and not the explosions themselves, were likely to be 
better candidate as sources of CRs. Supernova remnants are the expanding 
shocks that form in the ISM as a result of the stellar explosion, and the syn- 
chrotron radio emission detected from these objects implies the presence of 
magnetic fields and relativistic electrons there. Finally, in the seventies, the 
theory of diffusive shock acceleration was developed, independently, by four 


different teams 
[1978} [Bell] [1978), showing that strong shocks are indeed expected to accelerate 
particles with an universal spectral shape: a power law in energy of slope E~?. 

As the slope predicted from diffusive shock acceleration theory is quite close 
(though not equal) to what is needed to explain CR observations (slopes in 
the range 2.1 ... 2.4 are compatible with observations, see Sec. a scenario 
where CRs are accelerated at supernova remnant shocks became so popu- 
lar that is nowadays considered a paradigm (see and references 
therein). However, two things should be kept in mind: first of all, the supernova 
remnant paradigm for the origin of CRs remains, in fact, an hypothesis 
and, second, the energy-budget argument provided above holds 
for particles of energy around ~ 1 GeV, which are those who dominate the 
energetic. It follows that the origin of LECRs of energy significantly smaller 
than 1 GeV could be different than that of higher energy particles. We will 
discuss this issue in Sec. 


4.2 Leptonic gamma rays and synchrotron emission 
4.2.1 Inverse Compton scattering and relativistic Bremsstrahlung 


Energetic photons produced in interactions of CR electrons in the ISM also 
contribute to the diffuse gamma-ray emission from the Galactic disk. The two 
dominant radiation mechanisms are relativistic Bremsstrahlung and inverse 
Compton scattering [1970). Unlike hadronic (neutral 
pion decay) emission, which dominates the spectral energy distribution in the 
high energy gamma-ray domain, leptonic processes are characterised by much 
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broader energy spectra, extending down to the soft gamma-ray (MeV) and 
hard X-ray (keV) domains (Fig. (7). 

In inverse Compton scattering, a relativistic electron of mass Mme, Lorentz 
factor y, and energy Ee = ymec? scatters off a soft ambient photon of energy 
e. As a result, the latter is boosted to a much larger energy E. The energy 
of the soft photon seen from the rest frame of the electron can be computed 
performing a Lorentz transformation and is of the order of €e ~ ye. Then, 
provided that ye « mec?, in that rest frame the scattering takes place in the 
Thomson (elastic) regime, and the energy of the scattered photon remains 
e. Lorentz transforming back to the lab frame, one can see that the photon 
energy after the scattering is boosted to Ey ~ ye. In the expressions above, we 
neglected all the (small) factors depending on the angle between the velocity of 
the incident /scattered photons and that of the electron. Averaging over these 
angles one gets the exact result for the typical energy of the scattered photon: 

= (4/3)y%« (Blumenthal and Gould 197). 

In the Thomson regime, a relativistic electron moving through a radiation 
field made of n;agq photons per cubic centimeter would scatter off opnraac 
photons per second, where or is the Thomson cross section. Taking <e) to be 
the mean photon energy and Wrad = Nraa<€ to be the radiation field energy 
density, the rate at which energy is lost by the electron reads: 


dE. 4 4 
rae (OTNradC) (iro) = ZOTCY Wrad (21) 


Pic = 


The energy spectrum of the inverse Compton emission produced by a pop- 
ulation of electrons with a power law spectrum in energy dn./dE.«*cE;° (num- 
ber of particles per unit volume and unit energy) is itself a power law in photon 
energy qi? (E,)cE; $. The expected slope of the gamma-ray spectrum can be 
estimated from the approximate expression: 
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where we implicitly assumed that electrons of Lorentz factor y interacting 
with soft photons of mean energy <e) paer gamma-ray photons with a 
delta function spectrum centred at Ey = (4/3)y?(e). From Eq. [29] it follows 
that s = (6 + 1)/2 (see [Blumenthal and Gould [1970 1970) for a formal derivation 
of this result). 

The interstellar radiation field pervades the Galactic disk and provides a 
background of soft photons for inverse Compton scattering. It is dominated by 
three components: the starlight emission, characterised by an energy density 
in the local ISM equal to wy ~ 0.43 eV/cm® and an average photon energy 
of €, ~ 0.97 eV, the emission from dust heated by starlight, with an energy 
density of wq ~ 0.19 eV/cm® and average photon energy eq ~ 7.8 x 107? eV 
and of course the ubiquitous cosmic microwave Dackeround radiation (CMB), 


whose energy density is wowp ~ 0.26 eV/cm® and average photon energy 
€cmep ~ 6.3x 1074 eV (e.g./Vernetto Vernetto and Lipari Vernetto and Lipari [2016] 2016| and references therein). 
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If Ee = ymec? is the energy of the relativistic electron, interstellar photons 
y g 


undergoing inverse Compton scattering will be boosted to typical energies 


Ee A € 
B, =5.1( =] (=) MeV . (23) 


Then, electrons of energy in the range ~ 100 MeV-10 GeV (which is not probed 
by direct observations) will upscatter photons in a very broad band spanning 
four decades in energy, and centred around E, ~ 5 MeV, 40 keV, and 3 keV 
for starlight, dust emission, and CMB target photons, respectively. This shows 
that observations in the soft gamma-ray and hard X-ray domains can be used 
to probe the spectrum of CR electrons in the energy region between those 
probed by Voyager and AMS-02. 

Relativistic Bremsstrahlung radiation is produced when a relativistic elec- 
tron of Lorentz factor y is accelerated in the Coulomb field of an atomic nucleus 
of the ISM, that can be taken to be at rest, as thermal velocities are « c. If 
the nucleus belongs to a neutral or partially ionised atom, the screening effect 
of the atomic electrons has to be taken into account. It is convenient to look 
at the interaction from the frame where the incident electron is (initially) at 
rest, and the nucleus moves at a speed v ~ c towards the electron. The elec- 
tron feels the Lorentz-transformed electric potential of the incoming nucleus 
as a time varying electric field, which is strong for a time ~ b/yv, where b is 
the impact parameter of the collision. As this time interval is very short, the 
field of the nucleus appears to the electron as a pulse of radiation: a virtual 
quantum. Relativistic Bremsstrahlung can be seen as Compton scattering of 
such virtual quanta. The energy of the emitted photon can be obtained after 
Lorentz-transforming back to the lab frame. 

In relativistic Bremsstrahlung, emitted photons can have any energy Ey 
smaller than the electron energy Fe, with a probability that scales as oga1/E, 
for E, « Ee. This fact has two important consequences. First of all, the 
average energy of the radiated photons is of the order of the electron energy 
E, < Ee, ie. Bremsstrahlung is a catastrophic process, where the electron 
typically loses a significant fraction of its energy in a single encounter with a 
nucleus. Second, the energy spectrum of the emitted photons mimics to that of 
relativistic electrons, as can be easily seen for power law spectra of electrons 
of slope 6. The Bremsstrahlung emissivity is then: q? x Se, dyop pes Dd 


(Blumenthal and Gould] [1970 2004). 


CR electrons with energies in the range ~ 100 MeV-10 GeV produce 
Bremsstrahlung photons having approximatively the same energy. However, in 
this energy range neutral-pion-decay gamma rays dominate the diffuse emis- 
sion (Fig. [7p. It follows that Bremsstrahlung emission is a much less effective 
probe of the CR electron spectrum at energies which are not covered by direct 
observations. 

The expected contribution to the gamma-ray diffuse emission from the disk 
from the two leptonic processes described above is shown in Fig. |7| together 
with the hadronic contribution (red line), that was discussed already in the 
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previous Section. The predictions for the diffuse inverse Compton (green) and 
Bremsstrahlung (cyan line) emissions have been computed from a model of 
CR injection and transport in the entire Galaxy (2011). When com- 
paring predictions with data, one should keep in mind that the former did 
not include line (blue data points) and positronium emission, nor hard X-ray 
sources that dominate the signal below 100 keV. It can be seen that the inverse 
Compton emission largely dominates the diffuse flux below photon energies of 
few MeV, while relativistic Bremsstrahlung is subdominant at all energies (see 
also [Strong etal 200). 

For photon energies in the ~ 1-10 MeV range, expectations fail to fit 
COMPTEL measurements (green datapoint). The observed excess might be 
due to unresolved discrete sources (2011), but to test this hypothesis 
observations by superior sensitivity instruments would be badly needed (re- 
member that COMPTEL data are more than 20 years old). With this respect, 
the launch of the Compton Spectrometer and Imager (COSI, see{Tomsick et al.| 
(2021), expected in 2025, and the ongoing discussion about future mid-scale 
missions operating in the MeV domain will hopefully improve the situation in 


a not too distant future. Proposals include, among others, ee ASTROGAM (de 
Angelis et al.||2018), AMEGO (McEnery et al.{ |2019), and GECCO (Orlando 
et al. 2021). 


The large number of components (continuum, lines, positronium, sources, 
and the COMPTEL excess) contributing to the diffuse hard X-ray /soft gamma- 
ray emission makes it difficult to extract in a reliable way the CR electron 
spectrum from data. Luckily, electrons also interact with the interstellar mag- 
netic field to generate synchrotron photons. The observation of this additional 
radiation component constrains in a more straightforward way the leptonic 
spectrum of CRs. 


4.2.2 Synchrotron emission 


Synchrotron radiation is produced when a relativistic electron is accelerated 
due to the presence of a magnetic field. A relativistic electron of charge e and 
Lorentz factor y moving in a uniform magnetic field of strength B with pitch 
angle a (the angle between the electron velocity and the field) will spiral along 
the magnetic field with a characteristic gyration frequency wy = eB/ymec, 
independent on a. In the non-relativistic regime, the spectrum of the radiation 
emitted by the electron is close to a delta function at the cyclotron frequency 
Ve = W,/2n = eB/2nxm-c. For energetic particles, relativistic corrections apply: 
the spectrum of the emitted radiation is much broader than a delta function, 


and extends up to a characteristic frequency (Blumenthal and Gould} |1970): 


3a Ee \°/ B 
s=- A z4 — | MH 24 
Vs = 57 vesina s(a) Gra) Zz (24) 


4 which does depend on the pitch angle. For LECR electrons the synchrotron 
emission falls in the radio domain. The rightmost equality in the equation 
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above has been computed for the special case a = 7/2 and normalised to the 
typical value of the interstellar magnetic field 2001). Note that the 
energy of the synchrotron photon scales as the square of the electron energy, 
hv,cE?, as for inverse Compton scattering (see Eq. [23). 

The similarity with inverse Compton scattering holds also for the expres- 
sion of the power radiated by an electron: 

Pon = E = soreywp (25) 
which is identical to Eq. except for the fact that the energy density of the 
radiation field w;aq has been substituted with that of the magnetic field wg = 
B? /87. It follows, then, that the spectrum of synchrotron photons produced by 
relativistic electrons following a power law distribution in energy (or frequency, 


as this quantity is customarily used in radio astronomy) dne/dy = Key? 
is also a power law: v x q(v)cK_ BOtD/Ay-O-D/?2 = p-8 
1970}. 

A consequence of Eq. [24]is that observations in the MHz frequency range 
can be used to constrain the CR electron spectrum in the gap between Voy- 
ager and AMS-02 data. We will proceed in a similar fashion to what done in 
Sec. [4.1.1]and discuss first the observations of the diffuse synchrotron emission 
at large Galactic latitudes (b >10°), that can be used to constrain the local 
spectrum of CR electrons. Then, we will discuss the difficulties encountered in 
constraining the spatial distribution of CR electrons throughout the Galactic 
disk. 

The average spectral index of the high latitude diffuse synchrotron emission 
has been derived from various surveys performed between 22 and 408 MHz (see 
the compilation of data in [Strong et al. (2011). The values found in this way 
are all close to 8 ~ 0.5, which corresponds to a slope of the electron spectrum 
equal to 6 ~ 2 at ~ GeV energies. The spectrum of the diffuse synchrotron 
emission steepens at larger frequencies (see e.g. (2008), reaching 
8 ~ 1 at tens of GHz, though with significant scatter around this value 
and references therein). Note that 8 ~ 1 would correspond to a 
slope of the electron spectrum of 6 ~ 3, consistent with the one measured from 
direct (AMS-02) observations of CR electrons. 

constrained the local (~ 1 kpc) interstellar spectrum of CR 
electrons by combining observations of the diffuse radio and gamma-ray emis- 
sion from intermediate Galactic latitudes with direct observations (Voyager 
and AMS-02). The spectrum she obtained is shown as a red dotted line in the 
right panel of Fig.|8} together with the fitting formula we adopted in Eq. 
shown as a solid black line. The spectrum has a doubly broken power law 
shape and is mainly constrained by Voyager and AMS-02 direct observations 
at low (below tens of MeV) and high (above tens of GeV), respectively, and 
by radio observations in the intermediate region, where the spectral energy 
distribution is roughly flat, i.e. E2j¢,;5(Ee) ~ const. 

Finally, measuring the large scale spatial variation of the intensity of CR 
electrons throughout the Galactic disk is not a straightforward task. In prin- 
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ciple, observations of MCs in the gamma-ray and radio domain could be used 
to constrain the CR electron spectrum at specific locations in the disk (in a 
similar fashion to what done in Sec. [4.1.2] for CR protons). Unfortunately, this 
is problematic, for at least four reasons. 


1. As for the diffuse Galactic emission, the leptonic contribution to the gamma- 
ray flux from MCs is expected to dominate in the MeV domain (e.g. 
which is, to date, very poorly explored. 

2. The expected synchrotron radio fluxes from clouds are quite weak, and in 
several cases only upper limits have been obtained after subtracting the 
dominant thermal component 2011). Moreover, in the 
cases where some non-thermal emission has been detected, it was associ- 
ated to small volumes of the cloud, and its origin was likely unrelated to 


interstellar CRs (e.g. in situ acceleration of Soe ean 
of discrete non-thermal or foreground sources, [Rodriguez and Zapata|2013} 
but see also|Jones]2014). A possibly unique exception to this is the central 
molecular zone, a very massive complex of clouds located in the Galac- 
tic centre region, exhibiting a radio synchrotron spectrum 
which may be consistent with the emission of interstellar CR 
electrons penetrating the dense molecular region (2021). 

3. The interactions of CR protons with the gas produce not only neutral pions 
(Eq. [16), but also charged ones. This leads to the creation of electrons and 
positrons (hereafter referred to as secondary electrons) through the decay 
chain m= — pt — e+ (2006). Such secondary electrons are 
produced in large quantities in dense MCs, and the soft gamma-ray and 
radio emission that they generate may compete with that from interstel- 


lar primary electrons (Brown and Marscher\ |1977} |Marscher and Brown 
1978} [Crocker et al.||2007] |Protheroe et al.}|2008). Ultimately, the predom- 


inance of the radiation from primary or secondary electrons depends on 
the transport properties of these particles and on the cloud density. The 
former regulate the capability of energetic particles to penetrate into or 
escape from MCs, while the latter determines how many secondary elec- 
trons are produced per parent CR proton. Recent studies seem to indicate 
that, in fact, the contribution from secondary electrons to the synchrotron 
emission is subdominant (2018). Given 
the current observational difficulties, this issue will be hopefully settled by 
future observations by the Square Kilometer Array 
pois). 


4. The synchrotron emissivity depends on both CR electron density and 
strength of the cloud magnetic field. The former can be extracted from 
observations provided that an independent estimate of the field strength 
is available. Therefore, the uncertainty in the determination of the field 


strength (2012) affects the estimate of the CR electron intensity. 


In the absence of any reliable method to measure the large scale spatial 
distribution of CR electrons, we resort to theoretical predictions (Strong et al.| 


2011} [Orlando and Strong} (2013; |Orlando} (2018; [Di Bernardo et al. |2013). 
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Propagation codes describing the transport of CRs in the entire Galaxy can 
be used to this purpose, as their predictions can be adjusted in order to make 
them consistent with all available data (both direct and indirect) on CRs. 
In this way, the spatial distribution of CR electrons can be predicted and, 
though results depend on the assumptions made on the spatial distribution of 
CR sources, they show that large scale fluctuations in the CR electron density 
are limited to a factor of few (see e.g. Fig. 1 in{Di Bernardo et al.| (2015). 

We conclude that, as for CR protons, large scale spatial variations of CR 
electrons in the Galactic disk are likely to be quite mild. 


5 Integral constraints on the remote intensity of low energy cosmic 
rays: the cosmic-ray ionisation rate 


In this Section we describe a number of observations that can be used to 
measure the intensity of CRs in remote locations in the Galaxy, thereby con- 
straining their spatial distribution in the Galactic disk. The method in based 
on the measurement of the rate of ionisation of the dense interstellar gas, which 
under certain conditions is entirely due to CR ionisation, and it differs from 
that described in the previous Section, as it provides an integral constrain on 
the intensity of CRs, rather than their spectral energy distribution. 


5.1 Cosmic ray induced astrochemistry in interstellar clouds: formation and 
destruction of H} 


The most important photoionising agent in the ISM is UV radiation emitted 
by hot stars. However, UV photons of energy exceeding 13.6 eV (the ionisation 
potential of the very abundant hydrogen) are very quickly absorbed and do 
not pervade the interstellar gas. On the other hand, photons of energy smaller 
than 13.6 eV can penetrate significant column densities of gas, where they 
primarily ionise carbon, which is the most abundant species after hydrogen 
having a smaller ionisation potential (~ 11.3 eV). When a column density of 
gas larger than Ny ~ 8 x 1071 cm~? is traversed, the UV radiation is strongly 
attenuated and LECRs become the dominant ionising agents [1989). 

Following|Snow and McCall] (2006) (see their Fig. 1), we will call diffuse and 


dense MCs those characterised by a column density smaller than ~ 2 x 107! 
cm~? and larger than ~ 8 x 10?! cm~?, respectively. Clouds of intermediate 
column densities are called translucent, as they can still be studied by means of 
absorption spectroscopy of background stars (1989). 

In diffuse clouds, fully exposed to stellar radiation, carbon is predominantly 
singly ionised (C*), and its photoionisation constitutes the main source of free 
electrons. The typical ionisation fraction in these objects is of the order of 
~ 107+. Diffuse clouds are called atomic or molecular depending on the pre- 
dominant state in which hydrogen is found, with atomic hydrogen becoming 
less and less abundant with increasing column density. In translucent clouds, 
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mainly made of molecular hydrogen, photoionisation becomes less important 
and carbon is found primarily in atomic form (C). As a consequence, the elec- 
tron fraction decreases as the gas column density increases. Finally, in dense 
MCs carbon is predominantly locked up in carbon monoxide (CO), the second 
most abundant interstellar molecule after H2, and the ionisation fraction drops 
to ~ 107”. This classification of clouds also serves to describe the different lay- 
ers of a single object, i.e., dense molecular cores will be surrounded by regions 
of lower density having the same properties of translucent clouds, which are in 
turn embedded in envelopes of diffuse gas where hydrogen is present in both 
atomic and molecular form. 

A key point is that molecular hydrogen is primarily ionised by LECRs, 
and not by stellar photons. This is because its ionisation potential is 15.4 eV, 
larger than that of atomic hydrogen. Hg is therefore well protected from ionis- 
ing radiation, which is severely absorbed in the outer atomic layers of clouds. 
Remarkably, the ionisation of Hə triggers the formation of protonated hydro- 
gen, H4, which in turn starts a chain of chemical reactions that leads to the 
formation of a great variety of molecules. The importance of CR ionisation in 
interstellar chemistry was recognised in a number of pioneering papers 
1973), showing 
how the measurements of molecular abundances in clouds could lead to esti- 
mations of the ionisation rate and therefore to constraints on the spectrum of 
the LECRs responsible for the ionisation. 

The goal of this Section is to review briefly the available measurements 
of the CR ionisation rate based on observations of molecular lines from the 
direction of interstellar clouds. For more extended discussions the reader is 


refered to [Tielens] (2013), (2012), [Oka (2006), and 


(2006). A useful database of astrochemical reaction rates can be found in 
Let us start considering the ionisation of molecular hydrogen by a CR 
particle (either a nucleus or an electron): 


Hə + CR — H} +e +CR (26) 


The number of ionisations per unit second and unit volume of gas is (3?,n(Hy), 
where n(Hz) is the number density of molecular hydrogen in the cloud, and 

a is the CR ionisation rate per Hə molecule, that depends only on the 
intensity and spectrum of CR nuclei and electrons. The ionisation is followed 


by the very fast ion-neutral reaction: 
H} + H — H4 +H (27) 


that leads to the formation of protonated hydrogen, Hł , a pivotal molecule in 
interstellar chemistry. As the second reaction (Eq. is orders of magnitude 
faster than the first one (Eq. [26) for any plausible CR intensity and spectrum, 
the latter fixes the production rate of protonated hydrogen, which is then: 
C2.n(Hy). At equilibrium, the abundance of Hj in a cloud is obtained by 


8 KIDA: Kinetic Database for Astrochemistry: |nttps://kida.astrochem-tools.org/ 
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balancing its formation and destruction rates. The destruction of H} proceeds 
in a different way in diffuse and dense clouds. 

In diffuse clouds, the electron fraction is of the order of £e = n(e~)/ny ~ 
n(Ct)/ny, where ny = n(H) + 2n(H2). The abundance of Ct and of both 
atomic and molecular hydrogen have been measured directly for a limited 
number of clouds thanks to the observations of UV absorption lines in the 


spectra of background hot stars (Cardelli et al. |1996; |Bohlin et al.| |1978 
Savage et al.||1977). Based on these measurements, a reference value of x. = 
1074 


is commonly accepted. Remarkably, this value is large enough to make 
dissociative recombination the dominant channel for the destruction of A; 


Hj +e —>+H+H+H or H+H (28) 


The recombination rate per unit volume is then ken(Hz)n(e~), where ke is 
the rate coefficient. 
Balancing formation and destruction rates provides the equilibrium density 


of He, which reads: 
Hə n 2 
an- (8) p a 


Being interested here in the study of LECRs, it is more convenient to rewrite 
the expression above to obtain the CR ionisation rate as a function of the 
other physical quantities: 


Hə 2 kete 7 2 kheLe N(Hf) 


= + 

CR fa, n(Hz) ~ Ta, T (30) 
where we introduced the electron fraction ze, the fraction of hydrogen in molec- 
ular form fy, = 2n(H2)/ny, the HZ column density N(Hj), and the depth 
of the cloud along the line of sight L. Note that the latter quantity is dif- 
ficult to be constrained observationally. Imposing the approximate equality 
n(H}) ~ N(H3)/L, we assumed that Hj is uniformly distributed in the en- 
tire cloud. Making an assumption on the spatial distribution of Ae inside the 
cloud is needed as column densities, rather than volume densities, are more 
easily derived from observations. 

At this point one should notice that all the physical quantities on the right 
side of Eq. [30] could be determined from astronomical observations, except for 
the coefficient ke, which has to be measured in the lab. Unfortunately, for a 
very long time the actual value of ke was so uncertain that its determination 
was considered one of the most controversial issues in the field of electron- 
ion recombination, from both an experimental and theoretical point of view 
(2003). The breakthrough came when it became possible 
to perform laboratory measurements of this quantity under nearly interstel- 
lar conditions, providing a value of ke of the order of few times 1077 cm3/s 
for typical cloud temperatures T of the order of tens of degrees, with a mild 


~ T~°5 dependence (McCall et al.| |2003! |2004). Note that, as these mea- 


surements revealed that dissociative recombination is faster than previously 
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thought, earlier estimates of the CR ionisation rate in diffuse clouds are often 
too low (see for an historical perspective). The parameter- 
isation currently adopted in the KIDA database for the rate coefficient is 
ke = 6.7 x 1078 (7/200 K) ea" 7: 

A very similar line of reasoning can be pushed forward also for the case of 
dense MCs (fy, ~ 1), with one very important difference. Dense clouds are 
well shielded from interstellar radiation and therefore the electron fraction is 
so small (xe ~ 1077) that dissociative recombination of Ht becomes unimpor- 
tant. Due to the large abundance of the CO molecule, Hj is mainly destroyed 
through the proton hop reaction: 


Hi + CO — HCO* + H2 (31) 
characterised by a rate coefficient kco ~ 2 x 107°? cm3/s which is about 


two orders of magnitude slower than ke 2006). Imposing equilibrium 
between formation and destruction, an expression very similar to Eq. [29]can 


be obtained: 
{Gi n(Hə) 
— Co) col aa 


or, passing to column densities: 


CER = 2 kco [=] n(H3) ~ 2 kco | 


woj 


nH L 


which is analog to Eq. A typical value for n(CO)/ny in dense clouds is 
z 1074 1996), similar to the value of x. ~ n(Ct)/ny found in 
diffuse clouds. It follows, then, from Equations [80] and [33] that the density of 
Hj in interstellar clouds (both diffuse and dense) is mainly set by the CR 
ionisation rate G, 

Reversing this argument, one can see how the determination of the column 
density of H} is an essential step towards the determination of the CR ioni- 
sation rate, as GER LN (H4). L is in general unknown, but can be estimated 
when both the column density (Ny) and the average volume density (fix) of 
hydrogen are known, to give L ~ Ny /ii. It seems therefore useful to review, 
at this point, the available observations of He column densities in interstellar 
clouds. 


5.1.1 HZ in diffuse clouds 


The presence of He in interstellar clouds can be revealed in the infrared ab- 
sorption spectra of background stars. Under interstellar conditions, only the 
two lowest rotational levels of H} are expected to be significantly populated. 
Therefore, observations probing the transitions arising from these states can 
be used to probe the entire content of Hf along the line of sight (see e.g. 


Fig. 1 inJMcCall et al.) |1999). In particular, the R(1,0) and R(1,1)" orto-para 
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doublet can be conveniently observed in the IR band, near 3.67 pm P| After 
a long history of attempts, H4 was finally detected in interstellar space by 
o E caview ou I 
astronomy). 

As said above, the interpretation of early measurements of He column den- 
sities suffered from the large uncertainty in the knowledge of the dissociative 
recombination coefficient ke. After measuring this coefficient in the lab, [McCall] 
could estimate with unprecedented accuracy the CR ionisation 
rate in a diffuse cloud towards the star ¢ Persei. The column density of H} 
was derived from the observations of the R(1,0) and R(1,1)“ absorption lines, 
resulting in N(H}) ~ 8 x 1013 cm~?. Moreover, a gas temperature of 23 K 
was estimated from the ratio of the two lines. 

Remarkably, not only N(Hj), but all of the quantities appearing in the 
right hand side of Eq. [30] have been measured along the line of sight towards 
¢ Persei. Absorption lines in the UV domain provided us with an estimate of 
the column density of molecular hydrogen, N(H2) ~ 4.7 x 107° cm~? (Savage! 
fet al. {1977}, atomic hydrogen, N(H) ~ 6.3 x 102° cm~? (Bohlin et al.|/1978), 
and ionised carbon, N(C*) ~ 1.8 x 1017 cm~? (Cardelli et al.| |1996). Such 
values yield a total hydrogen column density Ny = N(H) + 2N(H2) ~ 1.6 x 
10?! cm~?, a fraction of hydrogen in molecular form fy, = 2N(H2)/Nu ~ 0.6, 
and an electron fraction ze ~ N(C+)/Ny ~ 107+. Substituting these values in 
Eq. and using the appropriate value for the dissociative recombination rate 
coefficient ke ~ 2.5 x 1077 cm/s (from KIDA), gives (32,L ~ 7 x 10° cm/s. 

A measurement of the average hydrogen volume density ny is needed to es- 
timate the absorption path length L. This can be obtained from the rotational 
excitation lines of the molecules CO and Cə (1986). 
The method is quite uncertain and gives values of ny of the order of hundreds 
atoms per cubic centimeter. Adopting a reference value of ñy ~ 250 cm~? 
yields a path length L = Ny/ñy ~ 2 pe from which one finally gets a CR 
ionisation rate equal to C42, ~ 10715 s7! [2003). This was a 
quite surprising result because, as we will see in the following, such a value 
for the CR ionisation rate is much larger than that derived from both earlier 
theoretical predictions and measurements of the HF column density in dense 
clouds. 

Triggered by this unexpected finding, surveys of H} in diffuse clouds were 
performed. In about a decade, the number of sight lines explored increased 


from 7 (McCall et al.||2002) to 50, resulting in 21 detections (Indriolo et al.| 
2007; |Indriolo and McCall} |2012). The CR ionisation rates derived in this way 


are shown in Fig.|10jas filled green square data points, together with 30 upper 
limits (green arrows). For some lines of sight, not all of the physical quantities 
needed to compute an from Eq. could be measured, and typical values 
where assumed! After performing a simple statistical analysis, 


McCall) (2012) concluded that the mean value of the CR ionisation rate in 
9 For a detailed description of Hy spectroscopy notation see|Lindsay and McCall) (2001). 


10 Assumed typical values are: ny = 200 cm~?, x3 = 1.5 x 1074, fu. = 0.67, T = 70 K. 
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Fig. 10 Measurements of the cosmic ray ionisation rate in diffuse and dense clouds. Data are 


from (2015) 
Indriolo and McCall] [2012 
open green points), (1998) (filled red circles), [Sabatini et al.[(2020) (filled red 
triangles), and (2007), [Hezarch ot al.| (2008 (2014), 
[Fuente et al] (open red points). The ionisation rate in the local interstellar medium 


is shown as a shaded yellow region for CR protons and with magenta lines for CR electrons 
(dotted and solid line refer to Emin = 3 MeV and 3 keV, respectively). 


diffuse clouds is ~ 3.5 x 10716 s~!, and that there is a probability of 68.3% 
(Lo equivalent) to find a diffuse cloud characterised by a ionisation rate in the 
range spanning from 5.0 x 10717 s~! to 8.8 x 10716 s71, In order to investigate 
the origin of this quite large scattering, correlations between the CR ionisation 
rate and the position in the sky were searched for, to no avail. This means that 
the origin of the scattering is most likely not related to spatial variations of 
the CR intensity on large spatial scales. However, it should be kept in mind 
that all of the clouds in the sample are nearby (background stars are all within 
a distance of ~ 4 kpc, and 60% of them are within 1 kpc), suggesting that 
variations of the CR intensity on quite small spatial scales might occur. 


As we will discuss in Sec. the CR ionisation rate has been predicted 


to decrease with the cloud column density (e.g. (2009). This 
is because LECRs, which are most effective at ionisation, will lose all their 
energy due to ionisation losses while traversing large column densities of gas 
(c.g.[Cravens and Dalgarno}[1978). This would make them less abundant inside 
clouds of larger column density, which would be therefore characterised by a 
smaller average ionisation rate. However, such a trend is not observed in diffuse 
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clouds, as no correlation between ce and N(H2) has been found. The lack 
of a correlation is indeed expected if LECRs permeate the entire volume of 
diffuse clouds, which might be a plausible hypothesis, given the characteristic 
low values of the hydrogen column densities of diffuse clouds. 


5.1.2 HZ in dense clouds 


Hi was first detected, and searched for, in dense clouds (Geballe and Oka 


1996) and, to date, the number of HF detections in dense clouds amounts 
to < 10 (McCall et al [1909 [Brittain et al) 2004) [Gibb et al] 2010). The 
reason why dense clouds were (erroneously) thought to be better targets for 
Hj searches can be understood by comparing Equations and As the 
ratios n(e—)/n(H2z) and n(CO)/n(H2) are, in most cases, of the same order 
(O(10-*)), it follows that, for a fixed value of the CR ionisation rate Caa the 
expected density of Hy in dense clouds is a factor of ke/kco ~ 100 larger than 
that in diffuse clouds (2006). In fact, as we saw in the previous Section, 
this is not the case, as observations demonstrated that HJ is abundant also in 
diffuse clouds. An immediate and important implication emerging from this 
observational finding is that the CR ionisation rate must be different in diffuse 
and dense clouds, and not a standard quantity, as previously assumed. 

In their pioneer paper, presented the first de- 
tection of interstellar H} from two lines of sight towards the young stellar 
objects GL2136 and W33A, both deeply embedded in dense MCs. Further in- 
vestigations on these two objects were performed by (1999), 
who found Hj column densities equal to N(H}) ~ 3.8 x 10% cm~? and 
~ 5.2 x 10'4 cm”, respectively. They also provided estimated temperatures 
equal to 47 and 36 K, inferred from HF line ratios. For such temperatures, the 
KIDA database suggests a value for the rate coefficient that appears in Eq. 
equal to kco ~ 2.1 x 107° cm/s. Assuming a typical abundance ratio for 
dense clouds n(CO)/n(H2) ~ 1.5 x 1074 [Lee et al, [1994 yields CALL ~ 120 
and 160 cm/s for the two lines of sight. 

The clouds’ hydrogen column density can be derived from the measured 
visual extinction and the dust-to-gas conversion factor, giving in Ny ~ 1.8 x 


1073 cm7? and ~ 2.8x 102% cm7? (McCall et al.||/1999). Taking a typical density 


ñy ~ 10° cm-%, as appropriate for the envelopes of deeply embedded stars 
(c.g. [van der Tak et al.|{2000), one obtains path lengths equal to L = Ny/ñy ~ 
0.6 and 0.9 pc, and values of the CR ionisation rate equal to ee ~ 7 and 
9x 107!" s7! for GL2136 and W33A, respectively. These rates are smaller than 
the mean value of ~ 3.5 x 10716 s71 found by [Indriolo and McCall] in 
diffuse clouds (see previous section). What we might be witnessing, here, is the 
predicted decrease of the CR ionisation rate with the cloud column density, due 
to ionisation losses suffered by LECRs while penetrating interstellar clouds. 
In fact, the difference between the values of the CR ionisation rate in dense 
and diffuse clouds might even be larger, as pointed out by[van der Tak and van] 


(2000). This is because a significant fraction of the H} observed along 
the line of sight might reside in the diffuse envelope surrounding the dense 
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cloud cores discussed here, or in foreground diffuse clouds that just happen to 
be located along the line of sight. In support of the latter hypothesis, 
provided evidence for a correlation between the 
measured H3 column density and the cloud heliocentric distance (see their 
Fig. 3), which points to an important role played by intervening clouds. To 
overcome this problem, they combined the observations of Hj with those of the 
HCO? line, which is observed in emission in the millimeter domain. HCOT 
is produced in the dense molecular gas only, and this strongly reduces the 
importance of intervening foreground diffuse clouds. Finally, its abundance is 
proportional to that of H} (see Eq. and therefore to the CR ionisation 
rate. The revised estimates of ee obtained in this way are shown as empty 
green squares in Fig. [10] and were found to be a factor of several smaller than 
those based on Hf observations only. 

As we will see in the next Section, alternative molecular tracers can be used 
to infer in an independent way the value of ce in dense clouds, confirming 
that its typical value is significantly smaller than that found in diffuse clouds. 


5.2 Other molecular tracers of the cosmic ray ionisation rate 


Before the detection of HF in space, constraints on the CR ionisation rate were 
usually obtained from the abundances of molecules like OH and HD (Black 
1906) 

Oxygen chemistry begins with the CR ionisation of atomic hydrogen, fol- 
lowed by charge exchange with oxygen to produce O* in a slightly endothermic 
reaction. This, in turn, starts a chain of exothermic reactions with Hə leading 
to the production of a number of molecular ions (hydrides, made of one sin- 
gle heavy atom and one or more hydrogen), according to the scheme: OF > 
OH* — H,0+t — H30°. In MC interiors, where molecular hydrogen dom- 
inates over atomic one, a very similar chain of reactions is initiated by H$, 
that can interact with oxygen to give either OH+ or HOt (Hollenbach et al.| 
(2012). Finally, OH is produced as a result of the dissociative recombination of 
H307. 

In an analogue manner, the production of HD is the result of CR ionisation 
of atomic hydrogen, followed by charge exchange first with deuterium, Ht + D 
< H + Dt, and then with molecular hydrogen, D+ + Hz e H+ + HD (Black 


and Dalgarno\|1973). If LECRs, rather than UV/X-ray photons, dominate the 


ionisation of hydrogen (as it is the case, unless the cloud is irradiated by an 
UV flux that largely exceeds the local interstellar one, [Hollenbach et al.]2012), 
the abundances of OH and HD can be used to constrain the CR ionisation 
rate. Early estimates of the CR ionisation rate in diffuse clouds based on the 
observations of these two molecules fell in the range of a few 10717 s~! (see 
e.g. and [Federman et al.J1996} but see [van Dishoeck] 
E i also larger values). As we saw in the previous 
Section, the detection of Hj in diffuse clouds led to an average value of ta 
about one order of magnitude larger than those early estimates. 
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It was later recognised that ionised atomic hydrogen can be neutralised via 
charge exchange with polycyclic aromatic hydrocarbons (PAHs}|"] or small 
grains. Therefore, the chain of reactions leading to OH and HD is “leaky”, 
as not all ionisations of H lead to the formation of one of such molecules. 
If this effect is not taken into account while interpreting the abundances of 
OH and HD, the CR ionisation rate is underestimated. This may explain the 
discrepancy with the ionisation rates derived from He observations 
por). 

The opposite might happen deeper into the cloud, when the CR ionisa- 
tion of molecular hydrogen starts the chemistry. There, PAHs can efficiently 
capture electrons, reducing their number and hence diminishing their role in 
the destruction of Hy via dissociative recombination. Then, the abundance 
of Hf and of the an bearing molecular ions increases in the presence of 
PAHs 1.} 2012). In this case, ignoring the presence of PAHs 
when interpreting the abundance of H} may lead to overestimating the CR 
ionisation rate, even though an accurate assessment of this effect is prevented 
due to the uncertainties in our knowledge of the PAH abundance in clouds 


(Shaw and Ferland| |2021). 


5.2.1 Hydride ions: OH*, H2O+, and ArH* 


Not only OH, but also the oxygen-bearing molecular ions formed in the chain of 
reactions that lead to its formation can be studied in order to constrain the CR 
ionisation rate. This fact became of great interest following the detection of the 
three molecular ions OH+, H2O*, and H30? in the submillimeter absorption 
spectra of sources observed Ca the HIFI instrument onboard of the Herschel 
Space Observatory (Germ et al] 2010} [Neufeld et a ROTO). 

In a fully molecular gas ( (ae = 1) each formation of OH™ is followed by 
the formation of H2O* through a a with Hy. Both OH* and H2O* are 
destroyed in interactions with H2, and the two reactions proceed at a similar 
rate. Therefore, the abundance ratio n(OH*)/n(H20+) is expected to be of 
the order of 1 in MCs. On the other hand, in atomic clouds, where fy, is signif- 
icantly smaller than 1, OH* can be destroyed by dissociative recombination 
before interacting with Hz to produce H,O*. As a consequence, the abun- 
dance of the latter molecular ion is suppressed. It follows that the abundance 
ratio n(OH*)/n(H2O+) can be used to estimate the fraction of hydrogen in 
molecular form (Gerin et al. [2012). 

(2015) performed a survey of OH*, H2O*, and H3O* using 
Herschel observations along 20 lines of sight. The first two molecular ions 
were detected in multiple velocity components along every sightline, while the 
third was only detected along 7 sightlines. In the vast majority of cases, the 
n(OH* )/n(H2O0+) abundance ratio was found to be « 1, pointing towards very 
small fractions of hydrogen in molecular form (fH, < 0.1). This fact makes 
Herschel observations unique probes of the CR ionisation rate in atomic clouds. 


11 PAHs are molecules made entirely of carbon and hydrogen arranged in cyclic (ring 


shaped) structures called aromatic rings 2008). 
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Among the oxygen-bearing molecular ions discussed here, OH* is the best 
tracer of CRs, as it is the first one to be produced after the CR ionisation of H, 
and is therefore less affected by “leaks” in the reaction chain. However, as seen 
above, not all CR ionisations of H lead to the production of OH*, but only a 
fraction € of them. Then, the ratio n(OH* )/n(H) constrains the product € (Up, 
and not the two quantities separately. were able to cali- 
brate observationally the efficiency parameter, by using data from a sightline 
were both OH+ and Hj were detected. The CR ionisation rate inferred from 
the H} column density allowed them to break the degeneracy and determine 
€ ~ 0.07, a value which is in broad agreement with theoretical expectations 
[2012h. Note that some of the physical parameters needed 
to determine the CR ionisation rate could not be measured, and typical values 
had to be assumed E2] 

It should be stressed that what is measured, in this case, is the CR ioni- 
sation rate of atomic hydrogen, Ga r, Which is different from that of molecular 
hydrogen, ae However, the two quantities differ by a small factor, as sug- 


gested by the approximate expression derived b (1974): 
1.5 03 = 2.308, which was adopted by (2015). The CR ioni- 
sation rates obtained in this way are shown as filled blue square data points in 
Fig. We plotted only points for which an estimate of the hydrogen column 
density is available, and we do not show measurements relative to lines of sight 
towards the Galactic centre, as this very peculiar region will be examined later 
in this review. All the measurements are characterised by quite low column 
densities of molecular hydrogen, reminding us that hydrides trace LECRs in 
atomic clouds, whose total column density is Ny > N(H2). We decided to plot 
these data versus the molecular hydrogen column density to ease the readabil- 
ity of the Figure (if plotted versus the total hydrogen column density the filled 
blue and green points would significantly overlap). 

Both the mean value and the standard deviation of the CR ionisation rates 


derived from Herschel observations were found to be in very good agreement 
with those derived from H? column densities (see Sec. and |Indriolo and 
McCall||2012). A comparison of the two datasets shows that no significant 


change in the CR ionisation rates is observed in clouds having total hydrogen 


column density in the range Ny ~ 0.7...20 x 10?! cm? (see Fig. 25 in 


Indriolo et al.|2015| where data are plotted versus the total, rather than the 


molecular hydrogen column density). The absence of a correlation between 
= and Ny for diffuse clouds means that they are exposed to the same flux 
of LECRs, regardless of their column density or, in alternative, that what we 
observe is in fact the superposition along the line of sight of many small clouds 
characterised by a very small hydrogen column density. In this latter scenario, 
all the gas probed by OH* or Hf observations would be exposed to the very 
same, almost unattenuated flux of Galactic LECRs. 
Finally, the abundances of OH*, and in some cases of OH too, can also be 
derived from UV absorption spectroscopy, providing alternative estimates of 


12 Assumed typical values: T = 100 K, re = 1.5 x 1074, ny = 35 cm7?. 
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the CR ionisation rate which are in agreement with the Herschel observations 
described above Bo). 

Before moving on to a discussion of dense clouds, it is worth discussing the 
role played by argonium (ArH*) in interstellar chemistry. It was discovered, 
unexpectedly, by Herschel, first as an ubiquitous and unidentified absorption 
feature 2013), later identified by|Barlow et al.](2013). As pointed 
out by |Schilke et al.| (2014), the existence of this noble gas molecular ions at 
detectable levels is the result of an astrochemical conspiracy: i) the ionisation 
potential of argon (15.76 eV) is larger than that of hydrogen, and therefore 
argon is shielded by the interstellar UV radiation and is mostly neutral in 
diffuse clouds; ii) argon is 10 times more easily ionised by LECRs than hydro- 
gen, and this enhances the abundance of Art with respect to that of Ht; iii) 
while other noble gas cations (Het or Net) are neutralised by a dissociative 
ionisation reaction with Hz, such a reaction is endothermic for argon cations, 
thus favouring the creation of argonium through Art + H —> ArHt + H; iv) 
finally, once formed, ArH* can survive for long times as it has an unusually 
slow rate of both dissociative recombination and photodissociation. 

ArH? has been detected from a number of the lines of sight in the survey 
of OH* and H2O* performed by (2015). Using a detailed diffuse 
cloud model, (2017) re-derived the CR ionisation rate by 
fitting simultaneously the column densities of OH+, H2O* and, when available, 
ArH*. To do so, they had to assume that two populations of clouds exist. 
In this scenario, ArHt absorption lines are produced in small atomic clouds 
along the line of sight, characterised by very small molecular fractions (fu, = 
1075...107?), while OH+ and H2O* absorption lines are produced in larger 
column density clouds, characterised by a larger molecular fraction (fu, ~ 
0.2). Small clouds may contribute up to the 50% of the total hydrogen column 
density along a line of sight. Assuming that both populations of clouds are 
irradiated by the very same LECR flux, they derived an average CR ionisation 
rate which is a factor of 2.6 larger than that obtained by (2015). 
This average value and its uncertainty are shown in Fig. |10|as a dashed blue 
line and shaded cyan region, respectively. This fits well with the fact that 
argonium is not protected by large column densities of gas, and therefore the 
CR ionisation rate derived from ArH* observations might be (hopefully?) 
considered a proxy of the actual interstellar value. 


5.2.2 Deuterium fractionation: the HCOt /DCO? ratio 


Deuterium fractionation is the process that enriches the abundance of deu- 
terium in interstellar molecules, leading to ratios between the abundance of 
deuterated molecules and their main isotopologue which are much larger than 
the cosmic elemental ratio D/H ~ 1.6 x 1075. It was soon recognised that 
measurements of the abundance of deuterated molecular ions could be used as 
tools to estimate the ionisation degree in dense clouds (1977), 
which is a crucial step in the determination of the CR ionisation rate. In the 
following, we provide a brief overview of this process, referring the interested 
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reader to the extended reviews by (2002) and [Ceccarelli et al.| (2014) 


and references therein. 

Deuterium fractionation starts with the formation of He a product of 
CR ionisation, followed by an isotope exchange reaction with HD, the major 
reservoir of deuterium in dense clouds: 


H} + HD e HD* + Hə (34) 


The forward reaction proceeds at a rate kp, while the reverse one is endother- 
mic, and is therefore inhibited at the cold temperatures typical of dense clouds. 
As a consequence, the abundance ratio n(H2D+)/n(Hł}) becomes much larger 
than the typical interstellar n(D)/n(H), and its actual level is determined 
by the destruction rate of HD+. Destruction of HəD* mainly proceeds via 
dissociative recombination, at a rate kf (analogue of Eq. or through the 
analogues of Eq. 


HəD* + CO — DCOt + Hy or HCOt + HD (35) 


The total rate of this two reactions is kag ~ kco, the latter being twice 
faster than the former, i.e., the production of DCO? proceeds at a rate kco/3 
1979). 

At this point, it is convenient to introduce the ratio between the abundance 
of DCO* and HCO*, which after some manipulations can be written as: 
_ n(DCOT) 1n(H2D*) 1 kpn(HD) 


~ n(HCO*) 3 n(Ht) ~ 3kin(e) + kgn(CO) 0) 


Rp 


where in the first equality we used the fact that HCOt and DCO® are de- 
stroyed at the same rate due to dissociative recombination, and in the second 
the equilibrium condition between formation and destruction of H2D™, i.e., 
kpn(H3 )n(HD) = kin(H2Dt)n(e7) + kagn(H2D* )n(CO). 

The expression above can be written in a more convenient form introducing 
abundances relative to H, zx = n(X)/ny, to show that the measurement of 
the deuteration fraction of HCO*, Rp, can be used to estimate the electron 


fraction (Caselli) (2002): 


te ~ 1077 E — 1.6 (=. + 5) | (37) 


The numerical values in the expression above were obtained by using the 
typical values for ayp ~ H/D ~ 1.6 x 1075 and zco ~ 1074 and the values of 
the reaction rates from e.g. [Vaupré et al.| (2014) with T = 10 K. The depletion 
factor fa accounts for the fact that in dense clouds a significant fraction of 
CO condenses onto grain surfaces (1/fa is the fraction of CO in gas phase), 
while ô is a correction term to account for the destruction of H2zDt with other 
neutrals (for example atomic oxygen). 

Once the electron fraction is constrained, the CR ionisation rate can be 
derived from the abundance of HCO*t. The equilibrium between production 
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(Eq. and destruction rate (dissociative recombination, at a rate k”) gives 
HCOF = kcon(H})/kn(e), which can be combined with Eq.[32|to obtain the 
observable abundance ratio: 


n(HCO*) kco | CER /n(H) 


CO). aeoe Cl fn) = a 


(38) 


which can be inverted to give the CR ionisation rate (Caselli) |2002): 


cll m ax 10" (18) Ja (545) + (r) ea o) 


Equations [37] and [39] can be used to compute the CR ionisation rate once 
the column densities of HCOt, DCO? and CO are measured from observations 
of emission lines in the millimeter domain, and provided that an estimate of 
the depletion factor is also available. This was done by (1998), 
whose findings are shown in Fig. [10] as filled red circles. 

A refined method has been proposed by (2020), who sug- 
gested to use observations of HD+ (and other deuterated isotopologues of 
H}) together with the DCO* /HCO™ ratio to estimate the column density 
of ‘at, and TE use 5 a [33] to infer ae R. This method has been exploited by 

Sabatini et al] Sabatini et al. , whose results are shown in Fig. |10 [10] as filled red triangles. 

Finally, an ne od circles show the measurement of by Sci i 
Bergin] (2007), [Hezareh et al] (2008), Morales Ortiz ct al.| (2014) and Fuente 
(2016), also based on the observation of various e lines including 
HCO?. 

The main conclusion that can be drawn after inspecting Fig. [10]is that the 
CR ionisation rate seems to stay quite constant, on average, in diffuse clouds 
(both atomic and molecular). The dispersion around the mean value, however, 
is large (at least one order of magnitude). On the other hand, the value of Ga 
found in dense clouds is significantly smaller, suggesting a possible correlation 


with the cloud column density. The interpretation of these trends is postponed 
to Sections [6] and 


5.3 Calculation of the ionisation rate from the cosmic ray spectrum 


The observations of molecular lines from interstellar clouds that we just re- 
viewed are summarised in Fig. where the total rate of ionisations due to 
LECRs per hydrogen molecule Con are plotted versus the cloud column den- 
sity. In fact, several definitions of the CR ionisation rate are possible. In the 
following, we will list them, and show how these rates are connected to the 
spectrum of LECRs. 

Let us start with the total CR ionisation rate per Hz molecule, representing 
the number of ionisation suffered by an hydrogen molecule at a given position 
in space and time (for simplicity the space and time dependences are omitted 
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in the following). This quantity can be written as (e.g. 2009 
2020): 


Emax Emas 


ch=4r]|>, Í dE jp(E)oiR E) [1+0 m(E)] + Í dE jp(E)og$i(E)| (40) 
krs) 0 


where j,(£) is the energy spectrum of CRs (particles per unit surface, solid 
angle, time and kinetic energy) of species k (k = p for CR protons, k = e 
for CR electrons, and k = He, etc., for heavier nuclei), I(H2) = 15.4 eV 
the ionisation potential of Hz and 07°F, () its ionisation cross section. Note 
that throughout this Section, F will represent the total (and not per nucleon) 
kinetic energy of the CR particle. The function ¢(£),H, accounts for the fact 
that electrons produced during the ionisation of Hz molecules can be energetic 
enough to ionise more molecules. The implicit assumptions made in Eq. [40] is 
that such secondary ionisations take place in the same place as the primary 
ones (i.e. the propagation of secondary electrons in space is negligible). The 
second term in Eq. represents the electron capture by a LECR protons 
in molecular hydrogen, characterised by a cross section 7 'jj,(#) (electron 
capture by heavier LECRs is usually neglected). 

Eq. immediately shows the main limitation of the constrains we can 
obtain from the measurement of the CR ionisation rate in interstellar clouds. 
As we deal with an integral quantity, the information on the spectrum of CR 
particles is lost. Also, it is impossible to evaluate from observations how large 
is the contribution from different species (protons, electrons, nuclei) to the 
total ionisation rate. 

The total CR ionisation rate is often expressed in terms of the CR ionisation 
rates due to protons, cae. p and electrons, Cae which are computed by taking 
only the terms k = p or k = e in the summation in Eq. resulting in the 
expression: 

Ch = "Chp + SCR (41) 
where the enhancement factor 7, of the order of ~ 1.5, accounts for the ioni- 
sations from CR nuclei heavier than hydrogen (see the appendix in [Padovani] 
[et al-]2009] for a derivation of n, or[Chabot[2016]for a detailed calculation of the 
ionisation rate of CR heavy nuclei). Finally, as we already saw in Sec. [5.2.1] the 
CR ionisation rate can also be defined per H atom (and not molecule), (Up. 
However, according to the estimate by (1974), the two 
rates are connected through a small multiplication factor, 1.5 Con = 2.3 Ca 
and therefore in the following we describe in some detail the derivation of os 
only. 

The cross sections for electron impact (e~ + Hə > Hj + 2 e7, red curve, 


Kim and Rudd/1994 2000 2021), electron capture (H+ 
+ H2 > Hj + H, dotted blue curve, 2009 2021), 


and proton impact (Ht + Hə > Hy + H + e`, dashed blue curve, 


1988; |Rudd et al.|}1992) on molecular hydrogen are plotted in the left panel 


of Fig. The solid blue curve represents the sum of the electron capture 
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Fig. 11 Left: ionisation cross sections for molecular hydrogen: electron impact (e7 + H2 > 
Hy + 2 e7, red curve), electron capture (Ht + H2 > Hy + H, dotted blue curve), proton 
impact (H+ + H2 > Hy + H + e7, dashed blue curve), and the sum of the last two 
(solid blue curve). Right: number of secondary ionisations per primary ionisation due to 
CR protons (blue) and electrons (red curve). 


and proton impact cross sections. In order to obtain the proton impact cross 
section at relativistic energies, we adjusted the non-relativistic model by 
to match the appropriate high energy scaling (see e.g. {Fano} 1954). 

The impact of the secondary ionisations can be estimated from the dif- 
ferential ionisation cross section doj°fj,(E)/dE, which describes the energy 
spectral distribution of electrons produced in primary ionisations. Here, Æ 
is the energy of the ionising primary particle of species k and Es that of 
the ejected secondary electron. If j,(£) is the intensity of primary parti- 
cles, the differential a rate of secondary electrons per unit volume 
is qs(Es) = 4rn(H2) | j (E) (doje, (E)/dEs)dE, where the integral has to be 
performed for energies E larger than J(H2). The equilibrium spectrum of sec- 
ondary electrons can be approximated as ns(Es) ~ qs(Es)r, where 7 is the 
residence time in the region of interest. If we assume that secondary electrons 
do not travel away from the production site, the residence time reduces to the 
energy loss time Tioss( Es). As the energy loss time is short in dense environ- 
ment, this approximation is in general not too bad (see however 
[2021] for a discussion of the general case). Finally, the ionisation rate (oa sec- 
ondary electrons can be computed as ¢12 k = 4r f| js(Es)o ih, (Es)dEs, where 
we introduced the intensity js(Es) = (vs/4)ns(Es), Us being the velocity of 
an electron of energy Es. Also in this case the integral has to be performed 
for E, > I(H2). At this point, after some manipulations and making use of 


Eq. [40] one obtains 2015): 


1 E dE oii (Es) Es doki, (E) 
IŒH2) L(Es) dE; 


Pk,He (E) = ion 


Tk, i (E) a 
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which depends only on cross sections and on an energy loss function Le( Es) = 
E,/n(H2)vsTioss, that will be introduced and discussed in Sec. [6.2] We antici- 
pate here that the loss function does not depend on the ambient gas density, 
as in dense environments Tigss¢1/n(Hz2). The functions ¢x,4,(£) represent 
the average number of ionisations due to secondary electrons per ionisation of 
Hə by a primary particle of energy E. They are shown in the right panel of 
ne Batic CR protons (k = p) and electrons (k = e). Note that, formally, 
Eq. is appropriate for a gas made uniquely of molecular hydrogen, but 
it provides a reasonably accurate approximation also for typical abundances 
found in clouds. 

At this point, we can attempt a first prediction of the ionisation rates 
expected in clouds, and compare them with observations. The most simplistic 
prediction can be obtained assuming that the local interstellar spectra of CRs 
that we discussed in Sections[3.4]and[4]are representative of the entire Galaxy. 
The ionisation rate can be computed by plugging in Eq. [40|the local interstellar 
spectra of CR protons and electrons, while Eq. [41] can be used to evaluate the 
contribution of heavier CR nuclei. In doing so, one should notice that local 
interstellar spectra of CRs are known only for particle energies larger than 
few MeV ona and (6), while ionisation cross sections peak at much lower 
energies (Fig.|11). Therefore, an extrapolation of LECR spectra to low energy 
might be required. However, as we saw in Sec. [3.4] (see the discussion on the top 
right panel of Fig. (6), performing a spectral gxiTapolationi is mandatory for CR 
electrons, as their contribution to ga Re 18 most likely dominated by sub-MeV 
particles, for which no observations exist. On the other hand, the contribution 
from CR protons to the total ionisation rate is likely dominated by particles in 
the sub-GeV domain, which is well constrained by observations, and therefore 
no extrapolation is needed (2009); (2018) 
other words, if we extrapolate the local spectra of CR protons (Eq. |14) and 
electrons (Eq. |1 [15)down to an energy Emin, the total CR ionisation rate will 
depend on this minimum energy as as ee Emin) ~ EA aT Can, e(Emin)- 

The ionisation rates estimated in this way are reported in Table [5.3] and 
plotted in Fig. as a yellow shaded region for CR. protons, and magenta 
lines for CR electrons. The dotted line refers to the guaranteed ionisation rate 
due to CR electrons, obtained setting Emin = 3 MeV (roughly corresponding 
to the minimum energy at which observations are available), while the solid 
line refers to Emin = 3 keV. The CR ionisation rate for local CR protons is 
at the same level of that found in dense clouds, while in order to match the 
larger ionisation rates observed in diffuse clouds an extrapolation of the CR 
electron spectrum of 3 orders of magnitude in energy (down to ~ 3 keV) is 
required. Such a comparison between observations and expectations should be 
taken with a grain of salt, as we are in fact mixing local expected ionisation 
rates with those measured in remote regions. Moreover, we are neglecting the 
suppression of the CR intensity (and therefore of the ionisation rate) induced 


13 This claim is true unless an unknown and steep spectral component pops up in the sub 
MeV domain. We will briefly discuss this possibility in Sec. 
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Cosmic-ray ionisation rates in the local ISM 


Emin (MeV) | GER p (10717 871) | CER e (10717 s71) | c3 (10717 s71) 
3 2.1 0.49 3.6 
0.3 2.3 T. 4.6 
0.03 2.4 4.6 8.2 
0.003 2.4 48 52 


Table 1 Cosmic-ray ionisation rates in the local interstellar medium obtained integrating 
Eq. [40] above a particle energy Emin. 


by energy losses during the penetration of such particles into MCs, a thing that 
will be discussed in detail in the coming Section. Nevertheless, the preliminary 
comparison between observations and expectations shown in Fig. [10] is useful, 
as it suggests that it might not be easy to explain the large ionisation rates 
observed in clouds, unless quite wild spectral extrapolations are invoked, or 
unless we give up the idea that the intensity of CRs is roughly uniform in the 
Galaxy. 

Before concluding we should remind that secondary electrons produced in 
ionisation events deposit a fraction of their energy into the cloud in form of 
gas heating. In molecular gas, CRs are the main heating agents. The heating 
happens through different channels, the most important being the dissociation 
of molecular hydrogen [2004). If, on average, each electron 
produced in a CR ionisation event releases an amount of heat Q, the gas 
is heated at a rate which is simply proportional to the CR ionisation rate: 
In = 6G2n(H2)Q. Typically, the value of Q is of the order of several eV. The 
balance between heating and cooling determines the thermal properties of the 
gas, and is therefore a crucial ingredient in the understanding of MCs. As we 
will not discuss this topic in this review, we refer the reader to the work by 


Glassgold et al.| (2012) and references therein. 


6 The transport of cosmic rays in and around MCs 


In the previous two Sections and [5) we described how the observations of 
MCs can be used to impose constraints on the intensity and on the spectrum of 
LECRs inside these objects. At low particle energies, the severe energy losses 
suffered by CRs in the dense and neutral gas may prevent them to penetrate 
freely into the clouds interior. As a result, the pristine interstellar spectrum 
of CRs will be modulated by energy losses, and the amount of modulation 
will depend on the properties of the transport of particles into clouds. The 
faster the penetration, the smaller the modulation, as particles might cross 
the cloud before losing significant energy. The fastest penetration is obtained 
assuming that particles penetrate into clouds moving along straight lines. How- 
ever, straight line propagation is not likely to happen, due to the ubiquitous 
presence of turbulent magnetic fields in the ISM, which deflect the trajectories 
of charged particles. On one side, this certainly complicates the interpretation 
of available observations, but on the other hand it provides us with a way to 
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confront with data the theoretical expectations derived from our (quite poor, 
in fact) knowledge of the transport of particles in turbulent fields. 

The transport of CRs in turbulent magnetic fields constitutes by itself 
an entire field of research. Therefore, we will limit ourselves to sketch the 
main aspects of this complex problem, and refer the interested reader to the 


) and |Schlickeiser| (2002), or to the 
recent reviews by (2021) and (2021). 


, let us start by considering the simplest problem 
of a LECR proton moving with a velocity v directed at an angle V (called pitch 
angle) with respect to an uniform magnetic field B. Subject to the Lorentz 
force, the LECR particle will move at a constant speed around magnetic field 
lines, following an helical trajectory of radius rg = p,c/eB (called Larmor 
radius) and gyration frequency 2, = v/rg = eB/ympe. Here, py and v] are 


the components of the particle momentum and velocity orthogonal to B , while 
y is the particle Lorentz factor. 

Perfectly uniform fields are of course an unrealistic idealisation, and there- 
fore we investigate now the effect of the presence of small scale fluctuations 
(MHD waves) superimposed to the uniform background field B. In particular, 
we limit our analysis to slab Alfvén waves, i.e., transverse, small amplitude 
magnetic disturbances (5B 1 Band ôB « B), which propagate along field 
lines at the Alfvén speed v4 = B/./470;, where o; is the ion mass density in 
the ISM, and oscillate with a frequency w ~ kva, k being the wavenumbey!4| 

The fluctuation in the magnetic field will add a y x 6B component to the 
Lorentz force that will perturb the trajectory of the CR particle. In fact, under 
most circumstances the interaction between a CR particle and a wave has very 
little effect, as v x 6B oscillates many times during one gyration of the particle 
around the magnetic field, so that the time integrated effect of the Lorentz force 
on the particle averages to zero. The particle trajectory is impacted only in the 
special case of resonant interactions, that occur when the particle interaction 
with a wave period lasts a time Te = 27/vk which is identical to the gyration 
time Tg = 27/92,. When the condition Te = Tg is satisfied the Lorentz force 
does not change direction during the entire interaction, and changes the pitch 
angle (but not the energy) of the CR particle in a systematic way (with this 
respect, Fig. 3 from [Wentzel] is particularly instructive). 

Neglecting factors of order unity, the resonance condition can be written as 
rg ~ 1/k, indicating that, in first approximation, CRs of a given energy interact 
only with waves of a given frequency. For the sake of clarity, we consider the 
repeated resonant scatterings between a particle and a train of uncorrelated 
wave packets of length 1/k (one period). The effect of one of such wave packets 
is to bend the field line by an angle 6B/B, and the particle pitch angle will 
be changed by this very same amount after scattering. As wave packets are 
uncorrelated (random phases), sometimes the change in the particle pitch angle 
will be positive, and some other times negative. Thus, after a time t the root 
mean square variation of the pitch angle will be of the order of 4/92) ~ 


14 See e.g. [Hollweg] (1974) for a formal approach to MHD waves. 
8 
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VN(6B;/B), where N = t/Te is the number of encounters with wave packets, 
and the subscript k indicates that only fluctuations at the resonant scale have 
to be considered. 

It follows that an initially anisotropic distribution of particles will be 
isotropised after a time given by the equation ./<0?) ~ 1 (a large angle), yield- 
ing t;,5 ~ 2,(5B,/B)?, where small numerical factors have been neglected. 
The isotropisation time tiso is the time it takes a CR particle to lose memory 
of its direction of motion, and the distance traveled by the particle during this 
time, Amfp ~ Uitiso, corresponds then to the mean free path for spatial scat- 
tering along the field line (i.e. after traveling a distance Amfp along a magnetic 
field line, a particle will have the same probability to continue to move forward 
or change direction and go backward). The spatial diffusion coefficient is then: 


4 TgV 


1 
Aag 3r 0B, / B) 


: (43) 


Amfp|| = 


where the subscript || indicates that the coefficient refers to the one-dimensional 
diffusion along magnetic field lines, and the numerical factors come from a 
more detailed calculation (e.g. [Blandford and Eichler} {1987). 

The derivation of the diffusion coefficient provided above, though admit- 
tedly oversimplified!| highlights a crucial aspect of particle transport: the 
energy dependence of the spatial diffusion coefficient is determined by the 
power spectrum of the magnetic perturbations, (6B;/B)?, through the reso- 
nance condition k ~ 1/r,. This quantity represents the amount of magnetic 
energy available in form of fluctuations at the scale k, and it is unfortunately 
very difficult to be determined, both observationally and theoretically. Note 
that, in this notation, the total magnetic energy in form of fluctuations is 
(6Brot/B)? = \(6B;,/B)?dk/k. 

In addition to parallel spatial diffusion, the scattering of CRs onto MHD 
waves induces also a diffusion perpendicular to field lines, which is believed 
to be much less effective than the parallel one in the limit of small amplitude 
fluctuations (e.g. (2001). Particles can also be transported per- 
pendicularly to the mean magnetic field due to the large scale wandering of 
magnetic field lines or to plasma motions (e.g. 1993). 
In his Figures 1 and 8, provided a compendium of the var- 


ious regimes of perpendicular transport, discussing them as a function of few 
parameters that determine the properties of the turbulent ambient magnetic 
field. He also suggested that magnetic mirroring might play a role when CRs 
encounter MCs, as magnetic field lines are focussed into their dense cores. 
On large Galactic scales, the superposition of all the effects mentioned 
above might lead to an effective three-dimensional isotropic spatial diffusion 
of CRs, as it is indeed often assumed in phenomenological studies (e.g. [Strong 
and references therein). However, at the smaller spatial scales 


15 For a detailed derivation of the spatial diffusion coefficient from the spectrum of magnetic 
fluctuations based on a state of the art theory of magnetic turbulence see for example{Fornieri| 


(2021) and the list of references therein. 
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of MCs the situation is probably different. The transport of particles in and 
around clouds is likely to be one-dimensional, as the size of such objects is 
roughly of the same order of the coherence length of the interstellar magnetic 
field. 

Finally, the origin of the magnetic turbulence needed to scatter CRs re- 
mains debated. However, an appealing scenario was proposed in the late 1960s, 
based on a plasma instability driven by the (small) anisotropy of CRs 
[Kulsrud and Pearce} 1969). As a result of the instability, 
CRs transfer momentum to MHD waves, enhancing their amplitude 6B;/B, 
and reducing correspondingly Dj. In other words, the transport of particles 
might be, in fact, a non-linear process, with CRs producing the turbulence 
that in turn scatters them. 


6.1 Nonlinear transport of cosmic rays: streaming instability 


In order to understand how CRs can generate magnetic turbulence, consider 
a slightly anisotropic distribution of energetic particles in a medium where a 
uniform background magnetic field is present. Due to the anisotropy, particles 
will move along the field lines at an average speed vq, called streaming or 
drift velocity. Assuming that energetic particles are protons of typical Lorentz 
factor y and volume density nor, we can compute the initial total momentum 
carried by the ensemble of such particles as P; = noryvam,p. We further 
assume that magnetic fluctuations of very small amplitude are superimposed 
to the background field. These fluctuations will scatter CRs and lead their 
angular distribution towards isotropy. This will reduce both the streaming 
velocity and the momentum carried by CR particles. 

Being a conserved quantity, the momentum lost by CRs has to be trans- 
ferred to something else. It can be shown that, due to resonant scattering, 
the momentum is transferred to resonant waves 2005). Momentum 
conservation also implies that only waves traveling in the same direction of 
the drift velocity will grow (while counterpropagating waves will be damped). 
The streaming velocity decreases until the CR particle distribution function 
becomes isotropic in the frame of the waves. Therefore, in the lab frame the 
final momentum carried by CRs will not be zero, but rather Pp = ncRYVAMp, 
as waves move at a speed vy. 

The rate at which momentum is transferred from CRs to waves is dPor/dt = 
—( P; — P;)/Tiso, where Tiso is the particle isotropisation time derived in the 
previous Section. This has to be equal to the momentum gained by resonant 
waves dP,,/dt = [7,5 Bz/8mv,, which is equal to the momentum of waves times 


a growth rate coefficient T}. Equating the two rates yields (Kulsrud| |2005): 


NCR [Va— VA Ud — VA 
Fy ~ amy See (MPA) = r, (44) (44) 


where a is a numerical factor of order unity, A; the mass number of the dom- 
inant ion in the gas, 2 = 72, is the nonrelativistic cyclotron frequency, and 
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I the growth rate corresponding to vą/va = 2. The term in the parenthe- 
sis shows that the instability is triggered when CRs stream faster than the 
Alfvén speed, and that it self-regulates, as a very large drift velocity vg > v4 
would induce a very fast growth of the waves, that would in turn reduce vq. 
Finally, the exponential growth of waves is very fast, being characterised by an 
e-folding time (6B,,«e’#") which is of the order of 1/I, ~ 5 x 10° yr for typical 
interstellar conditions (B ~ 3 uG, nor/ni ~ 1071, and (va — va)/va ~ 2). 
Eq. [44)has been derived in an approximate way, but is in fact quite accurate. 
The exact expression for a CR particle distribution function equal to a power 
law in momentum, f(p)ocp~°, is recovered after setting a = (7/4)(a—3)/(a—2) 
and substituting nog with ncr(> Pres), i.e. the total number of particles 


per unit volume with momentum exceeding the resonant one (Kulsrud and 
Cesarsky| |1971). An equivalent and perhaps more convenient way to express 
the growth rate was derived by|Skilling}(1975b), in terms of the spatial gradient 


of the CR particle distribution function: 


16r?vavpt of 
3 ôB? ôx 


Ty (45) 


where «x is the coordinate measured along the field line. This can be rewritten 


as: 
2 
2 (SF) pe = (46) 


B Ox 


where we have introduced the CR partial pressure Por = (47/3)up* f /(B?/87) 
normalised to the background field pressure. Eq. shows that the rate at 
which the energy density of Alfvén waves increases (left side) is equal to the 
rate at which CRs do work in scattering off the waves (right side). 

The mechanism for the amplification of magnetic fluctuations described 
here is called streaming instability, and has been recently reviewed by 


cowith et al.| (2021). In most circumstances, the instability proceeds until equi- 


librium is reached with competing mechanisms of wave damping. In this re- 
gard, noticed that Alfvén waves generated by 
the instability may be damped in interactions with waves from interstellar 
turbulence, and that this is likely to limit the viability of this mechanism to 
particle energies below ~ 100 GeV. As this is the range of energy considered 
in this review, we will assume in the following that the transport of CRs is 
regulated by self-generated turbulence. 

Another process of particular interest is the ion-neutral (or ambipolar) 
damping of Alfvén waves, which operates in partially ionised media 
202). Tt is due to the 
fact that ions follow the motion of the field lines more closely than neutrals. 
As a result of their relative motion, momentum-exchanging collisions between 
ions and neutrals take place at a rate Vin ~ Mn/(Mi + Mn)<TVYmeNn, Which 
depends mainly on the mass of the most abundant ion (m;) and neutral species 
(Mn) in the gas, on the number density of neutrals nn, and on the momentum 
transfer coefficient for ion-neutral collisions (avmt. The friction between ions 
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and neutrals subtracts energy to the waves, which are damped. For high fre- 
quency waves (which resonate with low energy particles) the damping rate is 
Tin = Vin /2 (See. Zweibel and Shall] 1952). 

When the growth of waves is balanced by damping (I, = Tin), making use 
of Eq. [44] one gets va/va = 14+ (Lin/I2), which means that CRs can stream 
much faster than the Alfvén waves if Ij, > I>. This can indeed happen in 
MCs, where the density of neutrals is large. An estimate of the damping rate 


for typical conditions found in diffuse cloudd!>] was given in|Pinto and Galli 
(2008) and (2021), who found a very fast damping time of 


1/Lin ~ 1 yr. It follows that streaming instability is inefficient in clouds, and 
that CRs stream almost freely in their interiors, unless their number density 
is of the order of nor ~ 1077 cm—3, which is quite large. 


6.2 Energy losses of cosmic rays in the interstellar medium 


An important aspect of the CR transport into clouds is the fact that LECRs 
lose energy when they interact with the dense molecular gas. In some cases 
(e.g. ionisation losses), the energetic particle loses a small fraction of its en- 
ergy in the interaction, and the process can be described by a continuous 
function called energy loss rate, dE/dt|;, where E is the energy of the incident 
particle and i refers to the interaction responsible for the loss of energy (e.g. 
i = ion for ionisation losses). Other processes, such as proton-proton interac- 
tions (i = pp) or Bremsstrahlung (i = B) are catastrophic, i.e. the energetic 
particle loses a significant fraction of its energy in a single interaction. How- 
ever, also in this case it is convenient to introduce a continuous function that 


has to be considered as an average energy loss rate (Berezinskii et al.{ |1990 
2004). This is called the continuous slowing down approximation 
and has been widely used in the literature (see e.g.|Padovani et al.||2009| and 


references therein). The energy loss rate can be defined in an alternative way 
as dE/da|; = (1/v)dE/dt|; to indicate the amount of energy lost in an in- 
finitesimal displacement of the particle dz = vdt, v being the velocity of the 
energetic particle. 

In the dense and neutral environments considered in this review, LECRs 
lose energy mainly via ionisation losses at low particle energies (both CR 
electrons and nuclei), and via Bremsstrahlung (CR electrons) or proton-proton 
interactions (CR nuclei) at large particle energies. Remarkably, for all of these 
processes the energy loss rate is proportional to the total (all elements) ambient 
gas density, dE/da|,ocn, and is therefore possible to build a quantity, called 
energy loss function, which is independent on the gas density, and reads (e.g. 


Padovani et al} 2009): 
1 /dE dE 
mæ =-= ($=) =- (47) 
k 


n \ dx dN 


16 Assumed values: nn = 300 cm~’, T = 50 K, xe ~ 1074, Hz and C+ dominant neutral 
and ion species, respectively. 
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Fig. 12 Left: average gas column density crossed by a CR proton (blue curve) or elec- 
tron (red curve) before losing its initial energy E. Energy losses are dominated by ionisa- 
tion, except at the highest energies where proton-proton interactions (dotted blue line) and 
Bremsstrahlung (dotted red line) dominate. The yellow shaded region marks the transition 
between diffuse and dense clouds. Right: CR ionisation rate versus cloud column density. 
Data as in Fig. [10] Blue lines show the expected ionisation rates for the free streaming trans- 
port model for CR protons (dotted), electrons (dashed), and the sum of the total, including 
heavy nuclei (solid). 


where k refers to the particle species, (dE/dx); is the sum of all the relevant 
loss rates, and N = Sndz is the total gas column density. From the second 
equality in the expression above, we see that the energy loss function represents 
the energy lost per gas column density traversed by the particle. 


(2018) calculated the energy loss function for CR protons 


and electrons. They considered an ISM composed predominantly by molecular 
hydrogen and containing traces of heavier elements according to typical inter- 
stellar abundances. Such a gas contains a fraction ey, = 0.835 of molecular 
hydrogen, and therefore the column density of molecular hydrogen, which has 
been most often used in the previous Sections, is N(H2) = cu, N. Starting 
from the values of the energy loss function given by es (2018), 
we built the quantity ¢y,£/L,, which has the dimensions of a column den- 
sity of Hə, and provides an approximate estimate of the amount of interstellar 
matter a CR particle can traverse before losing most of its energy. This quan- 
tity is shown in Fig. (left panel) for CR protons (solid blue curve) and 
electrons (solid red curve). For most particle energies, ionisation losses pro- 
vide the dominant contribution to the loss function, with the exception of the 
highest particle energies, where proton-proton interactions (dotted blue curve) 
or Bremsstrahlung (dotted red curve) dominate. 


The shaded yellow region in Fig. [12]marks the (blurred) boundary between 
column densities characteristic of diffuse and dense clouds. The intersection 
between the shaded region and the solid curves shows that, even in the most 
optimistic scenario in which CRs penetrate clouds along straight lines, CR 
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protons (electrons) of energy smaller than ~ 5 MeV (~ 0.1 MeV) would be 
excluded by dense clouds, as they would lose their energy before reaching the 
cloud’s centre. The exclusion of LECRs from cloud interiors suggests a possible 
explanation for the difference in the ionisation rates observed in diffuse and 
dense clouds: electrons characterised by particle energies in the keV domain 
might penetrate into diffuse clouds and produce there a large ionisation rate. 
On the other hand, they might be excluded from dense clouds, which would 
be then characterised by a much lower ionisation rate. The viability of this 
scenario will be discussed extensively in the coming Section. 


6.3 Modelling the penetration of low energy cosmic rays in molecular clouds 


We have now all the elements required to build models describing the transport 
of LECRs into clouds. These models can be used to make predictions of the 
CR ionisation rate as a function of the cloud column density. Unfortunately, 
our poor knowledge of the properties of the transport of CRs in turbulent 
fields will make predictions quite uncertain. For this reason, it is useful to 
consider two extreme scenarios for the propagation of LECRs, to obtain the 
most optimistic and pessimistic predictions on the CR penetration into clouds. 
It is very likely that neither scenario is realistic, but reality will certainly be 
somewhere between these two extremes. 

In the scenario that maximises the penetration of LECRs into clouds, and 
therefore also the CR ionisation rate, LECRs travel along straight lines from 
the edge to the centre of the cloud. We call this mode of propagation the free 
streaming transport, and we refer the reader to (2009), who 
described it in great detail. In a more conservative scenario, the streaming 
of LECRs into clouds triggers the plasma instability described in Sec. 
which in turn strongly suppresses the penetration of CRs as they scatter off 
self-generated Alfvén waves. As a consequence, the expected ionisation rate 
is reduced accordingly. We call this mode of propagation the self-regulating 
diffusive transport. The latter scenario was first suggested in a number of pi- 
oneering papers published in the late seventies/early eighties 


1976} |Cesarsky and Volk} |1978; |Morfill| |1982a]b), and recently recon- 
sidered by various authors (Everett and Zweibel} |2011; [Morlino and Gabici| 
2015} |Schlickeiser et al. |2016} |Ivlev et al.||2018) . 


If CRs travel along straight lines, the left panel of Fig. [12] can be used to 
estimate the minimum energy E% į» of particles (k = p,e for CR protons and 
electrons) that can penetrate into a cloud of column density Ny,. This can be 
done by solving the equation ey, E*,,,,/Ly(E* in) = Nu,. Once E*,,,.(Nu,) is 
derived, an approximate estimate of the CR ionisation rate can be obtained 
by means of Eq.|40| where the lower limit of the integration should be changed 
from I(H2) to E%,,,,(Nu,). In this Section, we will assume that the intensity 
of CR nuclei and electrons just outside of the cloud is identical to the one 
measured in the local ISM. The CR ionisation rate obtained in this way are 
plotted as blue lines in the right panel Fig. for CR protons (dotted), elec- 
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trons (dashed), and all particles, including heavy nuclei (solid line). Note that 
our approximate result is very close to that from the more accurate calculation 
from [Padovani etal (2022) 

The free streaming transport model is an idealisation, but it remains nev- 
ertheless quite useful as it provides a way to estimate the maximum possible 
CR ionisation rate for a given intensity of CRs just outside of a cloud. A more 
realistic description can be obtained by assuming that CRs move at a con- 
stant speed along magnetic field lines, which are in general turbulent (see e.g. 
Figures 7 to 12 from [Padovani et. al.||2013). Under these circumstances, the 
effective gas column density traversed by CRs to reach the cloud’s core would 
be larger than the observed line-of-sight cloud column density, and therefore 
smaller values of the CR ionisation rate would be expected in this case. In 
addition to that, the convergence of magnetic field lines towards cloud cores 
would focus CRs, enhancing their intensity, but, at the same time, some CRs 
would be excluded from the cloud’s core due to magnetic mirroring. In fact, it 
has been shown that the effects of focussing and mirroring cancel each other 
out almost perfectly, with the exception of magnetic pockets (local minima of 
the magnetic field strength), where the intensity of CRs can be significantly 
reduced (Padovani et al 2013) Silsbee oa 2018). 

In order to highlight the main limitation of the free streaming transport 
model, consider a cloud of column density Ny, embedded in a uniform mag- 
netic field, and assume that the flux of CRs of energy E impinging on each 
side of the cloud is FE (E). For E > E£ n, CRs cross the cloud without losing 
energy, and therefore the flux of CRs getting out from each side of the cloud 
will be F£ (E) = FE (E). However, if we consider lower energies, E = Ef in, 
energy losses become important, and the beam of CRs will be a eiaei while 
traversing the cloud. This will result in F* ,(E) < F8 (E), i.e. in an anisotropic 
distribution of CRs at the cloud boundary. The anisotropy will become of or- 
der unity for E < Ef in. As we saw in Sec. anisotropic distributions of 
particles trigger the streaming instability. As a result, the amplitude of mag- 
netic fluctuations increases, and the penetration of CRs into clouds becomes 
a non-linear process. As FE (E) < FE (E), the stream of particles is directed 
towards the cloud, and a converging flow of amplified Alfvén waves will form. 
Note that under most circumstances the growth of waves takes place just out- 
side of the cloud, in the mostly ionised diffuse ISM, and not inside it, where 
streaming instability is inhibited by ion-neutral damping (Skilling and Strong 
but see also for a scenario 
where streaming instability is excited in the diffuse envelope of the cloud). 

The emerging picture is that of a self regulating diffusive penetration of 
LECRs in clouds. The intensity of CRs, assumed to be roughly uniform in 
the Galaxy, is suppressed inside clouds due to severe ionisation energy losses. 
This generates a gradient in the spatial distribution of CRs around clouds, 
that triggers the amplification of Alfvén waves due to streaming instability. 
LECRs are advected towards the cloud by a converging flow of Alfvén waves, 
and at the same time undergo spatial diffusion along the background magnetic 
field due to resonant scattering off waves. Inside clouds, CRs do not diffuse 
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but they rather free stream along magnetic field lines, due to the absence of 
efficient scattering (ion-neutral friction is very effective in damping waves). 

Following/Morfill] (1982a), we simplify the problem by considering a uniform 
magnetic field oriented along the z axis, passing through a cloud of total 
column density N located at z = 0. In this idealised picture, the cloud has no 
thickness along z, i.e., CRs probe its entire column density every time they 
cross the point defined by z = 0. The one-dimensional version of the transport 
equation (Eq. |4) derived in Sec. can be used to study the penetration of 
LECRs in clouds. To do so, the energy loss term —(1/p?)0(p?pf)/ép has to 
be added to its right hand side, to account for ionisation losses in the cloud. 
The momentum loss rate can be written as p = L,.N0d(z), where Ly, is the loss 
function defined in the previous Section (2001). Such a choice 
of the particle momentum loss rate is equivalent to assuming that outside of 
the cloud (z # 0) the gas density is small enough to make ionisation losses 
negligible. Finally, the advection velocity uy that appears in Eq. [4] has to be 
substituted with v4 for z < 0 and —v4 for z > 0. 

The steady-state (0 f/0t = 0) solution of the problem can be obtained by 
integrating the transport equation first between z = 07 and 07, and then 
between z = —o and 0-. Then, noticing that the solution must be symmetric 
with respect to z = 0 (f(z) = f(—z)) one gets: 


2 0 N ô 
2va (fo — fo) = 3P VA 2 os 


(p° Li fo) (48) 


where fo is the CR distribution function inside the cloud and fæ that at 
z = +o. The first and second terms in the right hand side of the equation 
account for the acceleration of particles in the converging flow of Alfvén waves, 
and for ionisation energy losses inside the cloud, respectively. Remarkably, the 
diffusion coefficient disappeared from the equation, and only the advective 
transport term (left hand side of the equation) is present. We will further 
discuss this below. 

Neglecting numerical factors of order unity, the ratio between the second 


and third term in Eq. |48]is of the order of (Morlino and Gabici}|2015): 


VA\ €H LE 
ae ee Nala (l 


where Y(E) = 1 + mge? /(E + mgc?) is a slowly varying function of energy 
whose value changes between 1 and 2. To simplify the discussion, we set Y = 1 
in the following. When R < 1 energy losses dominate, and the intensity of CRs 
is suppressed inside the cloud (fo < fo). 

Note that, as we saw above when discussing the free streaming transport 
scenario, the condition to be satisfied to have a suppression of the CR intensity 
in one single cloud crossing is more stringent, and would correspond to R < 
va/v « 1 (because CRs move much faster than Alfvén waves). This means that 
in the self regulating diffusive scenario two characteristic energies exist: Emin 
and Esr, which correspond to the conditions R = v4/v and R = 1, respectively. 


Low energy cosmic rays 67 


Therefore, FE», is larger than Emin, and corresponds to the energy below which 
the CR intensity is suppressed due to ionisation losses accumulated due to 
repeated cloud crossings. Particles of energy in the range Emin < E < Ep, are 
trapped in the vicinity of the cloud by the converging flow of Alfvén waves. 
Before losing its energy, a particle of energy FE», crosses the cloud v/v, times. 
The number of crossings decreases for lower energy particles, and for E < 
Esr particles lose their energy before a single crossing. Below Emin, then, 
the suppression of the CR intensity becomes more pronounced 
Gabici 2015). 

performed a detailed modelling of the penetration pro- 
cess, including the effect of streaming instability operating outside of the cloud, 
and assuming free streaming of CRs inside of it. The intensity of CRs far 
away from the cloud was assumed to be identical to the one measured in 
the local ISM. The streaming instability is triggered by CR protons, which 
are most abundant, and CR electrons and nuclei behave as passive particles 
in the amplified turbulence generated by protons. Their study was focussed 
on diffuse clouds, and they found total CR ionisation rates decreasing from 
< 3 x 10717 s-! to ~ 8 x 10718 s~! when the cloud column density increases 
from Ny, = 107° to 10?? cm~?. Such ionisation rates are a factor of 4-5 be- 
low those predicted within the framework of the free streaming transport. In 
a study in preparation (Phan et al. 2022), it will be shown that a similar 
difference is found also for dense clouds. 

The important conclusion of the studies reviewed above is that the standard 
picture of CRs pervading the Galactic disk with roughly the same intensity fails 
to reproduce the measurements of the CR ionisation rates in diffuse clouds. 
On the other hand, for dense clouds, predictions match observations only if 
the most optimistic (free streaming) scenario for the CR penetration in clouds 
is invoked. To date, this is one of the most puzzling problems in the study of 
LECRs in clouds, and possible solutions will be discussed in Sec. 

We conclude this Section by commenting on the absence of the diffusion 
coefficient in Eq. which is used to derive the spectrum of CRs inside the 
cloud fo. In fact, the particle diffusion coefficient D does play a role, as it 
defines a characteristic spatial scale of the problem, lą ~ D/va. When R < 1, 
lą represents the extension of the region around the cloud where the inten- 
sity of CRs is suppressed with respect to the interstellar value fæ. The value 
of the typical Galactic diffusion coefficient for sub-GeV protons is not well 
constrained. However, taking a quite plausible value of Drsm < 1078 cm?/s 


(‘Tatischeff et al. |2021), we get: 


Dism VA E JB, a 
la ~ 3 x 10? = 50 
~— (i Sg) (00 ona) Cae fe 


where the last term in parenthesis represents the ratio between the energy 
density of amplified magnetic fluctuations and its typical interstellar value (see 


Eq.[43). As pointed out by|Morlino and Gabicil (2015), streaming instability is a 


necessary ingredient of the model, as its effect is to increase 6B;,/dBy.15u and 
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reduce lą to values smaller than the typical coherence length of the interstellar 
magnetic field (lg ~ 10-100 pc). Without streaming instability, we would likely 
have la > lg, and the solution of the transport equation obtained in Eq. 
would no longer be valid. In that case, indeed, due to the mixing of magnetic 
field lines on scales larger than lg, the boundary condition f = fæ should be 
imposed at |z| ~ lg and not at |z| = œ. It can be shown that smaller values 
of lg would correspond to a more effective penetration of CRs into the cloud, 
and that the free streaming scenario can be recovered as a limiting case of the 
self regulating diffusive scenario for lg — 0. 

Let us assume, for a moment, that the amplification needed to have lg < lg 
does take place. In this case, the flux of CR into one side of the cloud can be 
written then ai} 


fe 4 vafo = vafa (51) 


of 
= Dz; +vafo~ D 

z 
where we assumed that D is uniform in space. This explain why D disappears 
from the solution of the transport equation, making the solution of the problem 
somehow universal. For R « 1 losses are so important to make fo « fæ, and 


therefore Eq reduces to (Skilling and Strong} |1976; |Cesarsky and Volk\ 
1978} |Morlino and Gabici} |2015): 


N ô 
2va fo = = 


2 
Pop (p* Li. fo) (52) 
which can be interpreted in a very simple way. The equilibrium spectrum of 
CRs inside the cloud can be derived by balancing the rate at which CRs of a 
given momentum enter the cloud (left hand side, the factor of 2 indicates that 
the cloud has 2 sides) with the rate at which they are removed from it due to 
energy losses (right hand side). 


6.4 Beam meets target: supernova remnant/molecular cloud associations 


The measurements of the CR ionisation rates reported in Fig. refer to 
passive clouds. The term passive is widely used in gamma-ray astronomy to 
indicate a cloud located in a region where no evidence of current or recent 
CR acceleration activity has been reported. The CR ionisation rate and the 
gamma-ray flux from a passive cloud are therefore ideal probes of the inter- 
stellar intensity of CRs which characterises the region of the Galaxy where the 
cloud is located (Aharonian| {2004). 

The presence of current or recent CR acceleration can be inferred, for 
example, from an enhanced gamma-ray flux from the cloud, either in the GeV 
domain, explored by the Fermi space mission, or in the TeV domain, probed 
by means of ground based observations with Cherenkov telescopes (e.g. 
(2009). Interestingly, the CR ionisation rate has been measured in few 


17 We neglect the small correction due to the Compton-Getting effect. 
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Fig. 13 CR ionisation rates measured in the vicinity of the supernova remnant IC 443 


(green points, |Indriolo et al.}2010), W51C (blue point,|Ceccarelli et al.!2011), and W28 (red 
points, 2014). Black points are as in Fig. 


clouds that also exhibit an enhanced gamma-ray emission. All these clouds are 
located in the vicinity of supernova remnants, which are very often invoked as 
the sources of CRs. 

The CR ionisation rates measured from these gamma-ray bright clouds are 
shown in Fig. as coloured data points, together with the ionisation rates 
measured in passive clouds, shown as black data points. Green, blue, and 
red data points refer to measurements taken in the vicinity of the supernova 


remnants IC 443 (Indriolo et al.| {2010), W51C (Ceccarelli et al.| 2011), and 
W28 (Vaupré et al.| |2014), respectively. |Indriolo et al.| (2010) derived the 
ionisation rates from H} absorption lines, while and 
from DCOt /HCO* abundance ratios. 


Fig.|13/shows that the values of Che found in clouds located in the vicinity 
of supernova remnants are systematically larger than the typical values found 
in passive clouds of the same column density. They are of the order of ee z 


10715 s71 and do not show any evident scaling with the cloud column density. 
An ionisation rate significantly smaller than 10~1° s7! was measured only in 
one case, from a line of sight towards a cloud located at a projected distance 
of ~ 10 pe from the shock of the supernova remnant W28 (red data point 
with dashed error bar). This probably means that the cloud is too far from 
the supernova remnant to be irradiated by the cosmic rays accelerated there. 
This hypothesis is corroborated by the fact that all the other red data points 
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in the figure refer to projected distances from the shock of W28 of at most a 
couple of parsecs. 

studied in more detail the case of the supernova rem- 
nant W28. They constrained the spectrum of GeV and TeV CR protons and 
electrons in the region, using available gamma-ray and radio observations. Re- 
markably, they found that an extrapolation of the CR proton spectrum of 
only an order of magnitude in particle energy (down to ~ 30 — 300 MeV) 
would suffice to explain the enhanced ionisation rates measured in the region. 
It seems, then, that CR protons accelerated at the supernova remnant shock 
can explain the enhancement of both the gamma-ray flux and the ionisation 
rates measured in the region. 


7 Open questions on low energy cosmic rays 


The right panel of Fig. [12] illustrates one of the most pressing issues in LECR 
studies: observed CR ionisation rates in interstellar clouds are very large. They 
are much larger than our expectations for diffuse clouds, while for dense clouds 
expectations match observations only if the most optimistic conditions for CR 
transport (free streaming) are invoked. Thus, the main question we should try 
to answer is: 


1. What induces the large ionisation rates observed in MCs? 


Answering this question is difficult because the measurements of ce pro- 
vide us with an integral constrain on the CR spectrum, which tells us nothing 
about the relative contribution from particles of different energies, or differ- 
ent species (CR nuclei or electrons) to the total ionisation rate. Improving 
the modelling of the transport of CRs in and around clouds, though certainly 
possible and desirable, would not help, as we saw that the problem persists 
also in the idealised scenario where CRs penetrate into clouds moving along 
straight lines. 

proposed that the large ionisation rates could be 
due to an unseen component of CRs of very low energies, popping up in the 
sub-MeV domain, currently unconstrained by observations. However, we can 
see from the left panel of Fig. [12] that energy losses are particularly severe at 
these very low energies. For example, in an ISM of density nrsm ~ 1 cm7? 
the ionisation loss time for a proton (electron) of energy smaller than 1 MeV 
would be shorter than 7*,, ~ 2 x 104 yr (10° yr), see e.g. Fig. 2 in 
(2018). Moreover, inside clouds of density nea > nrg, the energy loss 
time would be a factor of na/nrsm shorter than that. It is clear, then, that 
maintaining such a population of very low energy CRs in the ISM and/or in 
clouds would require an excessively large amount of energy. This is because, to 
balance energy losses, freshly accelerated particles would have to be injected in 
the ISM at a very fast rate o1/r*,, (see [Recchia et al.| for quantitative 
estimates). 

A thing that should be kept in mind is that the discrepancy between ob- 
servations and expectations was found under the assumption that the spectra 
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of CRs measured in the local ISM are representative of the entire Galactic 
disk. Such an assumption is very well justified for CR protons in the GeV 
energy domain, as discussed in Sections [4.1.1] and and illustrated in the 
top panel of Fig. p but may not hold for lower energy CRs. Therefore, a very 
fundamental question emerges: 


2. Are the spectra of LECRs measured in the local ISM representative of the 
entire Galaxy? Or, what is the spatial distribution of LECRs throughout 
the Galactic disk, on both large and small scales? 


In an attempt to answer this questions, |Indriolo et al.| (2015) plotted the 


CR ionisation rate as a function of the galactocentric distance R (see their 
Fig. 23). The values of aan where obtained from the abundance of oxygen 
bearing ions in diffuse atomic clouds (see Sec. [5.2.1). They found no relation 
between the ionisation rate and the galactocentric distance for R > 5 kpc, and 
a moderate increase of Ce towards the Galactic centre for smaller values of 
R. Such a trend was predicted by [Wolfire et al.| (2003), as a result of the larger 
concentration of potential CR sources towards the inner Galaxy, coupled with 
the spatial distribution of the gas density which determines the effectiveness 


of energy losses. The prediction by |Wolfire et al.| (2003), rescaled to match 


available observations, is shown as a dashed line in the bottom panel of Fig. pl 
(2017) refined the analysis by considering only the lines of 
sights from for which also a detection of ArH* was avail- 
able. They then used a detailed diffuse cloud model to fit the abundances of 
both oxygen bearing ions and argonium. Their estimates of the CR ionisation 
rates are shown as blue data points in the bottom panel of Fig. pj and show 
no correlation with galactocentric distancdt’| Note that this is not necessarily 
in contradiction with the claim made by (2015), as only two 
clouds in the refined sample are located at R < 5 kpc. The yellow shaded 
region in the Figure indicates the expected level of the CR ionisation rate in 
diffuse clouds obtained by {Phan et al.| under the assumption that the 
intensity of CRs at any location in the Galactic disk is identical to the local 
one, and assuming a self regulating diffusive penetration of LECRs into clouds 
(see Sec. (6.3). The red bar marks the Sun’s location in the Galaxy. Data are 
a factor of ~ 10 — 100 above such predictions (the discrepancy would reduce 
to a factor of ~ 10 in the free streaming transport scenario, see e.g. 
or right panel of Fig. (12h. 

At this point, it is worth reminding that for a large number of diffuse clouds 
only an upper limit for the CR ionisation rate could be derived (see Fig. o), 
implying that we do not know how broad is the distribution of values of Ç} 
in such objects. That is to say: the blue points in the bottom panel of Fig. [9] 
might represent the upper envelope of the real distribution of CR ionisation 
rates, which may vary significantly from cloud to cloud. This suggests that the 
discrepancy between observations and expectations might be reduced if small 


18 The values of the primary CR ionisation rate per H atom reported by [Neufeld and 
T: : H2 
Wolfire| (2017) were multiplied by 2.3 to obtain CGR (Glassgold and Langer| |1974}. 
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scale fluctuations exist in the spatial distribution of LECRs, and the Solar sys- 
tem is located in a minimum of such fluctuating distribution. Such fluctuations 
are indeed expected, as LECRs do not travel far from their acceleration sites 
due to the severe ionisation losses suffered in the ISM. This idea, pioneered 
by [Cesarsky| (1975), was recently investigated in detail by {Phan et al.| (2021), 
who found that fluctuations might be relevant if LECRs are accelerated at 
supernova remnants (the rate of about 3 supernova explosions per century 
in the Galaxy determines both the average distance between sources and the 
CR power per source, which are the two key parameters of the problem). A 
comparison with measured ionisation rates is in progress (Phan et al. 2022, in 
preparation). 

The search for a dependence of a on galactocentric radius has been 
based on data on atomic diffuse clouds. This is because these objects are char- 
acterised by quite small gas column densities, and therefore the suppression 
of the CR intensity inside the clouds is likely to be quite moderate. As a con- 
sequence, the measured ionisation rate is expected to be quite close to the 
pristine interstellar one. However, Fig. 25 from |Indriolo et al.| shows 
that there is no appreciable difference in the ionisation rate of atomic and 
molecular diffuse clouds, for column densities varying over almost two decades 


(see also|Neufeld and Wolfire} (2017). Answering the question 


3. Why are diffuse atomic and MCs, despite their different column density, 
characterised by the same ionisation rate? 


will help us to better understand the transport of LECRs in and around clouds. 

Measurements of the CR ionisation rate in the Galactic centre region are 
shown in the bottom panel of Fig. [9| as a green data point and error bar 
au The ionisation rate is of the order 
of = 10-** s~“*, which is more than an order of magnitude larger than the 
ionisation rates measured in diffuse clouds (blue data points), and more than 
three orders of magnitude larger than the predictions from the self regulating 
transport scenario (yellow shaded region), leading to the question: 


4. Why is the ionisation rate so large in the Galactic centre region? 


The Galactic centre hosts a supermassive black hole, surrounded by the 
most massive MC complex in the Galaxy, the Central Molecular Zone, of 
radius ~ 100 pc. The region is characterised by an intense star formation 
rate, and a correspondingly enhanced supernova rate. It also shows signs of 
nuclear activity in the form of powerful outflows (e. g./Heywood et al.}[2022} and 
references therein). Given this rich phenomenology, one might expect to find 
in that region a largely enhanced intensity of CRs at all energies. Somewhat 
surprisingly, gamma-ray observations revealed that the excess in high energy 
(beyond GeV) CRs is present but moderate, pointing towards an effective 
rate of escape of CRs from the inner Galaxy (2017). Explaining 
the extremely large value of the ionisation rate would then require either a 
very fine tuned mechanism able to confine sub-GeV CRs (most effective in 
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ionising the gas) but not particles of higher energy (constrained by gamma- 
ray observations), or the presence of other ionising agents (UV/X-ray photons 
pervading the central molecular zone?). 

A possible way to explain the large discrepancy between observations and 
expectations is to assume that the LECRs responsible for the ionisation of 
clouds are produced in situ, inside the clouds. In this scenario, they would 
form a distinct population of particles, with no links with local interstellar 
LECRs. Ongoing particle acceleration in protostellar environments has been 
revealed by the presence of non-thermal synchrotron emission (and therefore 


of relativistic electrons) and of enhanced ionisation rates (e.g. 


Kamenetzky et al.||2017 2018). Shocks in protostellar jets seems to 


be the best places where to accelerate LECRs in these environments 
2016). Under certain conditions, LECRs accelerated in protostellar en- 
vironments might contribute to the ionisation rate in the cloud and reduce the 
discrepancy with observations (2019). To test this hypothesis, 
one should confirm observationally that protostellar particle accelerators are 
present in all clouds characterised by a large value of Ga 

In another scenario, ionisation losses suffered by LECRs could be balanced 
by a continuous acceleration of particles due to their interaction with turbulent 
reconnecting regions. As turbulence is ubiquitous in clouds, such mechanism 
is expected to produce homogeneous distributions of cosmic rays in the entire 
cloud volume [2021b). The aim of these lines of research is to 


answer the question 


5. What are the sites of acceleration of the LECRs responsible for the ionisa- 
tion of clouds? 


Reversing the argument, instead of invoking an enhanced intensity of CRs 
inside clouds, we might speculate on the existence of a mechanism able to 
suppress the local intensity of LECRs. In this case, the local CR. measurements 
would have to be demodulated to obtain the true interstellar spectra. This 
scenario is very appealing, as the Solar system is located within a cavity in the 
ISM, called the local bubble. Such a cavity has been inflated by the repeated 
supernova explosions that took place in a star cluster, and is surrounded by a 
shell of dense and cold gas 2022). It 
has been suggested that ionisation energy losses experienced by LECRs while 
crossing the shell could indeed suppress the spectrum of CRs in the MeV 


domain (Silsbee and Ivlev} [2019] 2020). This raises the question: 


6. Does the observed intensity of LECRs in the local ISM reflects the fact 
that we live in a special place in the Galaxy? 


Finally, we should also consider the possibility that the uncertainties in our 
knowledge of chemical reaction networks might impact on the observational 
estimates of the ionisation rates. As discussed in Sec. [5.2] 2 [5.2] ignoring the presence 


of PAH when defining the chemical network used to poe He observations 
might lead to significant errors in the derivation of eee Ser Ne 
2021 


[Liszt] [2003} [Hollenbach et al. 2003; [Hollenbach et al. et al. [2012| [Shaw and Ferland Shaw and Ferland We also saw 
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that, in dense clouds, estimates of en are affected by the poor knowledge of 
the fraction of carbon and oxygen depleted on dust grains (factors fg and 6 in 
Eq. [37). Moreover, also the assumption of a uniform CR ionisation rate inside 
clouds introduces errors, and astrochemical models should include the a 
of a position (column density) dependent CR ionisation rate 


[2012} [Gaches et al.| [2021a). The last question that we should ask is then: 


7. How reliable are the observational estimates of the CR ionisation rates in 
interstellar clouds? 


8 Future perspectives 


During the past two decades, a number of observational breakthroughs in the 
study of LECRs changed radically our view on these particles. Remarkably, 
progresses in this field of research spun across many branches of astrophysics, 
including space exploration (the Voyagers crossing the boundary of the helio- 
sphere in 2012), laboratory astrophysics (the measurement of the dissociative 
recombination rate of H} in 2003), gamma-ray astronomy (the unprecedented 
view of the Galaxy provided by Fermi/LAT, in orbit since 2008), CR detection 
from space (PAMELA operated until 2016, and AMS-02 since 2011), infrared 
and millimeter astronomy (the launch of Herschel in 2009, and the numerous 
ground based observational campaigns). Theoretical studies, as often happens, 
are trying to catch up with these observational advancements, but the open 
questions listed in the previous Section demonstrate that we are still far from 
a satisfactory understanding of LECRs. 

Luckily, the next two decades promise to be equally exciting. Moving from 
low to high energies, here is a list of some of the most important advancements 
we might expect to witness in the next 20 years. 

The Square Kilometer Array Observatory (SKAO) -Il will be the largest 
array of radio telescopes ever built. It is expected to be fully operational in 
the late twenties. As discussed in Sec. radio observations of clouds with 
instruments of superior sensitivity will help us in better characterise the pop- 
ulation of LECR electrons in such objects. 

The James Webb Space Telescope (JWST) - Launched on 2021’s Christ- 
mas day, it will dramatically improve our view on the infrared sky, which 
is one of the main sources of information we have on LECRs (see Sec. 5-1). 
Moreover, it has been recently proposed that JWST could detect the near- 
infrared photons produced in the decay of rotovibrational levels of molecular 
hydrogen. In interstellar clouds, the excitation of these levels is dominated by 
secondary CR electrons impacting on H2. This implies: that the observations of 
such transitions will provide us with an aaa of eee cpr Which is independent 


on any chemical network (Padovani et al. . Ina similar fashion, emission 


lines from excited H} might also be eee in the infrared band from clouds 


characterised by ites ionisation rates (Becker et al. Becker et al.|{2011). 


The Athena X-ray observatory — Expected to be launched in the mid thir- 
ties, it will explore with unprecedented sensitivity the X-ray sky. Of particular 
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relevance for LECR studies is the Fe Ka line at 6.4 keV, which is produced 
when a cold gas is irradiated with either X-rays or energetic particles. It re- 
sults from the removal of a K-shell electron followed by a transition from the 
L shell to fill the vacancy (e.g. (2012). The 6.4 keV line has 
been detected from clouds near the Galactic centre. However, its flux varies 
with time in a way which is not compatible with an origin due to LECRs, but 
is rather connected to photon irradiation from flares originating from the cen- 
tral supermassive black hole (e.g. and references therein). 
Nevertheless, upper limits on the steady component of such emission can be 
converted into upper limits on the intensity of LECRs in the Galactic centre 
region. The most stringent constraints on the intensity of LECRs obtained in 
this way are of the same order of the value of the CR proton intensity needed 
to explain the enhanced ionisation rate measured in the Galactic centre region 
(2021). Future observations are then of paramount importance, 
as a decrease in the line flux would rule out a CR proton origin of the ionisation 
rates. 

Tentative evidence for the presence of 6.4 keV line emission has been 
presented also for a handful of clouds interacting with supernova remnants 
(Nobukawa et al.||2018). The emission, if confirmed, could be due to the pres- 
ence of LECRs accelerated at the supernova remnant shocks. Also in this case, 
future observations will help in clarifying the role played by LECRs close to 
putative accelerators, and will possibly lead to further detections. 

A renaissance of MeV astronomy? — Gamma-ray lines in the 0.1-10 MeV 
energy domain are produced following the collisions of LECRs with interstel- 
lar matter, as a result of the de-excitation of the first nuclear level. They are 
expected to be particularly intense for interstellar 12C, 160, Ne, ?°Mg, ?°Si, 
and °°Fe [1979). Unfortunately, the MeV domain is poorly 
explored, but the situation might change in a near future. NASA’s Compton 
Spectrometer and Imager (COSI), expected to be launched in 2025, will re- 
open the MeV window, and hopefully pave the road to future mid-scale MeV 
missions. Remarkably, if the large intensity of LECRs needed to explain the 
CR ionisation rates of diffuse clouds is a proxy for the typical Galactic value, 
de-excitation nuclear line emission from the inner Galaxy would be within 
reach of any of such mid-scale future missions (2018). 

Theoretical/phenomenological studies will have to focus on both the very 
local and the remote interstellar medium. Future observations of clouds from 
the radio to the gamma-ray domain will help in constraining models of in situ 
acceleration of LECRs, while local observations on the composition of LECRs 
have been recently used to show that such particles are likely accelerate in 
a diffuse and hot medium, characteristic of superbubbles 
(2021). The latter result brings us back to the fundamental question: does our 
location inside the local bubble affects our view on LECRs? 


Acknowledgements I would like to thank the editors at A&A Rev, and in particular 
Luigina Feretti for her advices and infinite patience. I also thank Phan Vo Hong Minh, 
Sarah Recchia, and Thibault Vieu, who read and commented on the entire manuscript, 


76 Stefano Gabici 


Paola Caselli, with whom I had (almost 20 years ago) my first discussion on LECRs, and 
Cecilia Ceccarelli and Thierry Montmerle who gave me the opportunity to start working 
in this field. While writing this review I enjoyed discussing LECRs and related topics with 
G. Bernardi, D. Breitschwerdt, J. Duprat, C. Evoli, D. Gaggero, D. Galli, R. Giuffrida, D. 
Grasso, A. Ivlev, D. Maurin, P. Mertsch, M. Miceli, M. Padovani, G. Peron, L. Podio, S. 
Ravikularaman, J. Raymond, G. Sabatini, F. Schulze, A. Strong, V. Tatischeff, R. Terrier, 
and F. van der Tak. Finally, I acknowledge support from Agence Nationale de la Recherche 
(grants ANR- 17-CE31-0014 and CRitiLISM). 


References 


Acero F, Ackermann M, Ajello M, et al. (2016) Development of the Model of 
Galactic Interstellar Emission for Standard Point-source Analysis of Fermi 
Large Area Telescope Data. Astrophys J Suppl 223(2):26. https: //doi.org/ 


10.3847 /0067-0049 /223/2/26) [ar Xiv:1602.07246| [astro-ph.HE] 


Ackermann M, Ajello M, Atwood WB, et al. Raa Fermi- ne Observations 
of the Diffuse y-Ray Emission: Implications for Cosmic Rays and the Inter- 


stellar Medium. Astrophys J 750(1):3. https: //doi.org/10.1088/0004-637X / 
750/1/3) {arXiv:1202.4039] [astro-ph.HE] 


Adriani O, Barbarino GC, Bazilevskaya GA, et al. (2013) Time Dependence 
of the Proton Flux Measured by PAMELA during the 2006 July-2009 De- 
cember Solar Minimum. Astrophys J 765(2):91. https://doi.org/10.1088/ 
fastro-ph. HE) 


Adriani O, Barbarino GC, Bazilevskaya GA, et al. (2014) The PAMELA 
Mission: Heralding a new era in precision cosmic ray physics. Phys Rep 
544(4):323-370. https: / /doi.org/10.1016/j.physrep.2014.06.003 

Aguilar M, Aisa D, Alpat B, et al. (2015) Precision Measurement of the Proton 
Flux in Primary Cosmic Rays from Rigidity 1 GV to 1.8 TV with the Alpha 
Magnetic Spectrometer on the International Space Station. Phys Rev Lett 
114(17):171103. {https: //doi.org/10.1103/PhysRevLett.114.171103 

Aguilar M, Ali Cavasonza L, Ambrosi G, et al. (2016) Precision Measurement 
of the Boron to Carbon Flux Ratio in Cosmic Rays from 1.9 GV to 2.6 TV 
with the Alpha Magnetic Spectrometer on the International Space Station. 
Phys Rev Lett 117(23):231102. |https://doi.org/10.1103/PhysRevLett.117. 

Aguilar M, Ali Cavasonza L, Alpat B, et al. (2017) Observation of the Identical 
Rigidity Dependence of He, C, and O Cosmic Rays at High Rigidities by 
the Alpha Magnetic Spectrometer on the International Space Station. Phys 
Rev Lett 119(25):251101. https: //doi.org/10.1103 /PhysRevLett.119.251101 

Aguilar M, Ali Cavasonza L, Alpat B, et al. (2019) Towards Understanding 
the Origin of Cosmic-Ray Electrons. Phys Rev Lett 122(10):101101. 
7/doi.org/10.1103/PhysRevLett.122.101101 

Aguilar M, Ali Cavasonza L, Ambrosi G, et al. (2021) The Alpha Magnetic 
Spectrometer (AMS) on the international space station: Part II - Results 


from the first seven years. Phys Rep 894:1-116. https://doi.org/10.1016/j. 
physrep.2020.09.003 


Low energy cosmic rays 77 


Aharonian F, Yang R, de Oña Wilhelmi E (2019) Massive stars as major 
factories of Galactic cosmic I Nature AT, 3:561-567. https://doi. ftps //doi] 

An Dae Ge vans Ie F, Peron G, PeT: T aa T 5, ah R (2020) a the 
sea of galactic cosmic rays with Fermi-LAT. Phys Rev D 101(8):083018. 
https: / /doi.org/10.1103/PhysRevD.101.083018 

Aharonian FA (1991) Vary High and Ultra High Energy Gamma-Rays from 
Giant Molecular Clouds. Astrophys Space Sci 180(2):305-320. 


Aharonian FA (2004) Very high energy cosmic gamma radiation : a crucial 
window on the extreme Universe. https://doi.org/10.1142/4657 

Amato E, Casanova S (2021) On particle acceleration and transport in 
plasmas in the Galaxy: theory and observations. Journal of Plasma 


Physics 87(1):845870101. _https://doi.org/10.1017/50022377821000064 
arXiv:2104.12428) [astro-ph.HE] 


AMS Collaboration, Aguilar M, Alcaraz J, et al. (2002) The Alpha Magnetic 
Spectrometer (AMS) on the International Space Station: Part I - results 
from the test flight on the space shuttle. Phys Rep 366(6):331—405. 
//doi.org/10.1016/S0370-1573(02)00013-3 

Armitage PJ (2011) Dynamics of Protoplanetary Disks. Ann 
Rev Astron Ra TTR aT Seer TOS i 195- a https: on org/10.1146/ 


A Ee WI, Leer E, akon 5 (1977) The Se ee of aN Rays by 
Shock Waves. In: a nls Ray Conference. International Cos- 
mic Ray Conference, vol 11. p 132 

Baade W, Zwicky F (1934) Cosmic Rays from Super-novae. Proceedings of the 
National Academy of Science 20(5):259-263. https: //doi.org/10.1073/pnas. 

Bacalla XL, Linnartz H, Cox NLJ, et al. (2019) The EDIBLES survey. IV. 
Cosmic ray ionization rates in diffuse clouds from near-ultraviolet observa- 


tions of interstellar OH+. Astron Astrophys 622:A31. |https://doi.org/10. 


1051/0004-6361/201833039, arXiv:1811.08662 [astro-ph.GA] 


Barlow MJ, Swinyard BM, Owen PJ, et al. TAE Detection of a Noble Gas 
Molecular Ion, 36 ArH*, in the Crab Nebula. Science 342(6164):1343-1345. 
https://doi.org/10.1126/science.1243582 org/10. 1126/science.1243582 jarXiv:1312.4843 | [ arXiv:1312.4843 [astro-ph.GA] 
Becker JK, Black JH, Safarzadeh M, Schuppan F (2011) Tracing the Sources 
of Cosmic Rays with Molecular Ions. Astrophys J Lett 739(2):L43. 


//doi.org/10.1088/2041-8205/739/2/L43) arXiv:1106.4740) [astro-ph.GA] 


Bell AR (1978) The acceleration of cosmic rays in paper re - I. Mon Not 


Royal Astron Soc 182:147—-156. https: //doi-org/10.1093/mnras/182.2.147 


Berezinskii VS, Bulanov SV, Dogiel VA, Ptuskin VS (1990) Astrophysics of 
cosmic rays. Amsterdam: North-Holland 

Bergin EA, Tafalla M (2007) Cold Dark Clouds: The Initial See for Star 
Formation. Ann Rev Astron Astrophys 45(1):339-396. 


1146/annurev.astro.45.071206.100404, |jarXiv:0705.3765 aips m 


Beuermann K, Kanbach G, Berkhuijsen EM E Role structure of the 


78 Stefano Gabici 


Galaxy: thick disk and thin disk at 408 MHz. Astron Astrophys 153:17-34 
Bindi V, Corti C, Consolandi C, Hoffman J, Whitman K (2017) Overview of 
galactic cosmic ray solar modulation in the AMS-02 era. Advances in Space 


Research 60(4):865-878. https: //doi.org/10.1016/j.asr.2017.05.025 


Black JH, Dalgarno A (1973) The Cosmic Abundance of Deuterium. Astrophys 


J Lett 184:L101. https://doi.org/10.1086/181299 


Black JH, Dalgarno A (1977) Models of interstellar clouds. I. The Zeta Ophi- 
uchi cloud. Astrophys J Suppl 34:405—423. https: //doi.org/10.1086/190455 
Black JH, Fazio GG (1973) Production of Gamma Radiation in Dense In- 
terstellar Clouds by Cosmic-Ray Interactions. Astrophys J Lett 185:L7. 


https: //doi.org/10.1086/181310 


Blandford R, Eichler D (1987) Particle acceleration at astrophysical shocks: 
A theory of cosmic ray origin. Phys Rep 154(1):1-75. 
1016/0370-1573(87)90134-7 

Blandford RD, Ostriker JP (1978) Particle acceleration by astrophysical 


shocks. Astrophys J Lett 221:L29-L32. https://doi.org/10.1086/182658 


Blasi P (2013) The origin of galactic cosmic rays. Astron Astrophys Rev 21:70. 
https://doi.org/10.1007/s00159-013-0070-7, |arXiv:1311.7346| [astro-ph.HE] 

Blumenthal GR, Gould RJ (1970) Bremsstrahlung, Synchrotron ee 
tion, and Compton Scattering of High-Energy Electrons Traversing Dilute 
Gases. Reviews of Modern Physics 42(2):237-271. |https: //doi.org/10.1103/ 

Boezio M, Munini R, Picozza P (2020) Cosmic ray detection in space. Progress 
in Particle and Nuclear Physics 112:103765. https: //doi.org/10.1016/j.ppnp. 
2020.103765 

Bohlin RC, Savage BD, Drake JF (1978) A survey of interstellar H I from 
Lalpha absorption measurements. II. Astrophys J 224:132-142. |https://doi.] 

Bolatto AD, Wolfire M, Leroy AK (2013) The CO-to-Hə Conversion Fac- 
tor. Ann Rev Astron Astrophys 51(1):207-268. https://doi.org/10.1146/ 
astro-ph Gal 


Boschini MJ, Della Torre S, Gervasi M, La Vacca G, Rancoita PG (2019) The 
HELMOD model in the works for inner and outer heliosphere: From AMS 
to Voyager probes observations. Advances in Space Research 64(12):2459- 
2476. https: //doi.org/10.1016/j.asr.2019.04.007 

Bovino S, Ferrada-Chamorro S, Lupi A, Schleicher DRG, Caselli P (2020) A 
new proxy to estimate the cosmic ray ionization rate in dense cores. Mon Not 


Royal Astron Soc 495(1):L7—-L11. |https://doi.org/10.1093/mmnrasl/slaa048 
arXiv:2003.05416] [astro-ph.GA] 


Breitschwerdt D, Feige J, Schulreich MM, et al. (2016) The locations of 
recent supernovae near the Sun from modelling °°Fe transport. Nature 
532(7597):73-76. https: //doi.org/10.1038/naturel 7424 

Brittain SD, Simon T, Kulesa C, Rettig TW (2004) Interstellar H3+ Line 
Absorption toward LkHa 101. Astrophys J 606(2):911-916. 


10.1086 /383024 


Brown RL, Marscher AP (1977) Energetic secondary electrons in dense inter- 


Low energy cosmic rays 79 


stellar clouds. Astrophys J 212:659-663. https://doi.org/10.1086/155088 


Burbidge EM, Burbidge GR, Fowler WA, Hoyle F (1957) Synthesis of the 
Elements in Stars. Reviews of Modern Physics 29(4):547-650. 

Bykov AM (2014) Nonthermal particles and photons in starburst regions 
and superbubbles. Astron Astrophys Rev 22:77. https://doi.org/10.1007/ 
0159-014 0077-8) ar Xiv:151 1.04608 [astro-ph HE 


Caballero-Lopez RA, Moraal H (2004) Limitations of the force field equation 
to describe cosmic ray modulation. Journal of Geophysical Research (Space 


Physics) 109(A1):A01101. /https://doi.org/10.1029/2003A010098 


Cameron AGW (1957) Nuclear Reactions in Stars and Nucleogenesis. Publ 


Astron Soc Pac 69(408):201. https: //doi-org/10.1086/127051 


Cardelli JA, Meyer DM, Jura M, Savage BD (1996) The Abundance of Inter- 
stellar Carbon. Astrophys J 467:334. https://doi.org/10.1086/177608 

Casandjian JM (2015) Local H i Emissivity Measured with i and 
Implications for Cosmic-Ray Spectra. Astrophys J 806(2):240. 


org/10.1088/0004-637X /806/2/240) |arXiv:1506.00047| (astro-ph. TE 


Casanova S, Aharonian FA, Fukui Y, et al. ( a a Clouds as Wak 
Ray Barometers. Publ Astron Soc Japan 62:769. https: //doi.org/10.1093/ 
pasj /62.3.769, arXiv:0904.2887 [astro-ph.HE] 


Caselli P (2002) painoo nants Se as a eo of ionization fraction in dense 
interstellar clouds. Planet Space Sci 50(12-13):1133-1144. |https://doi.org/ 


10.1016 /S0032-0633(02)00074-0! arXiv:astro-ph/0204127) [astro-ph] 


Caselli P, Ceccarelli a abt} Our astrochemical heritage. ees Astrophys 
Rev 20:56. [https://doi.org/10.1007/s00159-012-0056-x, /arXiv:1210.6368) 
[astro-ph.GA] 

Caselli P, Walmsley CM, Terzieva R, Herbst E (1998) The Ionization Frac- 
tion in Dense Cloud Cores. Astrophys J 499(1):234-249. 

Casse F, Lemoine M, Pelletier G (2001) Transport of cosmic rays in 
chaotic magnetic fields. Phys Rev D 65(2):023002. https://doi.org/10.1103/ 
[astro-ph 


Ceccarelli C, Hily-Blant P, Montmerle T, et al. (2011) Supernova-enhanced 
Cosmic-Ray Ionization and Induced Chemistry in a Molecular Cloud of 


W51C. Astrophys J Lett 740(1):L4. [https://doi.org/10.1088/2041-8205/ 
740/1/L4| larXiv:1108.3600| [astro-ph.GA] 


Ceccarelli C, Caselli P, Bockelée-Morvan D, et al. (2014) Deuterium Frac- 
tionation: The Ariadne’s Thread from the Precollapse Phase to Mete- 
orites and Comets Today. In: Beuther H, Klessen RS, Dullemond CP, Hen- 


ning T (eds) Protostars and Planets VI. p 859. https://doi.org/10.2458/ 
azu_uapress_9780816531240-ch037, arXiv:1403.7143) [astro-ph.EP] 


Cesarsky CJ (1975) Interstellar Propagation of Low Energy Cosmic Rays. In: 
International Cosmic Ray Conference. International Cosmic Ray Confer- 
ence, vol 2. p 634 

Cesarsky CJ, Montmerle T (1983) Gamma-Rays from Active Regions in 
the Galaxy - the Possible Contribution of Stellar Winds. Space Sci Rev 


80 Stefano Gabici 


36(2):173-193. https: //doi.org/10.1007/BF 00167503 


Cesarsky CJ, Volk HJ (1978) Cosmic Ray Penetration into Molecular Clouds. 
Astron Astrophys 70:367 
Chabot M (2016) Cosmic-ray slowing down in molecular clouds: Effects 


of heavy nuclei. Astron Astrophys 585:A15. https://doi.org/10.1051/ 
0004-6361/201425441 


Chandran BDG (2000) Confinement and Isotropization of Galactic Cosmic 
Rays by Molecular-Cloud Magnetic Mirrors When Turbulent Scattering Is 
Weak. Astrophys J 529(1):513-535. https: //doi.org/10.1086/308232 

Cholis I, Hooper D, Linden T (2016) A predictive analytic ae for the solar 
modulation of cosmic rays. Phys Rev D 93(4):043016. 
1103/PhysTevD 93.043016 ar X151 L01507 fastro-ph-SR 


Chuvilgin LG, Ptuskin VS (1993) Anomalous diffusion of cosmic rays across 
the magnetic field. Astron Astrophys 279(1):278-297 

Clark GW, Garmire GP, Kraushaar WL (1968) Observation of High-Energy 
Cosmic Gamma Rays. Astrophys J Lett 153:L203. https: //doi.org/10.1086/ 


Coc A, Vangioni E (2017) Primordial nucleosynthesis. International 
Journal of Modern Physics E 26(8):1741002. https://doi.org/10.1142/ 
50218301317410026, arXiv:1707.01004| [astro-ph.CO] 


Compton AH, Getting IA (1935) An Apparent Effect of Galactic Rotation 
on the Intensity of Cosmic Rays. Physical Review 47(11):817-821. 
//doi.org/10.1103/PhysRev.47.817 

Corti C, Bindi V, Consolandi C, Whitman K (2016) Solar Modulation of the 
Local Interstellar Spectrum with Voyager 1, AMS-02, PAMELA, and BESS. 
Astrophys J 829(1):8. https://doi.org/10.3847/0004-637X/829/1/8 

Corti C, Potgieter MS, Bindi V, et al. (2019) Numerical Modeling of Galac- 
tic Cosmic-Ray Proton and Helium Observed by AMS-02 during the Solar 
Maximum of Solar Cycle 24. Astrophys J 871(2):253. 

Cox DP (2005) The Three-Phase Interstellar Medium Revisited. Ann Rev 
Astron Astrophys 43(1):337-385. https://doi.org/10.1146/annurev.astro.43. 
072103.150615 

Cravens TE, Dalgarno A (1978) Ionization, dissociation, and heating efficien- 
cies of cosmic rays in a gas of molecular hydrogen. Astrophys J 219:750-752. 
https: //doi.org/10.1086/155833 

Crocker RM, Jones D, Protheroe RJ, et al. (2007) The Cosmic Ray Distribu- 
tion in Sagittarius B. Astrophys J 666(2):934-948. https://doi.org/10.1086/ 
(518243) [astro-ph] 

Crutcher RM (2012) Magnetic Fields in Molecular Clouds. Ann Rev Astron 
Astrophys 50:29-63. https://doi.org/10.1146/annurev-astro-081811- 125514 

Cummings AC, Stone EC, Heikkila BC, et al. (2016) Galactic Cosmic Rays in 
the Local Interstellar Medium: Voyager 1 Observations and Model Results. 
Astrophys J 831(1):18. https://doi.org/10.3847/0004-637X /831/1/18 

Dalgarno A (2006) Interstellar Chemistry Special Feature: The galactic cos- 
mic ray ionization rate. Proceedings of the National Academy of Science 


Low energy cosmic rays 81 


103(33):12269-12273. {https: //doi.org/10.1073/pnas.0602117103 


Dartnell LR (2011) Ionizing Radiation and Life. Astrobiology 11(6):551—-582. 
https: //doi.org/10.1089/ast.2010.0528 

de Angelis A, Tatischeff V, Grenier IA, et al. (2018) Science with e- 
ASTROGAM. A space mission for MeV-GeV gamma-ray astrophysics. Jour- 


nal of High Energy Astrophysics 19:1—-106. https://doi.org/10.1016/j.jheap. 
2018.07.001, arXiv:1711.01265 [astro-ph.HE] 


Di Bernardo G, Evoli C, Gaggero D, Grasso D, Maccione L (2013) Cosmic ray 
electrons, positrons and the synchrotron emission of the Galaxy: consistent 
analysis and implications. J Cosm Astropart Phys 2013(3):036. 
org/10.1088/1475-7516/2013/03/036 [astro-ph.HE] 

Di Bernardo G, Grasso D, Evoli C, Gaggero D (2015) Diffuse synchrotron 
emission from galactic cosmic ray electrons. ASTRA Proceedings 2:21-26. 

https: //doi.org/10.5194/ap-2-21-2015 Be eel [astro-ph.HE] 

Dickinson C, Beck R, Crocker R, et al. (2015) SKA studies of in situ syn- 
chrotron radiation from molecular clouds. In: Advancing Astrophysics with 
the Square Kilometre Array (AASKA14). p 102. [astro- 
ph.GA] 

Dogiel VA, Chernyshov DO, Ivlev AV, Kiselev AM, Kopyev AV (2021) Self- 
modulation of Cosmic Rays in Molecular Clouds: Imprints in the Radio 


Observations. Astrophys J 921(1):43. https://doi.org/10.3847 /1538-4357 / 
ar Xiv:2108.07559] [astro-ph.GA] 


Dolginov AZ, Toptygin IN (1968) Cosmic rays in the interplanetary magnetic 
fields. Icarus 8(1-3):54—60. |https: //doi.org/10.1016/0019-1035(68)90061-4 
Donato F, Maurin D, Taillet R (2002) beta -radioactive cosmic rays in 
a diffusion model: Test for a local bubble? Astron Astrophys 381:539- 
559. |https://doi.org/10.1051/0004-6361:20011447 

[astro-ph] 

Drury LOC (2012) Origin of cosmic rays. Astroparticle Physics 39:52- 
60. [https://doi.org/10.1016/j.astropartphys.2012.02.006 
[astro-ph.HE] 

Everett JE, Zweibel EG (2011) The Interaction of Cosmic Rays with Diffuse 


Clouds. Astrophys J 739(2):60. https: //doi.org/10.1088/0004-637X/739/2/ 
astro-ph. GA 

Evoli C, Aloisio R, Blasi P (2019) Galactic cosmic rays after the AMS-02 ob- 
Hee pe Rev D Ee OA i 103023. Tepe ees Peer: 


nas U ( STS) Laine Call ee E Vi of a Fast Particles and the pes 
Strength of Optical Transitions. Physical Review 95(5):1198—1200. 
//doi.org/10.1103/PhysRev.95.1198 

Farmer AJ, Goldreich P (2004) Wave Damping by Magnetohydrodynamic 
Turbulence and Its Effect on Cosmic-Ray Propagation in the Interstel- 


lar Medium. Astrophys J 604(2):671-674. https: //doi.org/10.1086/382040 
[astro-ph] 

Favre C, Ceccarelli C, López-Sepulcre A, et al. (2018) SOLIS IV. Hydrocar- 
bons in the OMC-2 FIR4 Region, a Probe of Energetic Particle Irradiation 


82 Stefano Gabici 


of the Region. Astrophys J 859(2):136. |https: //doi.org/10.3847/1538-4357 / 
aabfd4, jar Xiv:1804.07825) [astro-ph.GA] 


Federman SR, Weber J, Lambert DL (1996) Cosmic Ray—induced Chemistry 


toward Perseus OB2. Astrophys J 463:181. https: //doi.org/10.1086/177233 

Ferrière KM (2001) The interstellar environment of our galaxy. Reviews 
of Modern Physics 73(4):1031-1066. https: //doi.org/10.1103/RevModPhys. 
astro-ph 

Field GB, Goldsmith DW, Habing HJ (1969) Cosmic-Ray Heating of the In- 
terstellar Gas. Astrophys J Lett 155:L149. https://doi.org/10.1086/180324 

Fisk LA, Forman MA, Axford WI (1973) Solar modulation of galactic cosmic 
rays: 3. Implications of the Compton-Getting Coefficient. J Geophys Res 
78(7):995. https://doi.org/10.1029/JA078i007p00995 

Forman MA (1970) The Compton-Getting effect for cosmic-ray particles and 
photons and the Lorentz-invariance of distribution functions. Planet Space 
Sci 18(1):25-31. https: //doi.org/10.1016/0032-0633(70)90064-4 

Fornieri O, Gaggero D, Cerri SS, De La Torre Luque P, Gabici S (2021) The 
theory of cosmic ray scattering on pre-existing MHD modes meets data. Mon 


Not Royal Astron Soc 502(4):5821-5838. https://doi.org/10.1093/mnras/ 
stab355\ arXiv:2011.09197) [astro-ph.HE] 


Fuente A, Cernicharo J, Roueff E, et al. (2016) Ionization fraction and the 
enhanced sulfur chemistry in Barnard 1. Astron Astrophys 593:A94. |h 


//doi.org/10.1051/0004-6361/201628285, arXiv:1605.04724 [astro-ph.GA] 


Gabici S, Aharonian FA, Blasi P (2007) aeos odrat rays Sas molecular 
clouds. Astrophys Space Sci a. 4): ve 371. https://doi.org/10.1007/ 
310509-007-9427-6) jarXiv:astro-ph/0610032) [astro-ph] 


Gabici S, Aharonian FA, Casanova S (2009) Broad-band non-thermal emission 
from molecular clouds illuminated by cosmic rays from gee supernova 
mtemi Mon Not E DONE Ge Lie Astron Soc ee aS BI act 1639. 

copa S, Evoli C, ae ee Gee D, et mI (2019) te origin A Galactic 
cosmic rays: Challenges to the HRS es International Jour- 


nal of Modern Physics D 28(15):1930022-339. https://doi.org/10.1142/ 
$0218271819300222) arXiv:1903.11584) [astro-ph.HE] 


Gaches Se ini Ome Offner Se Hee TC Bisbas TG ( So The ae Impact of 
Cosmic Rays in Protoclusters. I. Molecular Cloud Chemistry. Astrophys 
J 878(2):105. https: //doi.org/10.3847/1538-4357 /ab20c7| |arXiv:1905.02232| 
oes a GA] 

Gaches BAL, Bisbas TG, Bialy S (2021la) The impact of cosmic-ray atten- 
uation on the carbon cycle emission in molecular clouds. arXiv e-prints 
arXiv:2110.06958. [astro-ph.GA] 

Gaches BAL, Walch S, Lazarian A (2021b) Cosmic-Ray Acceleration from 
Turbulence in Molecular Clouds. Astrophys J Lett 917(2):L39. 
[astro-ph HE) 

Gaisser TK, Engel R, Resconi E (2016) Cosmic Rays and Particle Physics. 
Cambridge University Press 

Geballe TR, Oka T (1996) Detection of H*3 in interstellar space. Nature 


Low energy cosmic rays 83 


384(6607):334-335. https://doi.org/10.1038/384334a0 


Génolini Y, Maurin D, Moskalenko IV, Unger M (2018) Current status and 
desired precision of the isotopic production cross sections relevant to astro- 
physics of cosmic rays: Li, Be, B, C, and N. Phys Rev C 98(3):034611. 
//doi-org/10.1103/PhysRevC.98.034611 [astro-ph.HE] 

Génolini Y, Boudaud M, Batista PI, et al. (2019) Cosmic-ray transport from 
AMS-02 boron to carbon ratio data: Benchmark models and interpretation. 


Phys Rev D 99(12):123028. {https://doi.org/10.1103/PhysRevD.99.123028 
arXiv:1904.08917) [astro-ph.HE] 


Gerin M, de Luca M, Black J, et al. (2010) Interstellar OH+, HOF and 

H30+ along the sight-line to G10.6-0.4. Astron Astrophys 518:L110. https: ‘https: 
//doi.org/10.1051/0004-6361/201014576, |arXiv:1005.5653 | [astro-ph.GA] 

Ghelfi A, Barao F, Derome L, Maurin D (2016) Non-parametric as 
tion of H and He interstellar fluxes from cosmic-ray data. Astron Astrophys 
591:A94. |https://doi.org/10.1051/0004-6361/201527852) jarXiv:1511.08650| 
[astro-ph.HE] 

Gibb EL, AM SD, Rettig TW, et al. (2010) CO and H*3 Toward MWC 
1080, MWC 349, and LkHa 101. Astrophys J 715(2):757—766. 
org/10.1088/0004-637X/715/2/757 

Gieseler J, Heber B, Herbst K (2017) An Empirical Modification of the 
Force Field Approach to Describe the Modulation of Galactic Cosmic Rays 
Close to Earth in a Broad Range of Rigidities. Journal of Geophysical 
Research (Space Physics) 122(11):10,964-10,979. https://doi.org/10.1002/ 
2017JA024763) ar Xiv:1710-10834]| [physics space-ph 

Glassgold AE, Langer WD (1974) Model oouo for diffuse molecular 
clouds. Astrophys J 193: a 91. |https://doi.org/10.1086/153130. 

Glassgold AE, Galli D, Padovani M (2012) Cosmic-Ray and X-Ray Heating of 
Interstellar Clouds and De o Astrophys J 756(2):157. 


//doi.org/10.1088/0004-637X/756/2/157, arXiv:1208.0523| [astro-ph.GA] 


Gleeson LJ, Axford WI (1967) Cosmic Rays in the a peno Medium. 


Astrophys J Lett 149:L115. https://doi.org/10.1086/180070 


Gleeson LJ, Axford WI (1968) The Compton-Getting Effect. Astrophys Space 


Sci 2(4):431-437. |https:/ /doi.org/10.1007 /BF02175919 

Gloeckler G, Fisk LA (2015) More Evidence that Voyager 1 Is Still in the He- 
liosphere. Astrophys J Lett 806(2):L27.|https://doi.org/10.1088/2041-8205/ 

Guelin M, Langer WD, Snell RL, Wootten HA (1977) Observations of DCO*: 
the electron abundance in dark clouds. Astrophys J Lett 217:L165-L168. 
https: //doi.org/10.1086/182562 

Hanasz M, Strong AW, Girichidis P (2021) Simulations of cosmic ray propa- 
gation. Living Reviews in Computational Astrophysics 7(1):2. 


org/10.1007/s41115-021-00011-1, jarXiv:2106.08426) [astro-ph.HE] 


Hayakawa S (1952) Propagation of the Cosmic Radiation through Intersteller 
Space. Progress of Theoretical Physics 8(5):571-572. https://doi.org/10. 
1143/ptp/8.5.571 


Hayakawa S, Ito K, Terashima Y (1958) Origin of Cosmic Rays. Progress of 


84 Stefano Gabici 


Theoretical Physics Supplement 6:1-92. https: //doi.org/10.1143/PTPS.6.1 


Hayakawa S, Nishimura S, Takayanagi T (1961) Radiation from the Interstellar 
Hydrogen Atoms. Publ Astron Soc Jpn 13:184 
Hennebelle P, Falgarone E (2012) Turbulent molecular clouds. As- 


tron Astrophys Rev 20:55. https://doi.org/10.1007/s00159-012-0055-y 
[astro-ph.GA] 
Herbst E, Klemperer W (1973) The Formation and Depletion of Molecules 
in Dense Interstellar Clouds. Astrophys J 185:505-534. 
Heyer M, Dame TM (2015) Molecular Clouds in the Milky Way. 
Ann Rev Astron Astrophys  53:583-629. https://doi.org/10.1146/ 
Heywood I, Rammala I, Camilo F, et al. (2022) The 1.28 GHz MeerKAT 
Galactic Center Mosaic. Astrophys J 925(2):165. |https://doi.org/10.3847/ 


1538-4357/ac449a, arXiv:2201.10541) [astro-ph.GA] 


Hezareh T, Houde M, McCoey C, Vastel C, Peng R (2008) Simultaneous Deter- 
mination of the Cosmic Ray Ionization Rate and Fractional Ionization in DR 


21(OH). Astrophys J 684(2):1221-1227. https://doi.org/10.1086/590365 
ar Xiv:0805.4018) [astro-ph] 


Hollenbach D, Kaufman MJ, Neufeld D, Wolfire M, Goicoechea JR 
(2012) The Chemistry of Interstellar OH*t, H2O+, and H3O?t: Infer- 
ring the Cosmic-Ray Ionization Rates from Observations of Molecular 


Ions. Astrophys J 754(2):105. https://doi.org/10.1088/0004-637X/754/2/ 
arXiv:1205.6446) [astro-ph.GA] 


Hollweg JV (1974) Hydromagnetic Waves in Interplanetary Space. Pub Astron 


Soc Aust 86(513):561. |https://doi.org/10.1086/129646 
Hsieh KC, Mason GM, Simpson JA (1971) Cosmic-Ray {2}H from Satellite 
Measurements, 1965-1969. Astrophys J 166:221. https://doi.org/10.1086/ 
Indriolo N, McCall BJ (2012) Investigating the Cosmic-Ray Ionization Rate 
in the Galactic Diffuse Interstellar Medium through Observations of Ht 
3. Astrophys J 745(1):91. |https://doi.org/10.1088/0004-637X/745/1/91 
[astro-ph.GA] 


Indriolo N, Geballe TR, Oka T, McCall BJ (2007) H*3 in Diffuse Inter- 
stellar Clouds: A Tracer for the Cosmic-Ray Ionization Rate. Astrophys J 
671(2):1736-1747. |https://doi.org/10.1086/523036 [astro- 
ph] 

Indriolo N, Blake GA, Goto M, et al. (2010) Investigating the Cosmic-ray Ion- 
ization Rate Near the Supernova Remnant IC 443 through H* 3 Observa- 
n E ei J 724( E ls 1365. T //doi.org/10.1088/0004-637X/ 


MN N, N EN we DAs Gea Gain 7 et al Bi) Chemical Analysis of a Dif- 
fuse Cloud along a Line of Sight toward W51: Molecular Fraction and 


Cosmic-Ray Ionization Rate. Astrophys J 758(2):83. https://doi.org/10. 


1088 /0004-637X /758/2/83] {arXiv:1208.3202) [astro-ph.GA] 


Indriolo N, Neufeld DA, Gerin M, et al. ae Herschel i of Galac- 


Low energy cosmic rays 85 


tic OH+, HOt, and H30O*: Probing the Molecular Hydrogen Fraction 
and Cosmic-Ray Ionization Rate. Astrophys J 800(1):40. |https://doi.org/ 


10.1088 /0004-637X/800/1/40) {arXiv:1412.1106)| [astro-ph.GA] 


Issa MR, Wolfendale AW (1981) y Rays from the cosmic ray irradiation of 
local molecular clouds. Nature 292(5822):430—433. https://doi.org/10.1038/ 

Ivlev AV, Padovani M, Galli D, Caselli P (2015) Interstellar Dust Wee in 
Dense Molecular Clouds: Cosmic Ray Effects. Astrophys J 812(2):135.\h 

//doi.org/10.1088/0004-637X/812/2/135 peer woo wee ph.GA] 

Ivlev AV, Dogiel VA, Chernyshov DO, et al. (2018) Penetration of Cosmic 
Rays into Dense Molecular Clouds: Role of Diffuse Envelopes. Astrophys 
J 855(1):23. https: //doi.org/10.3847/1538-4357 /aaadb9| jarXiv:1802.02612| 
[astro- a HE] 

Ivlev AV, cee K, Padovani M, Galli D (2021) Rigorous Theory for Sec- 
ondary Cosmic-Ray Ionization. Astrophis J 909(2):107. 


3847 /1538-4357/abdc27, arXiv:2101.05803) [astro-ph.GA] 


Jokipii JR, Parker EN (1967) Energy changes of cosmic rays in the so- 
lar system. Planet Space Sci 15(9):1375-1386. https://doi.org/10.1016/ 
0032-0633(67)90111-0 

Jokipii JR, Parker EN (1970) on the Convection, Diffusion, and Adiabatic 
Deceleration of Cosmic Rays in the Solar Wind. Astrophys J 160:735. 
//doi.org/10.1086/150465 

Jones DI (2014) Prospects for Detection of Synchrotron Emission from 
Secondary Electrons and Positrons in Starless Cores: Application to 


G0.216+0.016. Astrophys J Lett 792(1):L14. |https://doi.org/10.1088/ 
2041-8205/792/1/L14, jarXiv:1408.1822) [astro-ph.HE] 


Jones DI, Protheroe RJ, Crocker RM ( aa aoe) Search a Synchrotron Emis- 
sion from Secondary Leptons in Dense Cold Starless Cores. Publ Astron 
Soc Aust 25(4):161-166. https: //doi-org/10.1071/AS08007 
[astro-ph] 

Jones DI, Crocker RM, Ott J, Protheroe RJ, Ekers RD (2011) Australia Tele- 
scope Compact Array Radio Continuum 1384 and 2368 MHz Observations 


of Sagittarius B. Astron J 141(3):82. |https://doi.org/10.1088/0004-6256 / 
141/3/82\ jarXiv:1011.1656] [astro-ph.GA] 


Jones FC, Lukasiak A, Ptuskin V, Webber W (2001) The Modified Weighted 
a RRI Models ae Results. Se STIG i aie i T 
pae P Dar i re ee eee tae R (2017) Does maero paan aN rate a the very 
high energy cosmic rays in an mae 200 pe of our Galaxy? Mon Not 


Royal Astron Soc 467(4):4622—4630. https: //doi.org/10.1093/mnras/stx361 
ar Xiv:1703.10398) [astro-ph.HE] 


Kafexhiu E, Aharonian F, Taylor AM, Vila GS (2014) Parametrization 
of gamma-ray production cross sections for p p interactions in a broad 
proton energy range from the kinematic threshold to PeV energies. 


Phys Rev D 90(12):123014. {https: //doi.org/10.1103/PhysRevD.90.123014 
arXiv:1406.7369) [astro-ph.HE] 


86 Stefano Gabici 


Kamae T, Karlsson N, Mizuno T, Abe T, Koi T (2006) Parameterization of 
y, et/~, and Neutrino Spectra Produced by p-p Interaction in Astronom- 
ical Environments. Astrophys J 647(1):692—708. https://doi.org/10.1086/ 
[astro-ph 


Kelner SR, Aharonian FA, Bugayov VV (2006) Energy spectra of gamma 
rays, electrons, and neutrinos produced at proton-proton interactions in 
the very high energy regime. Phys Rev D 74(3):034018. 
astro-ph 

Kim YK, Rudd ME (1994) Binary-encounter-dipole model for electron- 
impact ionization. Phys Rev A 50(5):3954-3967. https://doi.org/10.1103/ 


Kim YK, Santos JP, Parente F (2000) Extension of the binary-encounter- 
dipole model to relativistic incident electrons. Phys Rev A 62(5):052710. 
https://doi.org/10.1103/PhysRevA.62.052710 

Kirk JG (1994) Particle Acceleration (With 26 figures). In: Saas-Fee Advanced 
Course 24: Plasma Astrophysics. p 225 

Koldobskiy S, Kachelrie8 M, Lskavyan A, et al. (2021) Energy spectra of 
secondaries in proton-proton interactions. arXiv e-prints arXiv:2110.00496. 
astro-ph. HE] 

Krause J, Morlino G, Gabici S (2015) CRIME - cosmic ray interactions 
in molecular environments. In: 34th International Cosmic Ray Conference 
(ICRC2015). International Cosmic Ray Conference, vol 34. p 518 

Krimigis SM, Armstrong TP, Axford WI, et al. (1977) The Low Energy 
Charged Particle (LECP) Experiment on the Voyager Spacecraft. Space 
Sci Rev 21(3):329-354. |https:/ /doi.org/ 10.1007 /BF00211545 

Krimigis SM, Decker RB, Roelof EC, et al. (2013) Search for the Exit: Voy- 
ager 1 at Heliosphere’s Border with the Galaxy. Science 341(6142):144-147. 
https://doi.org/10.1126/science.1235721 

Krimigis SM, Decker RB, Roelof EC, et al. (2019) Energetic charged parti- 
cle measurements from Voyager 2 at the heliopause and beyond. Nature 
Astronomy 3:997—1006. https://doi.org/10.1038/s41550-019-0927-4 

Krumholz MR (2014) The big problems in star formation: The star formation 
rate, stellar clustering, and the initial mass function. Phys Rep 539:49- 


E Ni https: //doi.org/10.1016/j-physrep.2014.02.001, arXiv:1402.0867| [astro- 


eo GF (1977) A regular mechanism for the acceleration of charged parti- 
cles on the front of a shock wave. Akademiia Nauk SSSR Doklady 234:1306- 
1308 

Kuhlen M, Mertsch P (2019) Time-Dependent AMS-02 Electron-Positron 
Fluxes in an Extended Force-Field Model. Phys Rev Lett 123(25):251104. 
https: //doi.org/10.1103/PhysRevLett.123.251104) |arXiv:1909.01154| [astro- 
ph.HE] 

Kulsrud B Pearce WP (1969) The Effect of Wave-Particle Interactions on 


the Propagation of Cosmic Rays. Astrophys J 156:445. https://doi.org/10. 
1086/149981 


Kulsrud RM (2005) Plasma physics for astrophysics 


Low energy cosmic rays 87 


Kulsrud RM, Cesarsky CJ (1971) The Effectiveness of Instabilities for the 
Confinement of High Energy Cosmic Rays in the Galactic Disk. Astrophys 
Lett 8:189 

Larsson M (2000) Experimental studies of the dissociative recombination of 
H3+t. In: Astronomy, physics and chemistry of H*3. vol 3. pp 2433-2444. 
https: //doi.org/10.1098/rsta.2000.0658 

Larsson M, Geppert WD, Nyman G (2012) Ion chemistry in space. Reports on 
Progress in Physics 75(6):066901. https: //doi.org/10.1088/0034-4885 /75/6/ 
[066901] 

Le Petit F, Ruaud M, Bron E, et al. (2016) Physical conditions in the central 
molecular zone inferred by H3t. Astron Astrophys 585:A105. 
iastro-ph.Ga) 


Lee HH, Bettens RPA, Herbst E (1996) Fractional abundances of molecules 
in dense interstellar clouds: A compendium of recent model results. Astron 
Astrophys Suppl 119:111-114 

Lerche I (1967) Unstable Magnetosonic Waves in a Relativistic Plasma. As- 
trophys J 147:689. https: //doi.org/10.1086/149045 

Lindsay CM, McCall BJ (2001) Comprehensive Evaluation and Comipilarion of of 
H 3+ Spectroscopy. Journal of Molecular Spectroscopy 210(1):60-83. 
//doi.org/10.1006 /jmsp.2001.8444 

Lingenfelter RE (2018) Cosmic rays from supernova remnants and super- 
bubbles. Advances in Space Research 62(10):2750-2763. 


1016/j.asr.2017.04.006) |arXiv:1807.09726 [astro-ph.HE] 


Liszt H (2003) Gas-phase recombination, grain neutralization and cosmic-ray 
ionization in diffuse gas. Astron Astrophys 398:621—630. https: //doi.org/10. 


1051/0004-6361:20021660, arXiv:astro-ph/0211369) [astro-ph] 


Luhmann JG (1976) A quasi-linear kinetic equation for cosmic rays in the in- 
ee ae J Geophys Res 81(13):2089. https: //doi.org/10.1029/ 
JA081i013p02089 

Mac Low MM, Klessen RS (2004) Control of star formation by supersonic 


turbulence. Reviews of Modern Physics 76(1):125-194. https://doi.org/10. 
1103/RevModPhys.76.125) jarXiv:astro-ph/0301093) [astro-ph] 


Marcowith A, van Marle AJ, Plotnikov I (2021) The cosmic ray-driven stream- 
ing instability in astrophysical and space plasmas. Physics of Plasmas 
28(8):080601. |https://doi.org/10.1063/5.0013662 

Maret S, Bergin EA (2007) The Ionization Fraction of Barnard 68: Implications 
for Star and Planet Formation. Astrophys J 664(2):956-963. 
(astro-pl 


Marscher AP, Brown RL (1978) Energetic secondary electrons and the non- 
thermal galactic radio background: a probe of the magnetic field in inter- 
stellar clouds. Astrophys J 221:588-597. https://doi.org/10.1086/156061 

McCall BJ, Geballe TR, Hinkle KH, Oka T (1999) Observations of Ht; in 
Dense Molecular Clouds. Astrophys J 522(1):338-348. 


1086/307637 


McCall BJ, Hinkle KH, Geballe TR, et al. (2002) Observations of H3* in the 
Diffuse Interstellar Medium. Astrophys J 567(1):391—406. |https://doi.org/ 


88 Stefano Gabici 


10.1086 /338380) arXiv:astro-ph/0110674 [astro-ph] 


McCall BJ, Huneycutt AJ, Saykally RJ, et al. (2003) An enhanced cosmic- 
ray flux towards ¢ Persei inferred from a laboratory study of the H3*- 
e7 recombination rate. Nature 422(6931):500—502. https://doi.org/10.1038/ 
inatre(1 408) arXiv:astro-ph/0302106 [astro-ph 


McCall BJ, Huneycutt AJ, Saykally RJ, et al. (2004) Dissociative recombina- 
tion of rotationally cold Ht3. Phys Rev A 70(5):052716. https: //doi.org/10. 
1103/PhysRevA.70.052716 


McEnery J, van der Horst A, Dominguez A, et al. (2019) All-sky Medium 
Energy Gamma-ray Observatory: Exploring the Extreme Multimessenger 
Universe. In: Bulletin of the American Astronomical Society. vol 51. p 245. 


arXiv:1907.07558) [astro-ph.IM] 


McKee CF (1989) Photoionization-regulated Star Formation and the Structure 


of Molecular Clouds. Astrophys J 345:782. https: //doi.org/10.1086/167950 
McKee CF, Ostriker EC (2007) Theory of Star Formation. Ann Rev As- 
aie BRIBE TOSS 45( an 565- ae eae ee org/10.1146/annurev.astro.45. 


u M, eA A TE J, Ap o The production of the elements 
Li, Be, B by galactic cosmic rays in space and its relation with stellar ob- 
servations. Astron Astrophys 15:337 

Meng F, Sdnchez-Monge A, Schilke P, et al. (2019) The physical and chem- 
ical structure of Sagittarius B2. V. Non-thermal emission in the envelope 


of Sgr B2. Astron Astrphys 630:A73. https://doi-org/10.1051/0004-6361/ 
201935920) |arXiv:1908.07237 [astro-ph.GA| 


Menn W, RaW Heth Rome Oo M, Reimer O, et i (2000) ie Absolute Flux of Protons and 
Helium at the Top of the ‘Abropaphere Using IMAX. Astrophys J 533(1):281-— 


297. |https: //doi.org/10.1086/308645 


Mertsch P (2020) Test particle simulations of cosmic rays. Astrophys Space Sci 

365(8):135. |https: //doi.org/10.1007/s10509-020-03832-3) arXiv:1910.01172) 

jones ph.HE] 

Miller S, ee J, Geballe TR, Stallard T (2020) Thirty years of H3* 

astronomy. Reviews of Modern Physics 92(3):035003. 

Moraal H (2013) Cosmic-Ray Modulation Equations. Space Sci Rev 176(1- 

4):299-319. https: //doi.org/10.1007/s11214-011-9819-3 

Morales Ortiz JL, Ceccarelli C, Lis DC, et al. (2014) Ionization to- 
ward the high-mass star-forming region NGC 6334 I. Astron Astrophys 
563:A 127. |https://doi.org/10.1051/0004-6361/201322071) |arXiv:1306.3012| 
[astro-ph.GA] 

Morfill GE F Gamma ray emission from interstellar clouds - A plasma 

physical process capable of enhancing electron fluxes. Mon Not Royal Astron 

Soc 198:583-588. https://doi.org/10.1093/mnras/198.2.583 

Morfill GE (1982b) Secondary electron spectra in interstellar clouds, and the 

bremsstrahlung gamma-ray luminosity. Astrophys J 262:749-759. 


doi.org/10.1086/160470 


Mori M (2009) Nuclear enhancement factor in calculation of Galac- 


Low energy cosmic rays 89 


tic diffuse gamma-rays: A new estimate with DPMJET-3. Astroparti- 
cle Physics 31(5):341-343. https: //doi.org/10.1016/j.astropartphys.2009.03. 
astro-ph HE] 

Morlino G, Gabici S (2015) Cosmic ray penetration in diffuse clouds. Mon Not 
Royal Astron Soc 451:L100—L104. |https://doi.org/10.1093/mnrasl/slv074 
astro-ph. HE] 


Müller HSP, Endres CP, Stutzki J, Schlemmer S (2013) The CDMS view on 
molecular data needs of Herschel, SOFIA, and ALMA. In: Gillaspy JD, 
Wiese WL, Podpaly YA (eds) Eighth International Conference on Atomic 
and Molecular Data and Their Applications: ICAMDATA-2012. American 


Institute of Physics Conference Series, vol 1545. pp 96-109. https://doi.org/ 


10.1063/1.4815844) |arXiv:1304.1894 [astro-ph.GA] 


Neufeld DA, Wolfire MG ( TG The on Ray Ionization Rate in the Galac- 
tic Disk, as Determined from Observations of Molecular Ions. Astrophys 
J 845(2):163. |https: //doi.org/10.3847 /1538-4357 /aa6d68) arXiv:1704.03877, 
es - GA] 

Neufeld DA, Goicoechea JR, Sonnentrucker P, et al. (2010) Herschel/HIFI ob- 
servations of interstellar OH+ and H20O* towards W49N: a probe of diffuse 
clouds with a small molecular fraction. Astron Astrophys 521:L10. [https:] 

//doi.org/10.1051/0004-6361/201015077, [arXiv:1007.0987]| [astro-ph.GA] 

Nobukawa KK, Nobukawa M, Koyama K, et al. (2018) Evidence for a Neutral 

Iron Line Generated by MeV Protons from Supernova Remnants Interact- 


ing with Molecular Clouds. Astrophys J 854(2):87. https://doi.org/10.3847/ 
1538-4357 /aaa8dc; arXiv:1801.07881) [astro-ph.HE] 


Nozzoli F, Cernetti C (2021) Beryllium Radioactive Isotopes as a Probe to 
Measure the Residence Time of Cosmic Rays in the Galaxy and Halo Thick- 


ness: A “Data-Driven” Approach. Universe 7(6):183. https://doi.org/10. 
3390/universe7060183 


Oka T (2003) Help!!! Theory for H3+ Recombination Badly Needed. In: Dis- 
sociative Recombination of Molecular Ions with Electrons. p 209 

Oka T (2006) Interstellar Chemistry Special Feature: Interstellar H3+. Pro- 
ceedings of the National Academy of Science 103(33):12235-12242. 
//doi.org/10.1073/pnas.0601242103 

Oka T, Geballe TR, Goto M, et al. (2019) The Central 300 pc of the Galaxy 
Probed by Infrared Spectra of Hf and CO. I. Predominance of Warm and 
Diffuse Gas and High Hə Ionization Rate. Astrophys J 883(1):54 

//doi.org/10.3847/1538-4357 /ab3647, arXiv:1910.04762 [astro-ph. HE] 

Orlando E (2018) Imprints of cosmic rays in multifrequency observations i the 
eerie emission. Mon Not Royal Astron Soc 475(2):2724—2742. 
//doi.org/10.1093 /mnras/stx3280 [astro-ph.HE] 

Orlando E, Strong A (2013) Galactic synchrotron emission with cosmic ray 
propagation models. Mon Not Royal Astron Soc 436(3):2127-2142. 
//doi.org/10.1093/mnras/stt1718 [astro-ph.GA] 

Orlando E, Bottacini E, Moiseev A, et al. (2021) Exploring the MeV Sky with 
a Combined Coded Mask and Compton Telescope: The Galactic Explorer 
with a Coded Aperture Mask Compton Telescope (GECCO). arXiv e-prints 


90 Stefano Gabici 


arXiv:2112.07190. arXiv:2112.07190 [astro-ph.HE] 
Padovani M, Galli D (2018) Synchrotron emission in molecular cloud cores: the 


SKA view. Astron Astrophys 620:L4. https://doi.org/10.1051/0004-6361/ 
201834222) arXiv:1811.09698 [astro-ph.HE] 


Padovani M, Galli D, Glassgold AE (2009) Cosmic-ray ionization of molec- 
ular clouds. Astron Astrophys 501(2):619-631. https://doi.org/10.1051/ 
0004-6361/200911794, arXiv:0904.4149 [astro-ph.SR] 


Padovani M, Hennebelle P, Galli D ( aD aa) eae ray ionisation in collapsing 


clouds. Astron Astrophys 560:A114. https://doi.org/10.1051/0004-6361/ 
201522407, jarXiv:1510.2158 [astro-ph GA 
Padovani M, Marcowith A, Hennebelle P, Ferriére K (2016) Protostars: 
TOOL BST OA of cosmic 2u Astron TEON 590:A8. https://doi.org/10.1051/ Pipa.) 
pcan M, Ivlev ee coat fied Cael D e won Gand ray ionisation in 
ea discs. a oo m A111. Pepe: /Paok org/79.1061/ 


pecan M, Ivlev eee ae A Galli D, et al sae oe a aoa) a a Low-Energy Cosmic 
Rays on Star Formation. Space Sci Rev 216(2):29. https: //doi.org/10.1007/ 
$11214-020-00654- 1} arXiv:2002.10282 [astro-ph.GA] 


Padovani M, Bialy S, Galli D, et al. (2022) Cosmic rays in molecular clouds 
probed by Hə rovibrational lines — Perspectives for the James Webb Space 
Telescope. arXiv e-prints arXiv:2201.08457. [astro-ph.GA] 

Parker EN (1965) The passage of energetic charged particles through in- 
terplanetary space. Planet Space Sci 13(1):9-49. https://doi.org/10.1016/ 
0032-0633(65)90131-5 

Peron G, Aharonian F, Casanova S, Yang R, Zanin R (2021) Prob- 
ing the Cosmic-Ray Density in the Inner Galaxy. Astrophys J Lett 
907(1):L11. [https://doi.org/10.3847/2041-8213/abcaa9| [arXiv:2101.09510] 
n. ph.HE] 

Phan VHM con Cosmic ray interaction in molecular environment. PhD 
thesis, Université de Paris 

Phan VHM, Morlino G, Gabici S (2018) What causes the ionization rates 
observed in diffuse molecular clouds? The role of cosmic ray protons and 
aL Mon Not TS ae aS a on Soc EAA ):5167- ae 

PE eee T VHM, Gabici S, a G, et a Fee a aie the cosmic ray 
spectrum in the vicinity of the supernova remnant W28: from sub-GeV 


to multi-TeV energies. Astron Astrophys 635:A40. |https://doi.org/10.1051/ 
0004-6361/201936927, arXiv:1910.09987  [astro-ph.HE] 


Phan VHM, Schulze F, See aetes P, lea S, cae S (2021) Stochas- 
tic Fluctuations of Low-Energy Cosmic Rays and the Interpretation of 
Voyager Data. Phys Rev Lett 127(14):141101. https://doi.org/10.1103/ 
PhysRevLett 127 141701) arXiv:2103.00811] [astro-ph HE 


Pinto C, Galli D (2008) Three-fluid plasmas in star formation. II. Momentum 
transfer rate coefficients. Astron Astrophys 484(1):17—28. |https://doi.org/ 


10.1051 /0004-6361:20078819, arXiv:0804.2856 [astro-ph] 


Low energy cosmic rays 91 


Ponti G, Terrier R, Goldwurm A, Belanger G, Trap G (2010) Discovery of 
a Superluminal Fe K Echo at the Galactic Center: The Glorious Past of 
Sgr A* Preserved by Molecular Clouds. Astrophys J 714(1):732-747. 


//dci.org/10.1088/0004-637X /714/1/732) |arXiv:1003.2001) [astro-ph.HE] 


Porras AJ, Federman SR, Welty DE, Ritchey a (2014) on in N 
Molecular Clouds. Astrophys J Lett 781(1):L8. KA 7 /doi.org/10.1088/ 
2041-82057781/1/L8| orXiv:1312.5239| astro-ph. GA] 

Potgieter MS (2013) Solar Modulation of Cosmic Rays. Living Reviews in So- 
lar Physics 10(1):3. [https://doi.org/10.12942/lrsp-2013-3| [arXiv:1306.4421] 
[physics.space- a 

Pothast M, Gaggero D, Storm E, Weniger C (2018) On the progressive hard- 
ening of the cosmic-ray proton spectrum in the inner Galaxy. J Cosm 


Astropart Phys 2018(10):045. https://doi.org/10.1088/1475-7516/2018/10/ 
arXiv:1807.04554| [astro-ph. HE] 


Protheroe RJ, Ott J, Ekers RD, Jones DI, Crocker RM (2008) Interpretation 
of radio continuum and molecular line observations of Sgr B2: free-free and 
synchrotron emission, and implications for cosmic rays. Mon Not Royal As- 


tron Soc 390(2):683-692. https: //doi.org/10.1111/}.1365-2966.2008.13752.x 
ar Xiv:0807.0127) [astro-ph] 


Ptuskin VS, Soutoul A (1998) Cosmic-ray Clocks. Space Sci Rev 86:225-238. 


https: //doi.org/10.1023/A:1005071612125 


Quenby JJ (1984) The Theory of Cosmic-Ray Modulation. Space Sci Rev 
37(3-4):201-267. [https: //doi.org/10.1007 /BF 00226364 

Ramaty R, Kozlovsky B, Lingenfelter RE (1979) Nuclear gamma-rays from 
energetic particle interactions. Astrophys J Suppl 40:487-526. 


org/10.1086/190596 


Reames DV (1999) Particle acceleration at the Sun and in the heliosphere. 


Space Sci Rev 90:413-491. https: //doi.org/10.1023/A:1005105831781 


Recchia S (2020) Cosmic ray driven galactic winds. International Jour- 
nal of Modern Physics D 29(7):2030006-50. https://doi.org/10.1142/ 
Recchia S, Phan VHM, Biswas S, Gabici S (2019) Can a cosmic ray car- 
rot explain the ionization level in diffuse molecular clouds? Mon Not 


Royal Astron Soc 485(2):2276—2280. https: //doi-org/10.1093/mnras/stz555 
arXiv:1901.04912) [astro-ph.HE] 


Recchia S, Galli D, Nava L, et al. (2021) Grammage of cosmic rays in the 
proximity of supernova remnants embedded in a partially ionized medium. 


arXiv e-prints arXiv:2106.04948. jarXiv:2106.04948 [astro-ph.HE] 


Reeves H (1994) On the origin of the light elements (Zj6). Reviews of Modern 


Physics 66(1):193-216. |https://doi.org/10.1103/RevModPhys.66.193 
Reeves H, Fowler WA, Hoyle F (1970) Galactic Cosmic Ray Origin of Li, 
Be and B in Stars. Nature 226(5247):727-729. https://doi.org/10.1038/ 
[226727a0] 
Remy Q, Grenier IA, Marshall DJ, Casandjian JM (2017) Cosmic rays, gas and 
dust in nearby anticentre clouds. I. CO-to-Hə conversion factors and dust 


opacities. Astron Astrophys 601:A78. |https://doi.org/10.1051/0004-6361/ 


92 Stefano Gabici 


201629632, ar Xiv:1703.05237 [astro-ph.HE] 


Reville B, Giacinti G, Scott R (2021) Cosmic-ray current-driven in- 
stabilities - revisiting environmental conditions. Mon Not Royal 


Astron Soc 502(3):4137-4153. https://doi.org/10.1093/mnras/stab296 
ar Xiv:2102.00334)| [astro-ph.HE] 


Rimmer PB, Herbst E, Morata O, Roueff E (2012) Observing a column- 
dependent ¢ in dense interstellar sources: the case of the Horsehead nebula. 


Astron Astrophys 537:A7. |https://doi.org/10.1051/0004-6361/201117048 
arXiv:1110.2399| [astro-ph.GA] 


Rodriguez LF, Zapata LA n Star Formation in the Massive “Starless” 
Infrared Dark Cloud G0.253+0.016. Astrophys J Lett 767(1):L13. 
//doi.org/10.1088/2041-8205/767/1/L13 [astro-ph.SR] 

Rodriguez-Kamenetzky A, Carrasco-Gonzalez C, Araudo A, et al. (2017) The 
Highly Collimated Radio Jet of HH 80-81: Structure and Nonthermal Emis- 


sion. Astrophys J 851(1):16. https://doi.org/10.3847/1538-4357/aa9895 
arXiv:1711.02554| [astro-ph. HE] 


Rogers F, Zhang S, Perez K, Clavel M, Taylor A (2021) New Con- 
straints on Cosmic Particle Populations at the Galactic Center using X- 
ray Observations of the Molecular Cloud Sagittarius B2. arXiv e-prints 
arXiv:2108.13399. [astro-ph.HE] 

Rudd ME (1988) Differential cross sections for secondary electron production 
by proton impact. Phys Rev A 38(12):6129-6137. https://doi.org/10.1103/ 


PhysRevA.38.6129 
Rudd ME, Kim YK, Madison DH, Gay TJ (1992) Electron production in 


proton collisions with atoms and molecules: energy distributions. Reviews 
of Modern Physics 64(2):441—490. https://doi.org/10.1103/RevModPhys. 

Sabatini G, Bovino S, Giannetti A, et al. (2020) Survey of ortho-H2D* in high- 
mass star-forming regions. Astron Astrophys 644:A34. 
astro-ph. Gal 

Savage BD, Bohlin RC, Drake JF, Budich W (1977) A survey of interstellar 
molecular hydrogen. I. Astrophys J 216:291-307. https://doi.org/10.1086/ 

Scherer K, Fichtner H, Strauss RD, et al. (2011) On Cosmic Ray Modulation 
beyond the Heliopause: Where is the Modulation Boundary? Astrophys J 
735(2):128. https: //doi.org/10.1088/0004-637X/735/2/128 

Schilke P, Neufeld DA, Müller HSP, et al. (2014) Ubiquitous argonium (ArH*) 
in the diffuse interstellar medium: A molecular tracer of almost purely 


atomic gas. Astron Astrophys 566:A29. https: //doi.org/10.1051/0004-6361/ 
201423727, |arXiv:1403.7902  [astro-ph.GA|] 


Eoad R SF) Coa Cosmic hae oe 

Schlickeiser R, Webber WR, Kempf A (2014) Explanation of the Local Galactic 
Cosmic Ray Energy Spectra Measured by Voyager 1. I. Protons. Astrophys 
J 787(1):35. |https: //doi.org/10.1088/0004-637X/787/1/35 

Schlickeiser R, Caglar M, Lazarian A (2016) Cosmic Rays and MHD Tur- 
bulence Generation in Interstellar Giant Molecular Clouds. Astrophys J 


Low energy cosmic rays 93 


824(2):89. |https: //doi.org/10.3847/0004-637X /824/2/89 
Scoville NZ, Sanders DB (1987) Hə in the Galaxy. In: Hollenbach DJ, Thronson 
J Harley A (eds) Interstellar Processes. vol 134. p 21. 


Shaw G, Ferland GJ (2021) Role of Polycyclic Aromatic Hydrocarbons on the 

Cosmic-Ray Ionization Rate in the Galaxy. Astrophys J 908(2):138. 
//doi.org/10.3847/1538-4357 /abdab2| |arXiv:2101.03732| [astro-ph.GA] 

Shikaze Y, Haino S, Abe K, et al. (2007) Measurements of 0.2 20 vk 
cosmic-ray proton and helio spectra from 1997 through 2002 with the 
BESS spectrometer. Astroparticle Physics 28(1):154-167. 
n 

Shu FH, Adams FC, Lizano S (1987) Star formation in molecular clouds: 
observation and theory. Ann Rev Astron Astrophys 25:23-81. 

Silsbee K, Ivlev AV (2019) Diffusive versus Free-streaming Cosmic-Ray Trans- 
port in Molecular Clouds. Astrophys J 879(1):14. https://doi.org/10.3847/ 
astro-ph. HE] 


Silsbee K, Ivlev AV, Padovani M, Caselli P (2018) Magnetic Mirroring and 
Focusing of Cosmic Rays. Astrophys J 863(2):188. https://doi.org/10.3847/ 
1538-4997 acct eX 1807. 05025 [astro ph. HE] 

Skilling J (1975a) Cosmic ray streaming - I. Effect of = waves on particles. 

Mon Not Royal Astron Soc 172:557-566. https://doi.org/10.1093/mnras/ 

Skilling J (1975b) Cosmic ray streaming - III. Self-consistent solutions. Mon 

Not Royal Astron Soc 173:255-269. |https://doi.org/10.1093/mnras/173.2. 

Skilling J, Strong AW (1976) Cosmic ray exclusion from dense molecular 
clouds. Astron Astrophys 53(2):253-258 

Snow TP, McCall BJ (2006) Diffuse Atomic and Molecular Clouds. Ann Rev 
Astron Astrophys 44(1):367—414. https://doi.org/10.1146/annurev.astro.43. 
() () 0624 

Spitzer J Lyman, Tomasko MG (1968) Heating of H i Regions by Energetic 
Particles. Astrophys J 152:971. https://doi.org/10.1086/149610 

Stahler SW, Palla F (2004) The Formation of Stars 

Stecker FW (1971) Cosmic gamma rays, vol NASA-SP-249 

Stone EC, Vogt RE, McDonald FB, et al. (1977) Cosmic ray investiga- 
tion for the Voyager missions; energetic particle studies in the outer he- 
liosphere—And beyond. Space Sci Rev 21(3):355-376. 


Stone EC, Cummings AC, McDonald FB, et al. (2005) Voyager 1 Ex- 
plores the Termination Shock Region and the Heliosheath Beyond. Science 


309(5743):2017—2020. https: //doi.org/10.1126/science. 1117684 


Stone EC, Cummings AC, McDonald FB, et al. (2013) Voyager 1 Observes 
Low-Energy Galactic Cosmic Rays in a Region Depleted of Heliospheric 


Ions. Science 341(6142):150-153. https: //doi.org/10.1126/science. 1236408 


Stone EC, Cummings AC, Heikkila BC, Lal N (2019) Cosmic ray measure- 


94 Stefano Gabici 


ments from Voyager 2 as it crossed into interstellar space. Nature Astronomy 
3:1013-1018. |https: //doi.org/10.1038/s41550-019-0928-3 

Streitmatter RE, Stephens SA (2001) The abundance of radioactive nuclei 
under different physical environments around the solar system. Advances 
in Space Research 27(4):743-748. |https://doi.org/10.1016/S0273-1177(01) 
00115-6] 

Strong A (2015) Local interstellar cosmic-ray spectra derived from gamma-ray 
emissivities. In: 34th International Cosmic Ray Conference (ICRC2015). In- 
ternational Cosmic Ray Conference, vol 34. p 506. [astro- 
ph.HE] 

Strong AW (2011) Interstellar Gamma Rays and Cosmic Rays:. New Insights 
from Fermi-Lat and Integral. In: Giani S, Leroy C, Rancoita PG (eds) 
Cosmic Rays for Particle and Astroparticle Physics. pp 473-481. 
//doi.org/10.1142/9789814329033_0059 [astro-ph.HE] 

Strong AW, Moskalenko IV, Ptuskin VS (2007) Cosmic-Ray Propagation 
and Interactions in the Galaxy. Annual Review of Nuclear and Parti- 


cle Science 57(1):285-327. https://doi.org/10.1146/annurev.nucl.57.090506. 
123011, arXiv:astro-ph/0701517 [astro-ph] 


Strong AW, Porter TA, Digel SW, et al. (2010) Global Cosmic-ray- 
related Luminosity and Energy Budget of the Milky Way. Astro- 


phys J Lett 722(1):L58-L63. |https://doi.org /10.1088 /2041-8205/722/1/ 
ar Xiv:1008.4330) [astro-ph.HE] 


Strong AW, Orlando E, Jaffe TR (2011) The interstellar cosmic-ray elec- 
tron spectrum from synchrotron radiation and direct measurements. 


Astron Astrophys 534:A54. https://doi.org/10.1051/0004-6361/201116828 
arXiv:1108.4822) [astro-ph.HE] 


Tartari A, Zannoni M, Gervasi M, Boella G, Sironi G (2008) TRIS. II. The 
Diffuse Galactic Radio Emission at 6 = +42°. Astrophys J 688(1):32-42. 
https: //doi.org/10.1086/592135 [astro-ph] 

Tatischeff V, Gabici S (2018) Particle Acceleration by Supernova Shocks 
and Spallogenic Nucleosynthesis of Light Elements. Annual Review 
SSC TOIT ARA Nuclear and Particle sts I oi https: er org/10.1146/ 


TatischefF V, Doak N ee G eae Sena X- te from low- 
energy cosmic rays: applieatian to the 6.4 keV line emission from the 
aac cluster EAGER POE ue SR STOTT | 546:A88. https: //doi.org/10.1051/ pap aot org/70.1981/ 


Leki. eee V, g ao Duprat T aa S, eais S (2021) The 
origin of Galactic cosmic rays as revealed by their composition. Mon 
Not e A SI UU Soc aa a 1321- a https://doi.org/10.1093/mnras/ 

et Haar 3 a pa Sa and Stellar Energy. Reviews of 
Modern Physics 22(2):119-152. https://doi.org/10.1103/RevModPhys.22. 
119 

Tibaldo L, Gaggero D, Martin P (2021) Gamma Rays as Probes of Cosmic- 
Ray Propagation and Interactions in Galaxies. Universe 7(5):141. 


Low energy cosmic rays 95 


//doi.org/10.3390/universe7050141, jarXiv:2103.16423) [astro-ph.HE] 


Tielens AGGM (2008) Interstellar SS ek aromatic ee ais 
Ann Rev Astron Astrophys 46:289-337. |https://doi.org/10.1146/annurev. 
46.060407.14 
Tielens AGGM (2013) The molecular universe. Reviews of Modern Physics 
85(3):1021-1081. https://doi.org/10.1103/RevModPhys.85.1021 
Tomsick JA, Boggs SE, Zoglauer A, et al. (2021) The Compton Spectrometer 
and Imager Project for MeV Astronomy. arXiv e-prints arXiv:2109.10403. 


ar Xiv:2109.10403} [astro-ph.IM] 


Trimble V (1991) The origin and abundances of the chemical elements revis- 
ited. Astron Astrophys Rev 3(1):1—46. https: //doi.org/10.1007/BF00873456 

van der Tak FFS, van Dishoeck EF (2000) Limits on the cosmic-ray ion- 
ization rate toward massive young stars. Astron Astrophys 358:L79-L82. 
Pastro phi 

van der Tak FFS, van Dishoeck EF, Evans I Neal J, Blake GA (2000) Structure 
and Evolution of the Envelopes of Deeply Embedded Massive Young Stars. 


Astrophys J 537(1):283-303. https: //doi.org/10.1086/309011) jarXiv:astro-| 
ph/0001527 a rs 


van Dishoeck EF, Bergin EA (2020) Astrochemistry associated with planet 
formation. arXiv e-prints arXiv:2012.01472. [astro-ph. EP] 
van Dishoeck EF, Black JH (1986) Comprehensive Models of Diffuse Inter- 
stellar Clouds: Physical Conditions and Molecular Abundances. Astrophys 


J Suppl 62:109. |https: //doi.org/10.1086/191135 


van Dishoeck EF, Black JH (1989) Interstellar C 2, CH, and CN in Translucent 


Molecular Clouds. Astrophys J 340:273. https://doi.org/10.1086/167391 
Vangioni-Flam E, Cassé M, Audouze J (2000) Lithium-beryllium-boron: 

origin and evolution. Phys Rep 333:365-387. tips: / dot org/10.1016/ 

$0370-1573(00)00030-2 SOTO 1572 (00 00050 a 2007 Fall arXiv:astro-ph/9907171 [astro-ph] 


Vaupré S, eB P, Ceccarelli C, et al. (2014) Cosmic ray induced ion- 
isation of a molecular cloud shocked by the W28 supernova remnant. 


ae Astrophys 568:A50. |https://doi.org/10.1051/0004-6361/201424036, 
[astro-ph.GA] 

Vemerto 5, ier P (2016) sees of very high energy gamma rays in the 
Milky Way. Phys Rev D 94(6):063009. https: //doi.org/10.1103/PhysRevD. 
astro-ph. HE] 


von Rosenvinge TT, McDonald FB, Trainor JH, Webber WR (1979) The Mod- 
ulation of Galactic Cosmic Rays Between 1 and 17 AU. In: International 
Cosmic Ray Conference. International Cosmic Ray Conference, vol 12. p 
171 

Vos EE, Potgieter MS (2015) New Modeling of Galactic Proton Modulation 
during the Minimum of Solar Cycle 23/24. Astrophys J 815(2):119. 


//doi.org/10.1088/0004-637X /815/2/119 


Wakelam V, Loison JC, Herbst E, et al. (2015) The 2014 KIDA Network 
for Interstellar Chemistry. Astrophys J Suppl 217(2):20. |https://doi.org/ 


10.1088 /0067-0049/217/2/20) arXiv:1503.01594) [astro-ph.GA] 


Wallerstein G, Iben J Icko, Parker P, et al. (1997) Synthesis of the elements 


96 Stefano Gabici 


in stars: forty years of progress. Reviews of Modern Physics 69(4):995-1084. 


https: //doi.org/10.1103/RevModPhys.69.995 

Wardle M (2007) Magnetic fields in protoplanetary disks. Astro- 
phys Space Sci 311(1-3):35-45. |https: //doi.org/10.1007/s10509-007-9575-8 
astro-ph] 


Watson WD (1973) The Rate of Formation of Interstellar Molecules by Ion- 


Molecule Reactions. Astrophys J Lett 183:L17. https://doi.org/10.1086/ 


Webber WR, Golden RL, Stochaj SJ, Ormes JF, Strittmatter RE (1991) A 
Measurement of the Cosmic-Ray 2H and 3He Spectra and 2H/ 4He and 3He/ 


4He Ratios in 1989. Astrophys J 380:230. https: //doi.org/10.1086/170578 


Wentzel DG (1968) Hydromagnetic Waves Excited by Slowly Streaming Cos- 


mic Rays. Astrophys J 152:987. https://doi.org/10.1086/149611 


Wentzel DG (1972) Interstellar Motions: Minuet or Rock? Publ Astron Soc 


Pac 84(498):225. |https://doi.org/10.1086/129277 


Wentzel DG (1974) Cosmic-ray propagation in the Galaxy: collective effects. 
Ann Rev Astron Astrophys 12:71-96. https://doi-org/10.1146/annurev.aa. 
Wiedenbeck ME, Binns WR, Cummings AC, et al. (2007) An Overview of the 
Origin of Galactic Cosmic Rays as Inferred from Observations of Heavy Ion 
Composition and Spectra. Space Sci Rev 130(1-4):415-429. 
Wolfire MG, McKee CF, Hollenbach D, Tielens AGGM (2003) Neutral Atomic 
Phases of the Interstellar Medium in the Galaxy. Astrophys J 587(1):278- 
311. |https: //doi.org/10.1086/368016 [astro-ph] 
Woosley SE, Heger A, Weaver TA (2002) The evolution and explosion of mas- 
sive stars. Reviews of Modern Physics 74(4):1015-1071. 
Wootten A, Snell R, Glassgold AE (1979) The determination of electron abun- 
dances in interstellar clouds. Astrophys J 234:876-880. 


Yang R, Aharonian F, Evoli C (2016) Radial distribution of the diffuse y -ray 
emissivity in the Galactic disk. Phys Rev D 93(12):123007. 
astro-ph HE 

Yusef-Zadeh F, Hewitt JW, Wardle M, et al. (2013) Interacting Cosmic 
Rays with Molecular Clouds: A Bremsstrahlung Origin of Diffuse High- 
energy Emission from the Inner 2° x 1° of the Galactic Center. Astrophys 
J 762(1):33. https: //doi-org/10.1088/0004-637X/762/1/33 
[astro-ph.HE] 

Zank GP (1999) Interaction of the solar wind with the local interstellar 
medium: a theoretical perspective. Space Sci Rev 89:413-688. 


Zank GP (2014) Transport Processes in Space Physics and Astrophysics, vol 


877. {https: / /doi.org/10.1007/978-1-4614-8480-6 


Zank GP (2015) Faltering Steps Into the Galaxy: The Boundary Regions of 


the Heliosphere. Ann Rev Astron Astrophys 53:449-500. |https://doi.org/ 


a] 


Low energy cosmic rays 9 


10.1146 /annurev-astro-082214- 122254 

Zhao D, Galazutdinov GA, Linnartz H, Krelowski J (2015) Detection of OH+ 
in Translucent Interstellar Clouds: New Electronic Transitions and Probing 
the Primary Cosmic Ray Ionization Rate. Astrophys J Lett 805(2):L12. 
https: //doi.org/10.1088/2041-8205/805 /2/L12 

Zhao LL, Qin G, Zhang M, Heber B (2014) Modulation of galactic cosmic 
rays during the unusual solar minimum between cycles 23 and 24. Journal 
of Geophysical Research (Space Physics) 119(3):1493-1506. [h 


10.1002/2013JA019550) |arXiv:1310.7076 [astro-ph.SR] 


Zhong L, Wu B, Zheng S, Gu Q (2021) a ai rea of electron-impact ioniza- 
tion cross sections of molecules eee of H, C, N, O, and F. Physics of 


Plasmas 28(8):083505. https: / /doi.org/10.1063 /5.0054342 


Zucker C, Goodman AA, Alves J, et al. (2022) Star formation near the Sun is 
driven by expansion of the Local Bubble. Nature 601(7893):334-337. 


//doi.org/10.1038/s41586-021-04286-5, |jarXiv:2201.05124) [astro-ph.GA] 


Zweibel EG (2013) The microphysics and wee of cosmic rays. Physics 


of Plasmas 20(5):055501. https: //doi.org/10.1063/1.4807033 


Zweibel EG (2017) The basis for cosmic ray feedback: Written on the wind. 


Physics of Plasmas 24(5):055402. https: //doi.org/10.1063/1.4984017 


Zweibel EG, Shull JM (1982) Confinement of cosmic rays in molecular clouds. 


Actvopliys J 259:859--868. https: //doi.org/10.1086/160220 


